Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveforms”

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 667 records · Page 37

Major disruption of D" beneath Alaska

D'' represents one of the most dramatic thermal and compositional layers within our planet. In particular, global tomographic models display relatively fast patches at the base of the mantle along the circum-Pacific which are generally attributed to slab debris. Such distinct patches interact with the bridgmanite (Br) to post-bridgmanite (PBr) phase boundary to generate particularly strong heterogeneity at their edges. Most seismic observations for the D'' come from the lower mantle S wave triplication (Scd). Here, in this work, we exploit the USArray waveform data to examine one of these sharp transitions in structure beneath Alaska. From west to east beneath Alaska, we observed three different characteristics in D'': (1) the western region with a strong Scd, requiring a sharp δVs = 2.5% increase; (2) the middle region with no clear Scd phases, indicating a lack of D'' (or thin Br-PBr layer); and (3) the eastern region with strong Scd phase, requiring a gradient increase in δVs. To explain such strong lateral variation in the velocity structure, chemical variations must be involved. We suggest that the western region represents relatively normal mantle. In contrast, the eastern region is influenced by a relic slab that has subducted down to the lowermost mantle. In the middle region, we infer an upwelling structure that disrupts the Br-PBr phase boundary. Such an interpretation is based upon a distinct pattern of travel time delays, waveform distortions, and amplitude patterns that reveal a circular-shaped anomaly about 5° across which can be modeled synthetically as a plume-like structure rising about 400 km high with a shear velocity reduction of ~5%, similar to geodynamic modeling predictions of upwellings.

58 GEOSCIENCES↗

Characterization of induced seismicity patterns derived from internal structure in event clusters

Microseismicity induced by CO 2 injection at Decatur, Illinois, occurs in distinct clusters and shows no obvious correlation with the proceeding pressure front. Here, we analyze some of these clusters in more depth by using a waveform cross-correlation approach. With this approach we can associate about 1400 events from two clusters, with moment magnitudes between 1.1 and -1.7, with specific formations of much smaller vertical dimensions (tens of meters) than the depth resolution of traveltime-based event locations. The differentiation of reservoir and basement events, and the definition of subclusters by waveform correlation, rather than by location, helps to better analyze the spatiotemporal evolution of the events within a cluster. In the Decatur case, this is characterized by event migration from the reservoir into the adjacent basement. The spatial variation of Brune stress drop and Gutenberg b value exhibits signs of a fluid-driven triggering mechanism at the cluster level, revealing a punctual hydraulic connection between reservoir and basement, most likely associated with basement faults cutting into the reservoir. The observed clustering of microseismicity can thus be explained by the lateral heterogeneity of permeability and crustal strength and is overall consistent with a pressure-induced triggering mechanism. Hence, proper long-term risk mitigation for large-scale fluid injection close to the basement requires prior mapping of small subseismic basement-connected faults.

58 GEOSCIENCES↗

Simulations of Seismic Wave Propagation on Mars

In this paper, we present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust. For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE against the 2D axisymmetric wave propagation solver AxiSEM at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. Finally, we conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources.

58 GEOSCIENCES↗

Inversion of Meteor Rayleigh Waves on Earth and Modeling of Air Coupled Rayleigh Waves on Mars

Meteor impacts and/or meteor events generate body and surface seismic waves on the surface of a planet. When meteoroids burst in the atmosphere, they generate shock waves that subsequently convert into acoustic waves in the atmosphere and seismic waves in the ground. This effect can be modeled as the amplitude of Rayleigh and other Spheroidal modes excitation, due to atmospheric/ground coupling effects. First, an inversion of the seismic source of Chelyabinsk superbolide is performed. We develop an approach in order to model a line source in the atmosphere, corresponding to the consecutive generation of shock waves by the interaction with the atmosphere. The model is based on the known trajectory. We calculate the synthetic seismograms of Rayleigh waves associated with the event by the summation of normal modes of a model of the solid part and the atmosphere of the planet. Through an inversion technique based on singular value decomposition, we perform a full Rayleigh wave inversion and we provide solutions for the moment magnitude. SEIS will likely detect seismic waves generated by impacts and the later might be further located by remote sensing differential processing. In the case of Mars, we use the same method to obtain waveforms associated with impacts on the planetary surface or in low altitudes in the Martian atmosphere. We show that the contribution of the fundamental spheroidal solid mode is dominating the waveforms, compared to that of the first two overtones. We perform an amplitude comparison and we show that small impactors (diameter of 0.5 to 2 m), can be detected by the SEIS VBB seismometer of InSight mission, even in short epicentral distances, in the higher frequencies of the Rayleigh waves. We perform an analysis based on impact rate estimations and we calculate the number of detectable events of 1 meter diameter meteor impacts to be 6.7 to 13.4 per 1 Martian year for a Q=500 .

79 ASTRONOMY AND ASTROPHYSICS↗

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗

Simulations of Events for the LUX-ZEPLIN (LZ) Dark Matter Experiment

The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1–2) ×10−12 pb at a WIMP mass of 40 GeV/ c 2 . This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rolling hills on the core-mantle boundary

