Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave arrival time”

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 109 records · Page 6

Optical fiber sensors for damage analysis in aerospace materials

Under this grant, fiber optic sensors were investigated for use in the nondestructive evaluation of aging aircraft. Specifically, optical fiber sensors for detection and location of impacts on a surface, and for detection of corrosion in metals were developed. The use of neural networks was investigated for determining impact location by processing the output of a network of fiberoptic strain sensors distributed on a surface. This approach employs triangulation to determine location by comparing the arrival times at several sensors, of the acoustic signal generated by the impact. For this study, a neural network simulator running on a personal computer was used to train a network using a back-propagation algorithm. Fiber optic extrinsic Fabry-Perot interferometer (EFPI) strain sensors are attached to or embedded in the surface, so that stress waves emanating from an impact can be detected. The ability of the network to determine impact location by time-or-arrival of acoustic signals was assessed by comparing network outputs with actual experimental results using impacts on a panel instrumented with optical fiber sensors. Using the neural network to process the sensor outputs, the impact location can be inferred to centimeter range accuracy directly from the arrival time data. In addition, the network can be trained to determine impact location, regardless of material anisotropy. Results demonstrate that a back-propagation network identifies impact location for an anisotropic graphite/bismaleimide plate with the same accuracy as that for an isotropic aluminum plate. Two different approaches were investigated for the development of fiber optic sensors for corrosion detection in metals, both utilizing optical fiber sensors with metal coatings. In the first approach, an extrinsic Fabry-Perot interferometric fiber optic strain sensor was placed under tensile stress, and while in the resulting strained position, a thick coating of metal was applied. Due to an increase in the quantity of material, the sensor does not return to its original position upon removal of the applied stress, and some residual strain is maintained within the sensor element. As the metal thickness decreases due to corrosion, this strain is released, providing the sensing mechanism for corrosion detection. In the second approach, photosensitive optical fibers with long period Bragg gratings in the core were coated with metal. The Bragg gratings serve to couple core modes at discrete wavelengths to cladding modes. Since cladding modes interact with the metal coating surrounding the fiber cladding, the specific wavelengths coupled from core to cladding depend on the refractive index of the metal coating. Therefore, as the metal corrodes, the resulting change in index of the coating may be measured by measuring the change in wavelength of the coupled mode. Results demonstrate that both approaches can be successfully used to track the loss in metal coating on the optical fiber sensors due to corrosion.

Schindler, Paul↗

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗

Modified gravitational wave propagation with higher modes and its degeneracies with lensing

Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism — the WKB approach — for solving analytically wave propagation in the spatial domain with a modified dispersion relation (MDR). We compare this WKB approach to applying a stationary phase approximation (SPA) in the temporal domain with time delays associated to the group or particle velocity. To this end, we extend the SPA to generic signals with higher modes, keeping careful track of reference phases and arrival times. We find that the WKB approach coincides with the SPA using the group velocity, in agreement with the principles of wave propagation. We then explore the degeneracies between a GW propagation with an MDR and a strongly-lensed GW in GR, since the latter can introduce a frequency-independent phase shift which is not degenerate with source parameters in the presence of higher modes. We find that for a particular MDR there is an exact degeneracy for wave propagation, unlike with the SPA for particle propagation. For the other cases, we search for the values of the MDR parameters that minimize the χ 2 and conclude that strongly-lensed GR GWs could be misinterpreted as GWs in modified gravity. As a result, future MDR constraints with higher mode GWs should include the possibility of frequency-independent phase shifts, allowing for the identification of modified gravity and strong lensing distortions at the same time.

79 ASTRONOMY AND ASTROPHYSICS↗

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗

On the determination of shocked surface velocities

Laser Doppler velocimetry is commonly fielded in shock physics experiments to diagnose the times of arrival, magnitudes, and durations of shock impulses as they arrive at shock loaded surfaces. However, in the presence of ejecta, the resulting frequency noise often obscures the surface and subsequent impulses at later times. In this context, we evaluate a dominant wave analysis technique applied to velocimetry data from explosively shocked Sn subjected to a second shockwave, on which ejecta are emitted subsequent to both the first and subsequent shockwave impulses. The analysis reveals surface velocities, the times of arrival of subsequent impulses, and their durations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improved Ultrasonic Transducer For Measuring Cryogenic Flow

