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 739 records · Page 41

Gravitational wave tests of quantum modifications to black hole structure—with post-GW150914 update

Here a preliminary discussion is given of the prospects that gravitational-wave observations of binary inspiral of black holes (BHs) could reveal or constrain quantum modifications to BH dynamics, such as are required to preserve postulates of quantum mechanics. Different proposals for such modifications are characterized by different scales, and the size of these scales relative to those probed by observation of inspiral signals is important in determining the feasibility of finding experimental signatures. Certain scenarios with strong quantum modifications in a region extending well outside the horizon are expected to modify classical evolution, and distort the near-peak gravitational wave signal, suggesting a search for departures from waveforms predicted by general relativity. The near agreement of the GW150914 signal with such waveforms is discussed, and indicates constraints on some such scenarios. Important strategies for more precise future tests are (1) to develop more precise predictions from scenarios proposing quantum modifications, and (2) searching for observed deviations from numerical relativity predictions via analysis of gravity wave data, particularly focussing on the signal region corresponding to plunge and merger.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fluctuation statistics in the scrape-off layer of Alcator C-Mod

We study long time series of the ion saturation current and floating potential, sampled by Langmuir probes dwelled in the outboard mid-plane scrape-off layer and embedded in the lower divertor baffle of Alcator C-Mod. A series of ohmically heated L-mode plasma discharges is investigated with line-averaged plasma density ranging from ${{\bar{n}}_{\text{e}}}/{{n}_{\text{G}}}=0.15$ to 0.42, where $n_ G$ is the Greenwald density. All ion saturation current time series that are sampled in the far scrape-off layer are characterized by large-amplitude burst events. Coefficients of skewness and excess kurtosis of the time series obey a quadratic relationship and their histograms coincide partially upon proper normalization. Histograms of the ion saturation current time series are found to agree well with a prediction of a stochastic model for the particle density fluctuations in scrape-off layer plasmas. The distribution of the waiting times between successive large-amplitude burst events and of the burst amplitudes are approximately described by exponential distributions. The average waiting time and burst amplitude are found to vary weakly with the line-averaged plasma density. Conditional averaging reveals that the radial blob velocity, estimated from floating potential measurements, increases with the normalized burst amplitude in the outboard mid-plane scrape-off layer. For low density discharges, the conditionally averaged waveform of the floating potential associated with large amplitude bursts at the divertor probes has a dipolar shape. In detached divertor conditions the average waveform is random, indicating electrical disconnection of blobs from the sheaths at the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sequence modeling of higher-order wave modes of quasi-circular, spinning, non-precessing binary black hole mergers

Higher-order gravitational wave modes from quasi-circular, spinning, non-precessing binary-black-hole (BBH) mergers encode rich information about the nonlinear dynamics of strong-field gravity. We present a transformer-based sequence-completion surrogate that, given an early-inspiral segment, forecasts the subsequent late inspiral, merger, and ringdown. The intended applications are (i) patching or completing expensive or interrupted numerical-relativity (NR) simulations and (ii) providing late-time cross-checks and rapid hybridization studies. The training set is built from the NRHybSur3dq8 surrogate, which provides spherical-harmonic modes up to $\ell$ ≤ 4 (excluding (4, 0) and (4,±1), and including (5, 5)) for mass ratios q ≤ 8, dimensionless spin components s$^{z}_{1,2}$ ϵ[–0.8, 0.8], and inclination angles θ ϵ [0, π]. Waveforms are supplied on the interval t ϵ [–5000M, –100 M) and the model autoregressively generates the plus and cross polarizations (h + , h x ) on t ϵ [–100 M, 130M]. Training on the Delta supercomputer with 16 NVIDIA A100 GPUs required ~15 h on more than 14 million hybrid waveforms. Evaluation on a held-out test set of 840,000 samples yields mean and median overlaps of 0.996 and 0.997, respectively, with respect to the surrogate ground truth.

black-hole merger↗

Pulsed photoemission induced plasma breakdown

