Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wavelet Transform”

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 217 records · Page 12

Stratospheric Horizontal Wavenumber Spectra of Winds, Potential Temperature, and Atmospheric Tracers Observed by High-Altitude Aircraft

Horizontal wavenumber power spectra of vertical and horizontal wind velocities, potential temperatures, and ozone and N(2)O mixing ratios, as measured in the mid-stratosphere during 73 ER-2 flights (altitude approx. 20km) are presented. The velocity and potential temperature spectra in the 100 to 1-km wavelength range deviate significantly from the uniform -5/3 power law expected for the inverse energy-cascade regime of two-dimensional turbulence and also for inertial-range, three-dimensional turbulence. Instead, steeper spectra approximately consistent with a -3 power law are observed at horizontal scales smaller than 3 km for all velocity components as well as potential temperature. Shallower spectra are observed at scales longer than 6 km. For horizontal velocity and potential temperature the spectral indices at longer scales are between -1.5 and -2.0. For vertical velocity the spectrum at longer scales become flat. It is argued that the observed velocity and potential temperature spectra are consistent with gravity waves. At smaller scales, the shapes are also superficially consistent with a Lumley-Shur-Weinstock buoyant subrange of turbulence and/or nonlinear gravity waves. Contemporaneous spectra of ozone and N(sub 2)O mixing ratio in the 100 to 1-km wavelength range do conform to an approximately uniform -5/3 power law. It is argued that this may reflect interactions between gravity wave air-parcel displacements and laminar or filamentary structures in the trace gas mixing ratio field produced by enstropy-cascading two-dimensional turbulence.

TRACE GAS MIXING RATIO FIELDS↗

Evaluating Scalograms for Seismic Event Denoising

Denoising contaminated seismic signals for later processing is a fundamental problem in seismic signals analysis. The most straightforward denoising approach, using spectral filtering, is not effective when noise and seismic signal occupy the same frequency range. Neural network approaches have shown success denoising local signal when trained on short-time Fourier transform spectrograms (Zhu et al 2018; Tibi et al 2021). Scalograms, a wavelet-based transform, achieved ~15% better reconstruction as measured by dynamic time warping on a seismic waveform test set than spectrograms, suggesting their use as an alternative for denoising. We train a deep neural network on a scalogram dataset derived from waveforms recorded by the University of Utah Seismograph Stations network. We find that initial results are no better than a spectrogram approach, with additional overhead imposed by the significantly larger size of scalograms. A robust exploration of neural network hyperparameters and network architecture was not performed, which could be done in follow on work.

58 GEOSCIENCES↗

A Step-Down Test Procedure for Wavelet Shrinkage Using Bootstrapping

Wavelet thresholding (or shrinkage) attempts to remove the noises existing in the signals while preserving inherent pattern characteristics in the reconstruction of true signals. For data-denoising purpose, we present a new wavelet thresholding procedure which employs the step-down testing idea of identifying active contrasts in unreplicated fractional factorial experiments. The proposed method employs bootstrapping methods to a step-down test for thresholding wavelet coefficients. By introducing the concept of a false discovery error rate in testing wavelet coefficients, we shrink the wavelet coefficients with p -values higher than the error rate. The error rate controls the expected proportion of wrongly accepted coefficients among chosen wavelet coefficients. Bootstrap samples are used to approximate the p -value for computational efficiency. We also present some guidelines for selecting the values of hyper-parameters which affect the performance in the step-down thresholding procedure. Based on some common testing signals and an air-conditioner sounds example, the comparison of our proposed procedure with other thresholding methods in the literature is performed. The analytical results show that the proposed procedure has a potential in data-denoising and data-reduction in a variety of signal reconstruction applications.

42 ENGINEERING↗

Wavelet Analyses of F/A-18 Aeroelastic and Aeroservoelastic Flight Test Data