Improved ultrasonic transducer used to measure flow of cryogenic fluid. Includes wedge made nonintrusive by machining it out of bulk material of duct carrying fluid. Skewed surfaces of wedge suppress standing waves, thus reducing ringing and increasing signal-to-noise ratio. Increases accuracy of measurements of times of arrival of ultrasonic pulses, from which times flow inferred.

Barkhoudarian, Sarkis↗

A statistical study of Pc 3-5 pulsations observed by the AMPTE/CCE magnetic fields experiment. I - Occurrence distributions

A classification system based on dynamic power spectra of AMPTE/CCE magnetic field data has been developed which allows semiquantitative specification of pulsation activity at all times. Fifteen months of 6.24-s AMPTE/CCE magnetic field data have been analyzed using this scheme to generate a comprehensive data base of all pulsation activity observed from L = 5 to L = 9 at all local times in the frequency range 0-80 mHz. The pulsation divide naturally into five basic categories: toroidal fundamental mode resonances, toroidal harmonic resonances, compressional low-frequency waves, radially polarized waves, and disturbed intervals. Additional distinctions among members of these classes are made to arrive at 14 categories which describe 95 percent of the pulsations observed. Normal spatial distributions in local time and L shell are determined for each category. Magnetic latitude distributions are determined for four coherent categories of activity.

Anderson, B. J.↗

Solar Energetic Particle Variations

In the largest solar energetic-particle (SEP) events, acceleration occurs at shock waves driven out from the Sun by coronal mass ejections (CMEs). In fact, the highest proton intensities directly measured near Earth at energies up to approximately 1 GeV occur at the time of passage of shocks, which arrive about a day after the CMEs leave the Sun. CME-driven shocks expanding across magnetic fields can fill over half of the heliosphere with SEPs. Proton-generated Alfven waves trap particles near the shock for efficient acceleration but also throttle the intensities at Earth to the streaming limit early in the events. At high energies, particles begin to leak from the shock and the spectrum rolls downward to form an energy-spectral 'knee' that can vary in energy from approximately 1 MeV to approximately 1 GeV in different events. All of these factors affect the radiation dose as a function of depth and latitude in the Earth's atmosphere and the risk to astronauts and equipment in space. SEP ionization of the polar atmosphere produces nitrates that precipitate to become trapped in the polar ice. Observations of nitrate deposits in ice cores reveal individual large SEP events and extend back approximately 400 years. Unlike sunspots, SEP events follow the approximately 80-100-year Gleissberg cycle rather faithfully and are now at a minimum in that cycle. The largest SEP event in the last 400 years appears to be related to the flare observed by Carrington in 1859, but the probability of SEP events with such large fluences falls off sharply because of the streaming limit.

Reames, D. V.↗

An End-To-End Earthquake Detection Method for Joint Phase Picking and Association Using Deep Learning

Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. Here, in this study, we introduce a new end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations to associate and detect earthquakes across a seismic network. We use these sub-networks together with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the neural network model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches. Because our approach preserves information across tasks in the detection pipeline, it has the potential to outperform approaches that do not.

58 GEOSCIENCES↗

Wigner time delay in photoionization: a 1D model study

Abstract In scattering theory, the Wigner–Smith time delay, calculated through a phaseshift derivative or its multichannel generalization, has been demonstrated to measure the amount of delay or advance experienced by colliding particles during their interaction with the scattering potential. Fetic, Becker, and Milosevic argue that this concept cannot be extended to include photoionization, viewed as a half-scattering experiment. Their argument is based on the lack of information about scattering phaseshifts in the part of the wavefunction (satisfying the ingoing-wave boundary condition) going to the detector. This article aims to test this claim by examining a photodetachment process in a simple 1D model with a short-range symmetrical potential. Using time-dependent perturbation theory with a dipole interaction, the relevant wavepacket of the outgoing particle is analyzed and compared to the free wavepacket as a reference. Our findings confirm that, indeed, a time delay arises in the liberated fragmentation wavepacket, which is expressed as an energy derivative of the scattering phaseshift. Our study highlights that the source of the phaseshift content in the wavepacket arriving at the detector is the dipole matrix element, which is a direct consequence of imposing the ingoing-wave boundary condition. We illustrate our results through numerical simulations of both the non-free and free wavepackets. The amount of the observed time delay is found to be half of that appearing in a typical scattering experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗

A theoretical interpretation of Forbush decreases

