Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveform”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Spatial and Temporal Variations in the Moment Tensor Solutions of the 2008 Wenchuan Earthquake Aftershocks and Their Tectonic Implications

Following the mainshock of the 2008 M8 Wenchuan Earthquake, there were more than 300 M L ≥ 4.0 aftershocks that occurred between 12 May 2008 and 8 September 2010. We analyzed here the broadband waveforms for these events and found 160 events with sufficient signal-to-noise levels to invert for seismic moment tensors. Considering the length of the activated fault and the distances to the recording stations, four velocity models were employed to account for variability in crustal structure. The moment tensor solutions show considerable variations with a mixture of mainly reverse and strike-slip mechanisms and a small number of normal events and ambiguous events. We analyzed the spatial and temporal distribution of the aftershocks and their mechanism types to characterize the structure and the deformation occurring in the Longmen Shan fold and thrust belt. Our results suggest that the stress is very complex at the Longmen Shan fault zone. The moment tensors have both a spatial segmentation with two major categories of the moment tensor of thrust and strike slip; and a temporal pattern that the majority of the aftershocks gradually migrated to thrust-type events. The variability of aftershock mechanisms is a strong indication of significant tectonic release and stress reorganization that activated numerous small faults in the system.

58 GEOSCIENCES↗
DOE OSTI2018

Local Infrasound Variability Related to In Situ Atmospheric Observation

Abstract Local infrasound is widely used to constrain source parameters of near‐surface events (e.g., chemical explosions and volcanic eruptions). While atmospheric conditions are critical to infrasound propagation and source parameter inversion, local atmospheric variability is often ignored by assuming homogeneous atmospheres, and their impact on the source inversion uncertainty has never been accounted for due to the lack of quantitative understanding of infrasound variability. We investigate atmospheric impacts on local infrasound propagation by repeated explosion experiments with a dense acoustic network and in situ atmospheric measurement. We perform full 3‐D waveform simulations with local atmospheric data and numerical weather forecast model to quantify atmosphere‐dependent infrasound variability and address the advantage and restriction of local weather data/numerical weather model for sound propagation simulation. Numerical simulations with stochastic atmosphere models also showed nonnegligible influence of atmospheric heterogeneity on infrasound amplitude, suggesting an important role of local turbulence.

58 GEOSCIENCES↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2020

Pervasive low-velocity layer atop the 410-km discontinuity beneath the northwest Pacific subduction zone: Implications for rheology and geodynamics

Regional triplication waveforms of five intermediate-depth events are modeled to simultaneously obtain the compressional (P) and shear (SH) wave velocity structure beneath northwestern Pacific subduction zone. Both the P- and SH-wave velocity models for three different sub-regions show a low-velocity layer (LVL) with a thickness of ~55-80 km lying above the 410-km discontinuity with a ~900 km lateral extent from the Japan Sea to the northeastern Asian continental margin. With the dihedral angle approaching to zero around 400 km, a minute amount of melt atop the 410-km discontinuity caused by the hydrous slab might completely wet olivine grain boundaries and result in a low seismic velocity layer in this specific subduction zone. This mechanism suggests that the 410-LVL is a low viscosity zone that would partially decouple the upper mantle from the transition zone. We infer that the widespread 410-LVL provides evidence for a water-bearing mantle transition zone beneath the western Pacific subduction zone.

58 GEOSCIENCES↗
DOE OSTI2020

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Characterisation and testing of CHEC-M—A camera prototype for the small-sized telescopes of the Cherenkov telescope array

The Compact High Energy Camera (CHEC) is a camera design for the Small-Sized Telescopes (SSTs; 4 m diameter mirror) of the Cherenkov Telescope Array (CTA). The SSTs are focused on very-high-energy γ-ray detection via atmospheric Cherenkov light detection over a very large area. This implies many individual units and hence cost-effective implementation, as well as shower detection at large impact distance, and hence large field of view (FoV), and efficient image capture in the presence of large time gradients in the shower image detected by the camera. CHEC relies on dual-mirror optics to reduce the plate-scale and make use of 6 x 6 mm 2 pixels, leading to a low-cost (~150 k€), compact (0.5 m x 0.5 m), and light (~45 kg) camera with 2048 pixels providing a camera FoV of ~9 degrees. The CHEC electronics are based on custom TARGET (TeV array readout with GSa/s sampling and event trigger) application-specific integrated circuits (ASICs) and field programmable gate arrays (FPGAs) sampling incoming signals at a gigasample per second, with flexible camera-level triggering within a single backplane FPGA. CHEC is designed to observe in the γ-ray energy range of 1–300 TeV, and at impact distances up to 500 m. To accommodate this and provide full flexibility for later data analysis, full waveforms with 96 samples for all 2048 pixels can be read out at rates up to ~900 Hz. The first prototype, CHEC-M, based on multi-anode photomultipliers (MAPMs) as photosensors, was commissioned and characterised in the laboratory and during two measurement campaigns on a telescope structure at the Paris Observatory in Meudon. Here in this paper, the results and conclusions from the laboratory and on-site testing of CHEC-M are presented. They have provided essential input on the system design and on operational and data analysis procedures for a camera of this type. In conclusion, a second full-camera prototype based on Silicon photomultipliers (SiPMs), addressing the drawbacks of CHEC-M identified during the first prototype phase, has already been built and is currently being commissioned and tested in the laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Performance of the Muon g - 2 calorimeter and readout systems measured with test beam data