Time-frequency signal representations combined with subspace identification methods were used to analyze aeroelastic flight data from the F/A-18 Systems Research Aircraft (SRA) and aeroservoelastic data from the F/A-18 High Alpha Research Vehicle (HARV). The F/A-18 SRA data were produced from a wingtip excitation system that generated linear frequency chirps and logarithmic sweeps. HARV data were acquired from digital Schroeder-phased and sinc pulse excitation signals to actuator commands. Nondilated continuous Morlet wavelets implemented as a filter bank were chosen for the time-frequency analysis to eliminate phase distortion as it occurs with sliding window discrete Fourier transform techniques. Wavelet coefficients were filtered to reduce effects of noise and nonlinear distortions identically in all inputs and outputs. Cleaned reconstructed time domain signals were used to compute improved transfer functions. Time and frequency domain subspace identification methods were applied to enhanced reconstructed time domain data and improved transfer functions, respectively. Time domain subspace performed poorly, even with the enhanced data, compared with frequency domain techniques. A frequency domain subspace method is shown to produce better results with the data processed using the Morlet time-frequency technique.

Martin J Brenner↗

WigglyRivers: A tool to characterize the multiscale nature of meandering channels

Channel sinuosity is ubiquitous along river networks, producing complex patterns that encapsulate and influence morphodynamic processes and ecosystem services. Accurately characterizing these patterns is challenging with traditional curvature-based algorithms. Here, in this study, we present WigglyRivers, a Python package that builds on existing wavelet-based methods to create an unsupervised meander identification and characterization tool. The package uses planimetric information the user provides or from the USGS’s High-Resolution National Hydrography Dataset to characterize individual reaches or entire river networks. WigglyRivers also includes a supervised river identification tool for manually selecting individual meandering features. Here, we provide examples of idealized river transects and show the capabilities of WigglyRivers. We also use the supervised identification tool to validate the unsupervised identification on river transects across the continental US. WigglyRivers is a tool to understand better the multiscale characteristics of river networks and the link between river geomorphology and river corridor connectivity.

54 ENVIRONMENTAL SCIENCES↗

AI-enabled traveling wave protection for microgrids

Grid forming inverters provide voltage and frequency regulations for microgrids; in the meantime, new challenges are introduced for microgrid protections. For instance, inverters’ control strategies can affect protection behaviors, and low short-circuit ratios and bi-directional power flows also make protection operations complex. Protection schemes based on conventional principles such as overcurrent and distance relays do not always provide reliable, sensitive, or selective operations. Here, we devise a traveling wave protection approach for microgrids using a wavelet-driven deep neural network named WaveletKernelNet (WKN). Compared with conventional methods, the presented approach provides enhanced sensitivity, higher selectivity, and better identification of various faults in microgrids. Extensive case studies validate the efficacy and excellent performance of the devised approach.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Preliminary study of radio frequency waves in hypervelocity impact plasma

Orbital debris and meteoroids can cause mechanical and electrical anomalies by impacting spacecraft. Although mechanical damage tends to occur only from particles with masses >10 –3 g, electrical damage can be a result of hypervelocity (> 10 km/s) impacts (HVI) even with smaller particles. In HVI, the projectile speed exceeds the speed of sound in the material and has enough energy to create plasma by ionizing the material near the surface, which can then emit radio frequency (RF) waves. In this paper, we conduct a preliminary study on the characteristics of impact-generated RF emission with the ultimate goal of building an RF predictive model for risk assessment of HVI electrical effects that can damage spacecraft. First, we introduce wavelet analysis as a method for denoising and quantifying the total RF signal energy. An appropriate mother wavelet and the number of decomposition levels are selected using a synthetic signal. Second, the total energy of the RF emission is calculated using Parseval’s theorem. Results suggest that the total RF signal energy is proportional to the energy flux of the impactor. Lastly, spectrograms of the RF emission are analyzed to characterize features that can be used in the future to predict impacts that can cause electrical damage to satellites.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗

Stabilizing dynamic subsea power cables using Bi-stable nonlinear energy sinks