Abstract This article characterises the effects of cathode photoemission leading to electrical discharges in an argon gas. We perform breakdown experiments under pulsed laser illumination of a flat cathode and observe Townsend to glow discharge transitions. The breakdown process is recorded by high-speed imaging, and time-dependent voltage and current across the electrode gap are measured for different reduced electric fields and laser intensities. We employ a 0D transient discharge model to interpret the experimental measurements. The fitted values of transferred photoelectron charge are compared with calculations from a quantum model of photoemission. The breakdown voltage is found to be lower with photoemission than without. When the applied voltage is insufficient for ion-induced secondary electron emission to sustain the plasma, laser driven photoemission can still create a breakdown where a sheath (i.e. a region near the electrode surfaces consisting of positive ions and neutrals) is formed. This photoemission induced plasma persists and decays on a much longer time scale ( ∼ 10 s μ s) than the laser pulse length ( 30 ps). The effects of different applied voltages and laser energies on the breakdown voltage and current waveforms are investigated. The discharge model can accurately predict the measured breakdown voltage curves, despite the existence of discrepancy in quantitatively describing the transient discharge current and voltage waveforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ammonia generation in Ns pulse and Ns pulse/RF discharges over a catalytic surface

Plasma-catalytic ammonia synthesis in a ns pulse discharge and a 'hybrid' ns pulse/RF discharge in plane-to-plane geometry is studied by Fourier Transform infrared absorption spectroscopy. The data are taken in a preheated H 2 –N 2 mixture, with and without Ni/γ-Al 2 O 3 or Co/γ-Al 2 O 3 catalyst placed in the discharge section. The measurement results are taken using two different approaches. The first is a 'single-stage' process, where a ns pulse discharge in the H 2 –N 2 mixture is sustained continuously. In this case, the ammonia yield increases slowly, over a period of tens of minutes. The second is a 'two-stage' process, where the catalyst is first activated by the ns pulse discharge sustained in pure nitrogen, and then the activated catalyst is exposed to the H 2 –N 2 flow, with or without the discharge. In this case, a strong overshoot of the NH 3 number density at the reactor exit is detected, by over a factor of two compared to the single-stage process. After the initial overshoot, the ammonia yield gradually decreases to the 'single stage' value (with the discharge on), or to near zero (with the discharge off). The results demonstrate that the ammonia yield in the plasma-catalytic reactor is controlled by the N atom accumulation on the catalyst surface, which reacts with H atoms thermally dissociated on the catalyst or generated in the plasma. The results also show that the plasma-catalytic ammonia yield is significantly higher compared to that in the ns pulse discharge without the catalyst. The accumulation of H atoms on the catalyst, with their subsequent reactions with N atoms generated in the plasma, is of relatively minor importance at the present conditions. An additional series of measurements was made with a sub-breakdown RF waveform overlapped with the ns pulse discharge train, to enhance the vibrational excitation of nitrogen. The ammonia yield measured with the RF waveform added is approximately 20% higher compared to that at the baseline ns pulse discharge conditions, both with and without the catalyst. This effect is weaker compared to that of the catalyst activation by N atoms. Additional data are necessary to isolate the possible effect of the vibrationally excited N 2 molecules on the ammonia synthesis in the plasma catalytic reactions.

surface reactions↗

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

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

97 MATHEMATICS AND COMPUTING↗

The electronics and data acquisition system for the DarkSide-50 veto detectors

DarkSide-50 is a detector for dark matter candidates in the form of weakly interacting massive particles. It utilizes a liquid argon time projection chamber for the inner main detector, surrounded by a liquid scintillator veto (LSV) and a water Cherenkov veto detector (WCV). The LSV and WCV act as the neutron and cosmogenic muon veto detectors for DarkSide-50. This paper describes the electronics and data acquisition system used for these two detectors. The system is made of a custom built front end electronics and commercial National Instruments high speed digitizers. The front end electronics, the DAQ, and the trigger system have been used to acquire data in the form of zero-suppressed waveform samples from the 110 PMTs of the LSV and the 80 PMTs of the WCV. The veto DAQ system has proven its performance and reliability. This electronics and DAQ system can be scaled and used as it is for the veto of the next generation DarkSide-20k detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The electronics, trigger and data acquisition system for the liquid argon time projection chamber of the DarkSide-50 search for dark matter

The DarkSide-50 experiment at the Laboratori Nazionali del Gran Sasso is a search for dark matter using a dual phase time projection chamber with 50 kg of low radioactivity argon as target. Light signals from interactions in the argon are detected by a system of 38 photo-multiplier tubes (PMTs), 19 above and 19 below the TPC volume inside the argon cryostat. We describe the electronics which processes the signals from the photo-multipliers, the trigger system which identifies events of interest, and the data-acquisition system which records the data for further analysis. The electronics include resistive voltage dividers on the PMTs, custom pre-amplifiers mounted directly on the PMT voltage dividers in the liquid argon, and custom amplifier/discriminators (at room temperature). After amplification, the PMT signals are digitized in CAEN waveform digitizers, and CAEN logic modules are used to construct the trigger, the data acquisition system for the TPC is based on the Fermilab "artdaq" software. The system has been in operation since early 2014.