A single calorimeter station for the Muon $g-2$ experiment at Fermilab includes the following subsystems: a 54-element array of PbF$_{2}$ Cherenkov crystals read out by large-area SiPMs, bias and slow-control electronics, a suite of 800 MSPS waveform digitizers, a clock and control distribution network, a gain calibration and monitoring system, and a GPU-based frontend read out through a MIDAS data acquisition environment. The entire system performance was evaluated using 2.5 - 5 GeV electrons at the End Station Test Beam at SLAC. This paper includes a description of the individual subsystems and the results of measurements of the energy response and resolution, energy-scale stability, timing resolution, and spatial uniformity. All measured performances meet or exceed the $g-2$ experimental requirements. Based on the success of the tests, the complete production of the required 24 calorimeter stations has been made and installation into the main experiment is complete. Furthermore, the calorimeter response measurements determined here informed the design of the reconstruction algorithms that are now employed in the running $g-2$ experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2026

Designing the Protocols for Programmable Ammonia Catalysis

Programmable catalysis can provide a more energy-efficient and cost-effective route to enhancing commercial ammonia production, a key process in the advancement of renewable energy technologies and the manufacture of fertilizers and basic chemicals. This work explores the computational discovery of optimal forcing protocols to drive such dynamic catalysis models. By employing matrix-free time-stepper methods, coupled with an optimization approach, that integrates Bayesian optimization with a Bayesian continuation strategy to efficiently discover the periodic steady states of such periodically forced systems, we enable the discovery of complex optimal catalyst strain waveforms, while ensuring robust solver convergence. We demonstrate the flexibility of our approach to discover optimized forcing protocols under varying physical constraints on strain modulation or other catalyst operating parameters. We show that these can have a temporal structure more complex than simple step functions. In order to detect undesirable catalytic loops that may correlate with overall reduced performance, we perform a study using graph-theoretical analysis to investigate the dynamics of catalytic kinetic networks formed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Reverse battery model for anodic arc discharges near atmospheric pressure

In this study, the macroscopic electrical parameters of anodic arc discharges performed with solid carbon and molybdenum anodes are described in terms of an equivalent circuit. The proposed 1D model shows that the plasma volume behaves as a dissipative medium combined with self-generated voltage from the plasma sheaths. We test this model by fitting virtual circuits of resistors and batteries to the measured voltage–current (V–I) characteristics representative of atmospheric arc processes (300 and 500 Torr He), resulting in the identification of two phenomena relevant for arc technology. First, steady DC arc experiments show that a minimal arc voltage (V ≈ 20–30 V) or threshold voltage is necessary to keep the plasma arc on. Second, we prove that pulsed arc discharges excited with rectangular current waveforms (1–5 Hz) are never extinguished but remain ignited during the inactive time of the pulse as long as the aforementioned threshold voltage is maintained. In summary, an atmospheric anodic arc discharge with solid electrodes consumes electrical power as a voltage source opposed to the external voltage supply together with a resistor that accounts for the plasma impedance. This approach permits us to estimate basic plasma parameters from simple electrical diagnostics, as well as to predict how power consumption is distributed within an anodic arc and to understand dynamics of arc plasmas excited by periodic signals.

25 ENERGY STORAGE↗
DOE OSTI2026

Machine learning for seismic low-frequency extrapolation