This study investigates vibration mitigation of dynamic subsea cables through passive bi-stable nonlinear energy sinks (B-NESs). These devices suppress vibration energy in a broadband way, and can be regarded as extensions of classical linear tuned mass dampers (TMDs) which are narrowband devices. Through the open-source MoorDyn library, we simulated the vibrations of a vertical subsea cable equipped with a set of B-NESs. Multi-objective optimization was performed to detect the B-NES parameters for optimal mitigation of the cable vibrations. Advanced signal processing verified the efficacy of the optimized B-NESs not only to rapidly absorb and locally dissipate vibration energy, but also to nonlinearly scatter vibration energy from low to high frequencies within the cable itself. This last feature is especially beneficial for vibration mitigation of the undersea cable, as at higher frequencies the cable vibrations exhibit drastically reduced amplitudes and are more effectively dissipated by inherent structural damping and hydrodynamic radiation damping. This contrasts with traditional TMDs which can mitigate vibration only at a single frequency. Furthermore, our robustness study confirms the B-NES's effectiveness under even varying environmental conditions. Overall, the B-NES's capacity for broadband vibration mitigation renders it a promising retrofit solution for improving the performance and operational safety of dynamic power cables in offshore wind farms and other marine applications.

17 WIND ENERGY↗

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient and Flexible Hierarchical Data Layouts for a Unified Encoding of Scalar Field Precision and Resolution

To address the problem of ever-growing scientific data sizes making data movement a major hindrance to analysis, we introduce a novel encoding for scalar fields: a unified tree of resolution and precision, specifically constructed so that valid cuts correspond to sensible approximations of the original field in the precision-resolution space. Furthermore, we introduce a highly flexible encoding of such trees that forms a parameterized family of data hierarchies. We discuss how different parameter choices lead to different trade-offs in practice, and show how specific choices result in known data representation schemes such as zfp[52], idx[58], and jpeg2000 [76]. Lastly, we provide system-level details and empirical evidence on how such hierarchies facilitate common approximate queries with minimal data movement and time, using real-world data sets ranging from a few gigabytes to nearly a terabyte in size. Experiments suggest that our new strategy of combining reductions in resolution and precision is competitive with state-of-the-art compression techniques with respect to data quality, while being significantly more flexible and orders of magnitude faster, and requiring significantly reduced resources.

97 MATHEMATICS AND COMPUTING↗

Quantized Information in Spectral Cyberspace

The constant-Q Gabor atom is developed for spectral power, information, and uncertainty quantification from time–frequency representations. Stable multiresolution spectral entropy algorithms are constructed with continuous wavelet and Stockwell transforms. The recommended processing and scaling method will depend on the signature of interest, the desired information, and the acceptable levels of uncertainty of signal and noise features. Selected Lamb wave signatures and information spectra from the 2022 Tonga eruption are presented as representative case studies. Resilient transformations from physical to information metrics are provided for sensor-agnostic signal processing, pattern recognition, and machine learning applications.

74 ATOMIC AND MOLECULAR PHYSICS↗

Enhancing and Archiving the APS Catalog of the POSS I