Agnes, P.↗

Deep neural networks for energy and position reconstruction in EXO-200

Here, we apply deep neural networks (DNN) to data from the EXO-200 experiment. In the studied cases, the DNN is able to reconstruct the relevant parameters—total energy and position—directly from raw digitized waveforms, with minimal exceptions. For the first time, the developed algorithms are evaluated on real detector calibration data. The accuracy of reconstruction either reaches or exceeds what was achieved by the conventional approaches developed by EXO-200 over the course of the experiment. Most existing DNN approaches to event reconstruction and classification in particle physics are trained on Monte Carlo simulated events. Such algorithms are inherently limited by the accuracy of the simulation. We describe a unique approach that, in an experiment such as EXO-200, allows to successfully perform certain reconstruction and analysis tasks by training the network on waveforms from experimental data, either reducing or eliminating the reliance on the Monte Carlo.

97 MATHEMATICS AND COMPUTING↗

Multi-channel front-end ASIC for a 3D position-sensitive detector

Arrays of 3D position-sensitive detectors (3DPSD), operating at room temperature and using cadmium zinc telluride (CZT) and thallium bromide (TIBr) sensors, are suitable for gamma-ray spectrometry in many applications. One detector configuration, the 3D position-sensitive Virtual Frisch-Grid detector (VFG), is particularly advantageous for integrating into large area arrays. The signals generated inside each detector of the array are captured with the anode, cathode and four pads that enable the reconstruction of the position and energy of the ionizing interaction by measurements of amplitude and timing of the signals. For these applications, a low-noise front-end ASIC has been developed, capable of processing bipolar signals (needed because of AC-coupling of certain electrodes). The ASIC can be coupled to an ADC in order to form a compound “waveform digitizer” capable of post-processing the analog signals and determining amplitude and timing information. This paper describes a 32-channel front-end ASIC that is suitable for reading out a 3 × 3 or 4 × 4 element matrix in the VFG configuration. Each channel is composed of a low-noise charge amplifier with an adaptive continuous reset feedback circuit suitable for both positive and negative charge, a first order shaper and a single-to-differential converter output stage. Voltage and current references are all internally generated by 10-bit DACs and the chip is fully controllable with the I 2 C communication protocol. The readout channel response has been verified using the implemented injection circuit. Linear behavior up to ~75 ke ± with the gain of ~80 mV/fC, and up to ~200 ke ± with the gain of ~30 mV/fC was demonstrated. In conclusion, the first test result waveforms using a 137 Cs radioactive source on a 5 × 5 × 12 mm 3 TIBr crystal are reported.

47 OTHER INSTRUMENTATION↗

GRAPH — an readout ASIC for large MCP based detectors

We present a programmable 16 channel, mixed signal, low power readout ASIC, having the project historically named Gigasample Recorder of Analog waveforms from a PHotodetector (GRAPH). It is designed to read large aperture single photon imaging detectors using micro channel plates for charge multiplication, and measuring the detector's response on crossed strips anodes to extrapolate the incoming photon position. Each channel consists of a fast, low power and low noise charge sensitive amplifier, which provides a myriad of coarse and fine programmable options for gain and shaping settings. Further, the amplified signal is recorded using, to our knowledge novel, the Hybrid Universal sampLing Architecture (HULA) ADC. A kind of mixed signal double buffer memory, that enables concurrent waveform recording, and selected event digitized data extraction. The sampling frequency is freely adjustable between few kHz up to 125 MHz, while the chip's internal digital memory holds a history 2048 samples for each channel, with a digital headroom of 12 bits. An optimized region of interest sample-read algorithm allows to extract the information just around the event pulse peak, while selecting the next event, thus substantially reducing the operational dead time. The chip is designed in 130 nm TSMC CMOS technology, and its power consumption is around 47 mW per channel.

47 OTHER INSTRUMENTATION↗

Discrete solitons and vortices in anisotropic hexagonal and honeycomb lattices