The cycle-skipping problem that plagues full waveform inversion (FWI) can be at least partially mitigated if low frequencies (which encode the kinematics of wave propagation in seismic data) are recorded. However, seismic sources and receivers are band-limited, so seismic data does not generally include signals down to 0 Hz. To improve our ability to solve the seismic inverse problem, one can synthesize this missing low-frequency (LF) content from the recorded high-frequency (HF) data using machine learning (ML) models. Deep learning models such as convolutional neural networks (CNNs) demonstrate impressive ability to perform low frequency extrapolation. However, such models require powerful hardware (GPU machines) and careful training. We assess the extrapolation capabilities of three different ML models that do not require GPU machines, namely, random forest, Gaussian process regression and gradient boosting, on both synthetic and real data. Experimental results on two synthetic data sets (generated from a low velocity lens embedded in a homogeneous medium, and the Marmousi model) demonstrate that FWI applied to the extrapolated data consistently improves inversion accuracy relative to FWI applied to the original data sets that do not contain low frequencies. Application of low-frequency extrapolation to real data from the Northwest Shelf of Australia demonstrates that tree-based ML models such as gradient boosting can outperform CNNs in terms of both accuracy and computational cost on non-GPU architectures.

58 GEOSCIENCES↗
DOE OSTI2016

Global adjoint tomography: first-generation model

We present the first-generation global tomographic model constructed based on adjoint tomography, an iterative full-waveform inversion technique. Synthetic seismograms were calculated using GPU-accelerated spectral-element simulations of global seismic wave propagation, accommodating effects due to 3-D anelastic crust & mantle structure, topography & bathymetry, the ocean load, ellipticity, rotation, and self-gravitation. Fréchet derivatives were calculated in 3-D anelastic models based on an adjoint-state method. The simulations were performed on the Cray XK7 named ‘Titan’, a computer with 18 688 GPU accelerators housed at Oak Ridge National Laboratory. The transversely isotropic global model is the result of 15 tomographic iterations, which systematically reduced differences between observed and simulated three-component seismograms. Our starting model combined 3-D mantle model S362ANI with 3-D crustal model Crust2.0. We simultaneously inverted for structure in the crust and mantle, thereby eliminating the need for widely used ‘crustal corrections’. We used data from 253 earthquakes in the magnitude range 5.8 ≤ M w ≤ 7.0. We started inversions by combining ~30 s body-wave data with ~60 s surface-wave data. The shortest period of the surface waves was gradually decreased, and in the last three iterations we combined ~17 s body waves with ~45 s surface waves. We started using 180 min long seismograms after the 12th iteration and assimilated minor- and major-arc body and surface waves. The 15th iteration model features enhancements of well-known slabs, an enhanced image of the Samoa/Tahiti plume, as well as various other plumes and hotspots, such as Caroline, Galapagos, Yellowstone and Erebus. Furthermore, we see clear improvements in slab resolution along the Hellenic and Japan Arcs, as well as subduction along the East of Scotia Plate, which does not exist in the starting model. Point-spread function tests demonstrate that we are approaching the resolution of continental-scale studies in some areas, for example, underneath Yellowstone. Here, this is a consequence of our multiscale smoothing strategy in which we define our smoothing operator as a function of the approximate Hessian kernel, thereby smoothing gradients less wherever we have good ray coverage, such as underneath North America.

58 GEOSCIENCES↗
DOE OSTI2018

Analysis of seismic moment tensors and in situ stress during Enhanced Geothermal System development at The Geysers geothermal field, California

Seismicity in the vicinity of an Enhanced Geothermal System (EGS) demonstration project in the Northwest Geysers at the Prati 32 (P-32) injection well and the Prati State 31 (PS-31) production well is investigated to determine earthquake focal mechanisms, moment magnitudes and in situ stress of events that occur before and during reservoir stimulation starting 2011 October 6. We present a catalogue of 167 waveform-based seismic moment tensors ranging in moment magnitude from 0.6 to 3.9. The moment tensor catalogue is subsequently used to invert for the stress tensor and to investigate possible temporal stress changes resulting from fluid injection. We find an approximate 15-degree counter-clockwise rotation of the least compressive stress σ3 during injection. More remarkable is a change in orientation of the maximum compressive stress σ1 from subhorizontal to vertical when injection operations temporarily cease. The orientation of σ1 returns to subhorizontal when injection operations resume. It is found that there was a systematic reduction in the stress shape factor, R, as the injected water volume increases, indicating an evolution towards a more transtensional stress state.

58 GEOSCIENCES↗
DOE OSTI2017

Decentralized Interleaving of Paralleled Dc-Dc Buck Converters

We present a decentralized control strategy that yields switch interleaving among parallel-connected dc-dc buck converters. The proposed method is based on the digital implementation of the dynamics of a nonlinear oscillator circuit as the controller. Each controller is fully decentralized, i.e., it only requires the locally measured output current to synthesize the pulse width modulation (PWM) carrier waveform and no communication between different controllers is needed. By virtue of the intrinsic electrical coupling between converters, the nonlinear oscillator-based controllers converge to an interleaved state with uniform phase-spacing across PWM carriers. To the knowledge of the authors, this work presents the first fully decentralized strategy for switch interleaving in paralleled dc-dc buck converters.