We have worked on two different projects: 1) Archiving the APS Catalog of the POSS I for distribution to NASA's NED at IPAC, SIMBAD in France, and individual astronomers and 2) The automated morphological classification of galaxies. We have completed archiving the Catalog into easily readable binary files. The database together with the software to read it has been distributed on DVD's to the national and international data centers and to individual astronomers. The archived Catalog contains more than 89 million objects in 632 fields in the first epoch Palomar Observatory Sky Survey. Additional image parameters not available in the original on-line version are also included in the archived version. The archived Catalog is also available and can be queried at the APS web site (URL: http://aps.umn.edu) which has been improved with a much faster and more efficient querying system. The Catalog can be downloaded as binary datafiles with the source code for reading it. It is also being integrated into the SkyQuery system which includes the Sloan Digital Sky Survey, 2MASS, and the FIRST radio sky survey. We experimented with different classification algorithms to automate the morphological classification of galaxies. This is an especially difficult problem because there are not only a large number of attributes or parameters and measurement uncertainties, but also the added complication of human disagreement about the adopted types. To solve this problem we used 837 galaxy images from nine POSS I fields at the North Galactic Pole classified by two independent astronomers for which they agree on the morphological types. The initial goal was to separate the galaxies into the three broad classes relevant to issues of large scale structure and galaxy formation and evolution: early (ellipticals and lenticulars), spirals, and late (irregulars) with an accuracy or success rate that rivals the best astronomer classifiers. We also needed to identify a set of parameters derived from the digitized images that separate the galaxies by type. The human eye can easily recognize complicated patterns in images such as spiral arms which can be spotty, blotchy affairs that are difficult for automated techniques. A galaxy image can potentially be described by hundreds of parameters, all of which may have some relation to the morphological type. In the set of initial experiments we used 624 such parameters, in two colors, blue and red. These parameters include the surface brightness and color measured at different radii, ratios of these parameters at different radii, concentration indices, Fourier transforms and wavelet decomposition coefficients. We experimented with three different classes of classification algorithms; decision trees, k-nearest neighbors, and support vector machines (SVM). A range of experiments were conducted and we eventually narrowed the parameters to 23 selected parameters. SVM consistently outperformed the other algorithms with both sets of features. By combining the results from the different algorithms in a weighted scheme we achieved an overall classification success of 86%.

Humphreys, Roberta M.↗

Optical Implementation of Matching Pursuit for Image Representation

We have developed a technique for image analysis, representation, and decomposition. This technique was motivated by Stephane Mallat's matching-pursuit algorithm. We've altered and simplified the mechanics of his algorithm to enable an extremely fast implementation via optical processing. Initial computer simulations show that our algorithm is capable of decomposing and representing a 2-D image as a linear combination of basis images with both high speed and high fidelity.

matching pursuit adaptive wavelet transform image ↗

Wavelet compression techniques for hyperspectral data

Hyperspectral sensors are electro-optic sensors which typically operate in visible and near infrared bands. Their characteristic property is the ability to resolve a relatively large number (i.e., tens to hundreds) of contiguous spectral bands to produce a detailed profile of the electromagnetic spectrum. In contrast, multispectral sensors measure relatively few non-contiguous spectral bands. Like multispectral sensors, hyperspectral sensors are often also imaging sensors, measuring spectra over an array of spatial resolution cells. The data produced may thus be viewed as a three dimensional array of samples in which two dimensions correspond to spatial position and the third to wavelength. Because they multiply the already large storage/transmission bandwidth requirements of conventional digital images, hyperspectral sensors generate formidable torrents of data. Their fine spectral resolution typically results in high redundancy in the spectral dimension, so that hyperspectral data sets are excellent candidates for compression. Although there have been a number of studies of compression algorithms for multispectral data, we are not aware of any published results for hyperspectral data. Three algorithms for hyperspectral data compression are compared. They were selected as representatives of three major approaches for extending conventional lossy image compression techniques to hyperspectral data. The simplest approach treats the data as an ensemble of images and compresses each image independently, ignoring the correlation between spectral bands. The second approach transforms the data to decorrelate the spectral bands, and then compresses the transformed data as a set of independent images. The third approach directly generalizes two-dimensional transform coding by applying a three-dimensional transform as part of the usual transform-quantize-entropy code procedure. The algorithms studied all use the discrete wavelet transform. In the first two cases, a wavelet transform coder was used for the two-dimensional compression. The third case used a three dimensional extension of this same algorithm.

Evans, Bruce↗

An Investigation into the Potential Application of Wavelets to Modal Testing and Analysis

The analysis of transient data of the type found in vibrating mechanical systems has been greatly improved through the use of modern techniques such as Fourier analysis. This is especially true when considered in conjunction with the development of the so-called Fast Fourier Transform algorithm by Cooley and the tremendous strides in computational power of the last several decades. The usefulness of the discrete Fourier Transform is its ability to transform sampled data from the "time-domain" to the "frequency domain," thereby allowing the analyst to decompose a signal into its frequency content. More recent developments have led to the wavelet transform. The strength of wavelet analysis is its ability to maintain both time and frequency information, thus making it an attractive candidate for the analysis of non-stationary signals. This report is an overview of wavelet theory and the potential use of the wavelet transform as an alternative to Fourier analysis in modal identification.

A. Fort Gwinn, Jr.↗