Recent results suggest that an iron-rich oxide may have fractionally crystallized from a primordial magma ocean and settled on the core–mantle boundary (CMB). Based on experimental results, the presence of only a few percent of Fe-rich oxide could slow seismic waves down by several percent. This heavy layer can become highly undulating as predicted from dynamic modeling but can remain as a distinct structure with uniform velocity reductions. Here, we use the large USArray seismic network to search for such structures. Strong constraints on D" are provided by the core-phase SKS where it bifurcates, containing a short segment of P-wave diffractions (P d ) when crossing the CMB, called SKS d . Synthetics from models with moderate velocity drops (less than 10%) involving a layer with variable thickness, perhaps a composite of sharp small structures, with strong variation in thickness can explain both the observed SKS d waveforms and large scatter in differential times between SKKS and SKS. A smooth 3D image is obtained from inverting SKS d waveforms displaying rolling-hills with elongated dome-like structures sitting on the CMB. The most prominent one has an 80-km height, ~8° length, and ~4° width, thus adding still more structural complexity to the lower mantle. We suggest that these results can be explained by a dynamically-stabilized material containing small amounts (~5%) iron-rich (Mg,Fe)O providing a self-consistent physical interpretation.

Sun, Daoyuan↗

A digital seismogram archive of nuclear explosion signals, recorded at the Borovoye Geophysical Observatory, Kazakhstan, from 1966 to 1996

Seismologists from Kazakhstan, Russia, and the United States have rescued the Soviet-era archive of nuclear explosion seismograms recorded at Borovoye in northern Kazakhstan during the period 1966–1996. The signals had been stored on about 8000 magnetic tapes, which were held at the recording observatory. After hundreds of man-years of work, these digital waveforms together with significant metadata are now available via the project URL, namely http://www.ldeo.columbia.edu/res/pi/Monitoring/Data/ as a modern open database, of use to diverse communities. Three different sets of recording systems were operated at Borovoye, each using several different seismometers and different gain levels. For some explosions, more than twenty different channels of data are available. A first data release, in 2001, contained numerous glitches and lacked many instrument responses, but could still be used for measuring accurate arrival times and for comparison of the strengths of different types of seismic waves. The project URL also links to our second major data release, for nuclear explosions in Eurasia recorded in Borovoye, in which the data have been deglitched, all instrument responses have been included, and recording systems are described in detail. This second dataset consists of more than 3700 waveforms (digital seismograms) from almost 500 nuclear explosions in Eurasia, many of them recorded at regional distances. It is important as a training set for the development and evaluation of seismological methods of discriminating between earthquakes and underground explosions, and can be used for assessment of three-dimensional models of the Earth’s interior structure.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Underground hydrogen storage leakage detection and characterization based on machine learning of sparse seismic data

Underground hydrogen storage (UHS) is considered as a scalable approach for massive storage and seasonal extraction of hydrogen (H 2 ). Although conventional leakage detection and characterization methods based on time-lapse seismic imaging and inversion generally apply to H 2 leakage detection problem, a high-fidelity yet cost effective geophysics approach is still missing to reliably inform leakage location and properties based on very sparse data. In response, we develop a novel supervised machine learning method to detect and characterize H 2 leakage from UHS. The input to our neural network are sparse time-lapse seismic waveforms, while the output from the neural network includes the spatial location and physical properties of a H 2 leakage. Here, we generate high-quality time-lapse waveforms using the elastic-wave equations to train the neural network. We train and validate our machine learning model and find that it attains high accuracy in using extremely sparse time-lapse seismic data to detect and characterize H 2 leakage. Our investigation is the first systematic study that focuses on applying machine learning to subsurface H 2 leakage detection and characterization and could potentially serve as a cost-effective geophysical tool for underground hydrogen leakage detection and characterization with high fidelity.

08 HYDROGEN↗

In situ monitoring of microstructure evolution during thermal processing of uranium-zirconium alloys using laser-generated ultrasound

In this study, laser-generated ultrasound was used to monitor microstructure evolution during thermal processing of as-cast, polycrystalline binary uranium-zirconium metallic fuel alloys with compositions of U-20wt.%Zr (U-20Zr), U-50wt.%Zr (U-50Zr) and U-80wt.%Zr (U-80Zr). Ultrasonic waveforms were recorded during heating and cooling the samples from room temperature to >973 K and back. A phase transition temperature for all three compositions was estimated from the temperature at which an abrupt and rapid reduction in ultrasonic velocities was observed. Microstructural features on the length scale of tens of micrometers were inferred from the observation of scattering of ultrasonic waves by elastic heterogeneities above ~823 K in U-20Zr, while a hysteresis in the ultrasonic velocities of U-80Zr upon cooling was attributed to a partial retention of the high temperature phase following thermal annealing. The U-50Zr alloy exhibited a reversible viscoelastic response above 933 K, as evidenced by the observation of high frequency attenuation of the shear component of the waveforms at high temperature. Ultrasonic measurements were supplemented by in situ transmission electron microscopy (TEM). The TEM images revealed that the δ-U-Zr matrix in the three compositions underwent a spinodal decomposition above ~823 K into nanoscale regions. The ultrasonic measurements revealed larger, micron-scale structure evolution in the U-20Zr alloy at the same temperature. This large-scale structure is associated with heterogeneous regions having different Zr content. Our findings show the potential heating rate dependence of microstructural evolution in U-Zr alloys and highlight differences in the thermomechanical response and associated length scales during thermal annealing between single- and dual-phase compositions. These results demonstrate the utility of laser ultrasonics to rapidly and efficiently scan phase boundaries and monitor micrometer-scale structure evolution in metallic fuel alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Tectonic tremor characterized by principal-component analysis in the vicinity of central Chile and Argentina