We consider the self-focusing discrete nonlinear Schrödinger equation on hexagonal and honeycomb lattice geometries. Our emphasis is on the study of the effects of anisotropy, motivated by the tunability afforded in recent optical and atomic physics experiments. We find that multi-soliton and discrete vortex states undergo destabilizing bifurcations as the relevant anisotropy control parameter is varied. Furthermore, we quantify these bifurcations by means of explicit analytical calculations of the solutions, as well as of their spectral linearization eigenvalues. Finally, we corroborate the relevant stability picture through direct numerical computations. In the latter, we observe the prototypical manifestation of these instabilities to be the spontaneous rearrangement of the solution, for larger values of the coupling, into localized waveforms typically centered over fewer sites than the original unstable structure. In weak coupling, the instability appears to result in a robust breathing of the relevant waveforms.

97 MATHEMATICS AND COMPUTING↗

A HYBRID MODE MODEL OF THE BLAZHKO EFFECT, SHOWN TO ACCURATELY FIT KEPLER DATA FOR RR Lyr

The waveform for Blazhko stars can be substantially different during the ascending and descending parts of the Blazhko cycle. A hybrid model, consisting of two component oscillators of the same frequency, is proposed as a means to fit the data over the entire cycle. One component exhibits a sawtooth-like velocity waveform while the other is nearly sinusoidal. One method of generating such a hybrid is presented: a nonlinear model is developed for the first overtone mode, which, if excited to large amplitude, is found to drop strongly in frequency and become highly non-sinusoidal. If the frequency drops sufficiently to become equal to the fundamental frequency, the two can become phase locked and form the desired hybrid. A relationship is assumed between the hybrid mode velocity and the observed light curve, which is approximated as a power series. An accurate fit of the hybrid model is made to actual Kepler data for RR Lyr. The sinusoidal component may tend to stabilize the period of the hybrid which is found in real Blazhko data to be extremely stable. It is proposed that the variations in amplitude and phase might result from a nonlinear interaction with a third mode, possibly a nonradial mode at 3/2 the fundamental frequency. The hybrid model also applies to non-Blazhko RRab stars and provides an explanation for the light curve bump. A method to estimate the surface gravity is also proposed.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes ($\ell$, |m|) = {(2,2), (2,1), (3,3), (3,2), (4,4)}, and mode mixing effects in the $\ell$ = 3, |m| = 2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

79 ASTRONOMY AND ASTROPHYSICS↗

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

Modelling infrasound signal generation from two underground explosions at the Source Physics Experiment using the Rayleigh integral

For this study, we use the Rayleigh integral (RI) as an approximation to the Helmholtz–Kirchoff integral to model infrasound generation and propagation from underground chemical explosions at distances of 250 m out to 5 km as part of the Source Physics Experiment (SPE). Using a sparse network of surface accelerometers installed above ground zero, we are able to accurately create synthetic acoustic waveforms and compare them to the observed data. Although the underground explosive sources were designed to be symmetric, the resulting seismic wave at the surface shows an asymmetric propagation pattern that is stronger to the northeast of the borehole. This asymmetric bias may be attributed to the subsurface geology and faulting of the area and is observed in the acoustic waveforms. We compare observed and modelled results from two of the underground SPE tests with a sensitivity study to evaluate the asymmetry observed in the data. This work shows that it is possible to model infrasound signals from underground explosive sources using the RI and that asymmetries observed in the data can be modelled with this technique.

58 GEOSCIENCES↗

An Adaptable Seismic Data Format

We present ASDF, the Adaptable Seismic Data Format, a modern and practical data format for all branches of seismology and beyond. The growing volume of freely available data coupled with ever expanding computational power opens avenues to tackle larger and more complex problems. Current bottlenecks include inefficient resource usage and insufficient data organization. Properly scaling a problem requires the resolution of both these challenges, and existing data formats are no longer up to the task. ASDF stores any number of synthetic, processed or unaltered waveforms in a single file. A key improvement compared to existing formats is the inclusion of comprehensive meta information, such as event or station information, in the same file. Additionally, it is also usable for any non-waveform data, for example, cross-correlations, adjoint sources or receiver functions. Last but not least, full provenance information can be stored alongside each item of data, thereby enhancing reproducibility and accountability. Any data set in our proposed format is self-describing and can be readily exchanged with others, facilitating collaboration. The utilization of the HDF5 container format grants efficient and parallel I/O operations, integrated compression algorithms and check sums to guard against data corruption. To not reinvent the wheel and to build upon past developments, we use existing standards like QuakeML, StationXML, W3C PROV and HDF5 wherever feasible. Usability and tool support are crucial for any new format to gain acceptance. We developed mature C/Fortran and Python based APIs coupling ASDF to the widely used SPECFEM3D_GLOBE and ObsPy toolkits.

58 GEOSCIENCES↗