14 SOLAR ENERGY↗
DOE OSTI2022

Making Invisible Visible: Data-Driven Seismic Inversion With Spatio-Temporally Constrained Data Augmentation

Deep learning and data-driven approaches have shown great potential in scientific domains. The promise of data-driven techniques relies on the availability of a large volume of high-quality training datasets. Due to the high cost of obtaining data through expensive physical experiments, instruments, and simulations, data augmentation techniques for scientific applications have emerged as a new direction for obtaining scientific data recently. However, existing data augmentation techniques originating from computer vision yield physically unacceptable data samples that are not helpful for the domain problems that we are interested in. In this article, we develop new data augmentation techniques based on convolutional neural networks. Specifically, our generative models leverage different physics knowledge (such as governing equations, observable perception, and physics phenomena) to improve the quality of the synthetic data. To validate the effectiveness of our data augmentation techniques, we apply them to solve a subsurface seismic full-waveform inversion using simulated CO 2 leakage data. Our interest is to invert for subsurface velocity models associated with very small CO 2 leakage. We validate the performance of our methods using comprehensive numerical tests. Here via comparison and analysis, we show that data-driven seismic imaging can be significantly enhanced by using our data augmentation techniques. Particularly, the imaging quality has been improved by 15% in test scenarios of general-sized leakage and 17% in small-sized leakage when using an augmented training set obtained with our techniques.

58 GEOSCIENCES↗
DOE OSTI2020

Time Resolution Measurements of EJ-232Q With Single- and Dual-Sided Readouts

In view of its importance to applications in high energy physics, nuclear security, and time-of-flight PET, the time resolution of 0.5% benzophenone-quenched EJ-232Q plastic scintillators of various lengths was studied using a fast waveform digitizer and fast PMTs. The scintillator sizes were 5 mm × 5 mm × 10, 25, 50, 75, or 100 mm. First, the timing measurements were performed by using a single PMT reading out one end. The degradation in time resolution with increasing scintillator length was studied. The irradiation geometry was also observed to have a non-negligible impact on the time resolution due to the variation of depth-of-interaction positions, especially for the scintillators with a high aspect ratio. Second, two PMTs were employed for a dual-sided readout. This readout together with the mean time method was demonstrated to be an effective approach to improve the timing performance since the depth-of-interaction-related effect was considerably mitigated. A time resolution of 79 ± 1 ps (full-width at half-maximum) has been achieved for the 5 × 5 × 10 mm 3 scintillator. The 5 × 5 × 100 mm 3 scintillator can still provide a time resolution as good as 172 ± 1 ps. Here, the energy range used for each was between 340 keVee and 450 keVee. In addition, the time information from individual PMTs in a dual-sided readout was able to show a time resolution almost identical to that determined with the mean-time-based technique after correcting for the depth-of-interaction variation for the 100 mm long scintillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Multistep inversion workflow for 3D long-offset damped elastic waves in the Fourier domain

We present a new workflow for imaging damped three-dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short-offset and low-frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine-scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite-difference frequency-domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full-waveform inversion workflow is applied to the 3D Marmousi-2 and SEG/EAGE Salt models with long-offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one-third of the maximum offset length at a resolution consistent with the inverted frequencies.

3D full waveform inversion↗
DOE OSTI2022

Flicker: A review of temporal light modulation stimulus, responses, and measures

Flicker has been an important lighting system consideration for over a century. More precise terms are temporal light modulation (TLM) as the stimulus, and responses to TLM as the unwanted visual, cognitive, or physiological consequences. As lighting technology evolved, different forms of TLM emerged, and so did responses to them. Today, some LED systems – encompassing the LED, driver, and control – can result in TLM causing severe unwanted effects, while other LED systems produce no unwanted effects at all. LED systems can deliver a much wider range of luminous waveforms than conventional lighting systems, some exhibiting very high modulation depths. More than any light source before, they can elicit perceptions of the phantom array. Direct flicker effects at modulation frequencies less than about 80 Hz and the stroboscopic effect at frequencies greater than 80 Hz are fairly well understood, but the phantom array effect needs more exploration and characterisation. This review focuses on the technology and research history that led to current metrics for quantifying TLM and human responses to TLM. Notably, visually impaired individuals may exhibit alterations in their response to TLM, but such a discussion is beyond the intent of this review. Thus, the focus is on individuals with normal visual function.

60 APPLIED LIFE SCIENCES↗