An analysis is presented of the intensity variations of relativistic cosmic rays, at 1 AU, which are associated with the passage of a flare shock wave. The magnitude and time profile obtained is similar to that observed in flare-induced Forbush decreases, and the principal cause of the energy reduction is an increase in adiabatic cooling of the arriving particles due to prolonged containment behind the compressed field of the flare shock wave. The large decreases calculated for the smaller diffusion coefficients are inconsistent with observation, implying larger mean free paths than those traditionally assumed at these rigidities. The method presented may be used to study intensity variations associated with other large-scale structures of the interplanetary field, such as corotating interaction regions.

Gall, R.↗

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↗

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗

Multiplexing Technology for Acoustic Emission Monitoring of Aerospace Vehicles

The initiation and propagation of damage mechanisms such as cracks and delaminations generate acoustic waves, which propagate through a structure. These waves can be detected and analyzed to provide the location and severity of damage as part of a structural health monitoring (SHM) system. This methodology of damage detection is commonly known as acoustic emission (AE) monitoring, and is widely used on a variety of applications on civil structures. AE has been widely considered for SHM of aerospace vehicles. Numerous successful ground and flight test demonstrations have been performed, which show the viability of the technology for damage monitoring in aerospace structures. However, one significant current limitation for application of AE techniques on aerospace vehicles is the large size, mass, and power requirements for the necessary monitoring instrumentation. To address this issue, a prototype multiplexing approach has been developed and demonstrated in this study, which reduces the amount of AE monitoring instrumentation required. Typical time division multiplexing techniques that are commonly used to monitor strain, pressure and temperature sensors are not applicable to AE monitoring because of the asynchronous and widely varying rates of AE signal occurrence. Thus, an event based multiplexing technique was developed. In the initial prototype circuit, inputs from eight sensors in a linear array were multiplexed into two data acquisition channels. The multiplexer rapidly switches, in less than one microsecond, allowing the signals from two sensors to be acquired by a digitizer. The two acquired signals are from the sensors on either side of the trigger sensor. This enables the capture of the first arrival of the waves, which cannot be accomplished with the signal from the trigger sensor. The propagation delay to the slightly more distant neighboring sensors makes this possible. The arrival time from this first arrival provides a more accurate source location determination. The multiplexer also identifies which channels are acquired by encoding TTL logic pulses onto the latter portion of the signals. This prototype system was demonstrated using pencil lead break (Hsu-Neilsen) sources on an aluminum plate. It performed as designed providing rapid low noise trigger based switching with encoded channel identification. this multiplexing approach is not limited to linear arrays, but can be easily extended to monitor sensors in planar ot three dimensional arrays. A 32 channel multiplexing system is under development that will allow arbitrary sensor placement. Another benefit of this multiplexing system is the reduction in the expense of data acquisition hardware. In addition, the reduced weight and power requirements are of extreme importance for proposed AE systems on aerospace vehicles.

Prosser, William↗

Current Closure in the Auroral Ionosphere: Results from the Auroral Current and Electrodynamics Structure Rocket Mission

The Auroral Current and Electrodynamics Structure (ACES) mission consisted of two sounding rockets launched nearly simultaneously from Poker Flat Research Range, AK on January 29, 2009 into a dynamic multiple-arc aurora. The ACES rocket mission was designed to observe electrodynamic and plasma parameters above and within the current closure region of the auroral ionosphere. Two well instrumented payloads were flown along very similar magnetic field footprints, at different altitudes, with small temporal separation between both payloads. The higher altitude payload (apogee 360 km), obtained in-situ measurements of electrodynamic and plasma parameters above the current closure region to determine the input signature. The low altitude payload (apogee 130 km), made similar observations within the current closure region. Results are presented comparing observations of the electric fields, magnetic components, and the differential electron energy flux at magnetic footpoints common to both payloads. In situ data is compared to the ground based all-sky imager data, which presents the evolution of the auroral event as the payloads traversed through magnetically similar regions. Current measurements derived from the magnetometers on the high altitude payload observed upward and downward field-aligned currents. The effect of collisions with the neutral atmosphere is investigated to determine it is a significant mechanism to explain discrepancies in the low energy electron flux. The high altitude payload also observed time-dispersed arrivals in the electron flux and perturbations in the electric and magnetic field components, which are indicative of Alfven waves.

Kaeppler, S. R.↗