No conclusive evidence has been presented to date for tectonic tremor (TT) in the vicinity of central Chile, where the Nazca Plate is subducting beneath the South American Plate. Subduction in our experimental location (roughly 35.5° S, 70.5° W) is steep and fairly unobstructed compared to the flattened and more seismogenic behavior to the north. We seek to identify TT in our experimental area, whose geodynamics are comparable to tremor-rich subduction zones such as Cascadia and the Nankai Trough. Our method combines time-series visual inspection, frequency-spectrum analysis, waveform cross-correlation, and 3-component (3C) signal covariance to explore the presence of TT in this region. We have identified TT using stations in central Chile and the Malargüe region, Argentina. The TT exhibits similar features to other TT observations worldwide. Waveform characteristics for the TT in our study, particularly dimension of the 3C signal covariance, vary as a function of apparent source location. The duration of one episode of identified TT was about 10 h, which may indicate that the plate interface where tremor generates is strongly coupled. Here, we conclude that our observations reflect features of the local propagation, rather than the tremor source itself.

58 GEOSCIENCES↗

Improved neutron-gamma discrimination for a 3 He neutron detector using subspace learning methods

3 He gas based neutron linear-position-sensitive detectors (LPSDs) have been applied for many neutron scattering instruments. Traditional Pulse-Height Analysis (PHA) for Neutron-Gamma Discrimination (NGD) resulted in the neutron-gamma efficiency ratio on the orders of 10 5 -10 6 . The NGD ratios of 3 He detectors need to be improved for even better scientific results from neutron scattering. Digital Signal Processing (DSP) analyses of waveforms were proposed for obtaining better NGD ratios, based on features extracted from rise-time, pulse amplitude, charge integration, a simplified Wiener filter, and the cross-correlation between individual and template waveforms of neutron and gamma events. Fisher linear discriminant analysis (FLDA) and three multivariate analyses (MVAs) of the features were performed. The NGD ratios are improved by about 10 2 -10 3 times compared with the traditional PHA method. Finally, our results indicate the NGD capabilities of 3 He tube detectors can be significantly improved with subspace-learning based methods, which may result in a reduced data-collection time and better data quality for further data reduction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Belle II imaging Time-of-Propagation (iTOP) Detector

High precision flavor physics measurements are an essential complement to the direct searches for new physics at the LHC. Such measurements will be performed using the upgraded Belle II detector that will take data at the SuperKEKB accelerator. With 40x the luminosity of KEKB, the detector systems must operate eciently at much higher rates than the original Belle detector. A central element of the upgrade is the barrel particle identication system. Belle II has built and installed an imaging-Time-of-Propagation (iTOP) detector. The iTOP uses quartz optics as Cerenkov radiators. The photons are transported down the quart bars via total internal reection with a spherical mirror at the forwardend to reflect photons to the backward end where they are imaged onto an array of segmented Micro-Channel Plate Photo-Multiplier Tubes. The system is readout using gigsample per second waveform sampling Application-Specific Integrated Circuits that provide precise photon timing. The combined timing and spatial distribution of the photons for each event are used to determine particle species. This paper provides an overview of the iTOP system.

Fast, James E.↗

Calibration of a compact ASIC-based data acquisition system for neutron/$γ$ discrimination and spectroscopy with organic scintillators

Segmented neutron detectors that use silicon photomultipliers (SiPMs) are receiving significant attention in nuclear security applications. Some of these detectors employ hundreds of channels and would therefore benefit from the use of high-channel-density and low-cost-per-channel data acquisition (DAQ) systems. Candidate DAQ systems that meet these requirements exist, but few perform full waveform digitization, which permits neutron and gamma-ray interaction discrimination via pulse shape. In this work, we study the performance of the TOFPET2 (PETsys Electronics), an ASIC-based DAQ designed for positron emission tomography, which has been adapted to provide sensitivity to pulse shape by the use of variable-period charge integration. We use a light-emitting diode to calibrate a combination of an ON Semiconductor (SensL) 60035-64P J-Series SiPM and TOPFET2 DAQ and evaluate the linearity of their response and dynamic range. Here, the calibration curve was obtained by comparing the DAQ response to that measured with a photodiode interfaced with a traditional waveform digitizer. This calibration process was used for rudimentary spectroscopy of various neutron and sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The data acquisition system of the LZ dark matter detector: FADR

The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. Furthermore, the hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