Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wavelet functions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Scalability of a Low-Cost Multi-Teraflop Linux Cluster for High-End Classical Atomistic and Quantum Mechanical Simulations

Scalability of a low-cost, Intel Xeon-based, multi-Teraflop Linux cluster is tested for two high-end scientific applications: Classical atomistic simulation based on the molecular dynamics method and quantum mechanical calculation based on the density functional theory. These scalable parallel applications use space-time multiresolution algorithms and feature computational-space decomposition, wavelet-based adaptive load balancing, and spacefilling-curve-based data compression for scalable I/O. Comparative performance tests are performed on a 1,024-processor Linux cluster and a conventional higher-end parallel supercomputer, 1,184-processor IBM SP4. The results show that the performance of the Linux cluster is comparable to that of the SP4. We also study various effects, such as the sharing of memory and L2 cache among processors, on the performance.

Kikuchi, Hideaki↗
DOE OSTI2020

Software-Defined Ultrasonic Communication System Based on Time-reversal Signal Processing

Ultrasonic communication through solid channels suffers from channel dispersion, attenuation, and multipath effect. For example, in elbow shape pipes or complex structures, the elastic waves are reflected and reverberated in multipath resulting in a major dispersion of the signals. In solid channels such as flat bar, multipath propagation is the main reason for bit-wavelet spreading and intersymbol interference (ISI). The conventional modulation method like OOK (on-off-keying) or PSK (phase-shift-keying) fails to provide a useful Signal-to-Noise Ratio (SNR). Time-reversal (TR) signal processing method is a promising technique that can compress the spreading of the bit-wavelet and strengthen the peak signal energy on the receiver. In this study, Finite Element Analysis (FEA) using ABAQUS software tools is used to characterize the bit-wavelet spreading effect caused by dispersive communication channels. Also, this study presents the impact of TR ultrasonic communication as a function of excitation pulse width, channel length, and carrier frequency through a flat-bar channel. A Software-Defined System-on-Chip (SD-SoC) platform is examined as a viable solution for the low-complexity and energy-efficient ultrasonic communication system. The proposed TR based SD-SoC ultrasonic communication system can make full use of signal dispersion for peak performance by minimizing the ISI within a dispersive and reverberant metal channel.

FEM simulation↗
NASA NTRS1996

Perceptually Lossless Wavelet Compression

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp -1), where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We propose a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a 'perceptually lossless' quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗
NASA NTRS1995

Visibility of Wavelet Quantization Noise

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp)-L , where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We describe a mathematical model to predict DWT noise detection thresholds as a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗
NASA NTRS2004

Complexity Induced Anisotropic Bimodal Intermittent Turbulence in Space Plasmas

The "physics of complexity" in space plasmas is the central theme of this exposition. It is demonstrated that the sporadic and localized interactions of magnetic coherent structures arising from the plasma resonances can be the source for the coexistence of nonpropagating spatiotemporal fluctuations and propagating modes. Non-Gaussian probability distribution functions of the intermittent fluctuations from direct numerical simulations are obtained and discussed. Power spectra and local intermittency measures using the wavelet analyses are presented to display the spottiness of the small-scale turbulent fluctuations and the non-uniformity of coarse-grained dissipation that can lead to magnetic topological reconfigurations. The technique of the dynamic renormalization group is applied to the study of the scaling properties of such type of multiscale fluctuations. Charged particle interactions with both the propagating and nonpropagating portions of the intermittent turbulence are also described.

Chang, Tom↗
DOE OSTI2024

Denoising Seismic Waveforms Using a Wavelet-Transform-Based Machine-Learning Method

Seismic waveform data recorded at stations can be thought of as a superposition of the signal from a source of interest and noise from other sources. Frequency‐based filtering methods for waveform denoising do not result in desired outcomes when the targeted signal and noise occupy similar frequency bands. Recently, denoising techniques based on deep‐learning convolutional neural networks (CNNs), in which a recorded waveform is decomposed into signal and noise components, have led to improved results. These CNN methods, which use short‐time Fourier transform representations of the time series, provide signal and noise masks for the input waveform. These masks are used to create denoised signal and designaled noise waveforms, respectively. However, advancements in the field of image denoising have shown the benefits of incorporating discrete wavelet transforms (DWTs) into CNN architectures to create multilevel wavelet CNN (MWCNN) models. The MWCNN model preserves the details of the input due to the good time–frequency localization of the DWT. In this report we use a data set of over 382,000 constructed seismograms recorded by the University of Utah Seismograph Stations network to compare the performance of CNN and MWCNN‐based denoising models. Evaluation of both models on constructed test data shows that the MWCNN model outperforms the CNN model in the ability to recover the ground‐truth signal component in terms of both waveform similarity and preservation of amplitude information. Model evaluation of real‐world data shows that both the CNN and MWCNN models outperform standard band‐pass filtering (BPF; average improvement in signal‐to‐noise ratio of 9.6 and 19.7 dB, respectively, with respect to BPF). Evaluation of continuous data suggests the MWCNN denoiser can improve both signal detection capabilities and phase arrival time estimates.

58 GEOSCIENCES↗
DOE OSTI2025

The hidden influence of terrestrial groundwater on salt marsh function and resilience

Salt marshes are hotspots of nutrient processing and carbon sequestration. So far, studies addressing spatiotemporal variability in and drivers of salt marsh biogeochemical function, carbon storage and resilience have focused on ocean-driven surface hydrologic influences, neglecting effects of terrestrial hydrology through subsurface connections. Here we evaluate drivers of salt marsh redox potential, a proxy for biogeochemical state, through wavelet analyses and information theory using data from seven marshes. The results point to terrestrial groundwater level as a dominant control on redox variability across all sites. Because redox is a key driver of biogeochemical processes, and specifically oxidation of organic matter that sequesters carbon and maintains marsh elevation, these terrestrial influences are critical to understanding marsh function and evolution. Further, the newly identified links between onshore groundwater levels and marsh redox conditions shift the traditional paradigm and suggest that terrestrial hydrology is a primary control on salt marsh carbon sequestration potential and resilience.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Constraining the Multiscale Structure of Geophysical Fields in Machine Learning: The Case of Precipitation

The use of deep-learning algorithms for estimating the value of geophysical variables from remotely sensed information is rapidly expanding. The typical objective function minimized in such algorithms is the mean square error (MSE), which is known to lead to smooth estimates with compressed dynamical range as compared to the true distribution of the variable of interest. Here, we introduce and evaluate alternative objective functions, focusing on the retrieval of precipitation rates from satellite passive microwave radiometric measurements using a deep convolutional neural network. For this testbed application, the results show that explicitly imposing the preservation of the statistical distribution and spatial wavelet power spectrum of the target variable allows to accurately reproduce extreme values and sharp gradients across multiple scales in this study.

58 GEOSCIENCES↗
NASA NTRS2010

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI↗
NASA NTRS1997

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↗
NASA NTRS2008

Image Compression Algorithm Altered to Improve Stereo Ranging

A report discusses a modification of the ICER image-data-compression algorithm to increase the accuracy of ranging computations performed on compressed stereoscopic image pairs captured by cameras aboard the Mars Exploration Rovers. (ICER and variants thereof were discussed in several prior NASA Tech Briefs articles.) Like many image compressors, ICER was designed to minimize a mean-square-error measure of distortion in reconstructed images as a function of the compressed data volume. The present modification of ICER was preceded by formulation of an alternative error measure, an image-quality metric that focuses on stereoscopic-ranging quality and takes account of image-processing steps in the stereoscopic-ranging process. This metric was used in empirical evaluation of bit planes of wavelet-transform subbands that are generated in ICER. The present modification, which is a change in a bit-plane prioritization rule in ICER, was adopted on the basis of this evaluation. This modification changes the order in which image data are encoded, such that when ICER is used for lossy compression, better stereoscopic-ranging results are obtained as a function of the compressed data volume.

Kiely, Aaron↗
NASA NTRS1985

Specification of Surface Roughness Over Oceans in the GLAS Fourth Order GCM

The surface roughness height z sub 0 is a parameter which measures the effect of surface irregularities on the mean profiles which occur in the atmospheric surface layer. Over land surfaces, z sub 0 is a function of soil composition, plant canopy type and distribution and other factors such as buildings and structures. Over open water, z sub 0 depends only on the distribution of waves, wavelets and other surface irregularities. These, in turn, are thought to be determined by the mean wind within the surface layer. Large and Pond have derived an empirical relationship between the mean wind speed U sub 10 at the ten meter level over open water and the neutral drag coefficient C sub DN. Their relationship covers the range U sub 10 or = M/S. Earlier, Kondo had computed a relationship between U sub 10 and C sub DN on the basis of laboratory observations which extended down to the range of infinitely small wind speeds. An interpretation between the relationships of Large and Pond for large wind speed and of Kondo for small wind speed is presented.

Helfand, H. M.↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗
NASA NTRS2018

The Impact of Stellar Variability on the Detection of Transiting Earth-like Planets

Stellar variability represents a key limitation on the detectability of weak transit signatures of small, Earth-size planets with space-based missions such as NASA's Kepler Mission or ESA's PLATO Mission. The expectations for the performance of the Kepler Mission in the face of solar-like variability were predicated on observations of the Sun with the Variability of Solar Irradiance and Gravity Oscillations (VIRGO) instrument aboard NASA and ESA's joint Solar and Heliospheric Observatory (SOHO) mission, which was launched in December 1995. Analyses of SOHO data indicated that solar-like variability would contribute approx.10 ppm to Kepler's noise budget at 6.5 ours, a typical grazing transit timescale. This proved to be optimistic as the typical 12th mag G2V star in the Kepler field of view exhibited 20 ppm of stellar variability-related noise at 6.5 hours, making the detection of true Earth-Sun analogs more difficult. While the consequent increase in the total noise at 6.5 hours from 20 ppm to 30 ppm could have been overcome by extending the mission to a total duration of 7 to 8 years. Unfortunately, a second reaction wheel failed in May 2013, ending the Kepler primary mission. Nevertheless, Kepler proved to be a veritable cornucopia of science results, both for exoplanets and for astrophysics. The phenomenal photometric precision and continuous observations required in order to identify small, rocky transiting planets enables the study of a large range of phenomena contributing to stellar variability for many thousands of solar-like stars in Kepler's field of view in exquisite detail. These effects range from less than 1 ppm acoustic oscillations on timescales from a few minutes and longward, to flares on timescales of hours, to spot-induced modulation on timescales of days to weeks to activity cycles on timescales of months to years. Kepler discovered over 2600 validated and confirmed exoplanets, and measured the pressure-mode oscillations of over 15000 stars over the course of its mission and the reconstituted K2 mission. We present the adaptive, wavelet-based matched filter used for both the Kepler and the TESS missions, and discuss the connections between this detection algorithm and fundamental detection theory, which allows the detector to function as a noise characterization engine. This provides a dynamical measurement of the photometric noise at transit timescales and thus, a key performance metric for missions like PLATO and Kepler. We give an overview of the stellar variability we see across the full range of spectral types observed by Kepler, from the cool, small red M stars to the hot, large late A stars, both in terms of amplitude as well as timescale.

transit surveys↗
NASA NTRS2018

The Impact of Stellar Variability on the Detection of Transiting Earth-Like Planets

Stellar variability represents a key limitation on the detectability of weak transit signatures of small, Earth-size planets with space-based missions such as NASA's Kepler Mission or ESA's PLATO Mission. The expectations for the performance of the Kepler Mission in the face of solar-like variability were predicated on observations of the Sun with the Variability of Solar Irradiance and Gravity Oscillations (VIRGO) instrument aboard NASA and ESA's joint Solar and Heliospheric Observatory (SOHO) mission, which was launched in December 1995. Analyses of SOHO data indicated that solar-like variability would contribute ~10 ppm to Kepler's noise budget at 6.5 ours, a typical grazing transit timescale. This proved to be optimistic as the typical 12th mag G2V star in the Kepler field of view exhibited 20 ppm of stellar variability-related noise at 6.5 hours, making the detection of true Earth-Sun analogs more difficult. While the consequent increase in the total noise at 6.5 hours from 20 ppm to 30 ppm could have been overcome by extending the mission to a total duration of 7 to 8 years. Unfortunately, a second reaction wheel failed in May 2013, ending the Kepler primary mission. Nevertheless, Kepler proved to be a veritable cornucopia of science results, both for exoplanets and for astrophysics. The phenomenal photometric precision and continuous observations required in order to identify small, rocky transiting planets enables the study of a large range of phenomena contributing to stellar variability for many thousands of solar-like stars in Kepler's field of view in exquisite detail. These effects range from <1 ppm acoustic oscillations on timescales from a few minutes and longward, to flares on timescales of hours, to spot-induced modulation on timescales of days to weeks to activity cycles on timescales of months to years. Kepler discovered over 2600 validated and confirmed exoplanets, and measured the pressure-mode oscillations of over 15000 stars over the course of its mission and the reconstituted K2 mission. We present the adaptive, wavelet-based matched filter used for both the Kepler and the TESS missions, and discuss the connections between this detection algorithm and fundamental detection theory, which allows the detector to function as a noise characterization engine. This provides a dynamical measurement of the photometric noise at transit timescales and thus, a key performance metric for missions like PLATO and Kepler. We give an overview of the stellar variability we see across the full range of spectral types observed by Kepler, from the cool, small red M stars to the hot, large late A stars, both in terms of amplitude as well as timescale.

exoplanets↗
NASA NTRS2001

Normal and abnormal tissue identification system and method for medical images such as digital mammograms

A system and method for analyzing a medical image to determine whether an abnormality is present, for example, in digital mammograms, includes the application of a wavelet expansion to a raw image to obtain subspace images of varying resolution. At least one subspace image is selected that has a resolution commensurate with a desired predetermined detection resolution range. A functional form of a probability distribution function is determined for each selected subspace image, and an optimal statistical normal image region test is determined for each selected subspace image. A threshold level for the probability distribution function is established from the optimal statistical normal image region test for each selected subspace image. A region size comprising at least one sector is defined, and an output image is created that includes a combination of all regions for each selected subspace image. Each region has a first value when the region intensity level is above the threshold and a second value when the region intensity level is below the threshold. This permits the localization of a potential abnormality within the image.

Heine, John J.↗
NASA NTRS2016

Unsteady Heat-Flux Measurements of Second-Mode Instability Waves in a Hypersonic Boundary Layer

In this paper we report on the application of the atomic layer thermopile (ALTP) heat- flux sensor to the measurement of laminar-to-turbulent transition in a hypersonic flat plate boundary layer. The centerline of the flat-plate model was instrumented with a streamwise array of ALTP sensors and the flat-plate model was exposed to a Mach 6 freestream over a range of unit Reynolds numbers. Here, we observed an unstable band of frequencies that are associated with second-mode instability waves in the laminar boundary layer that forms on the flat-plate surface. The measured frequencies, group velocities, phase speeds, and wavelengths of these instability waves are in agreement with data previously reported in the literature. Heat flux time series, and the Morlet-wavelet transforms of them, revealed the wave-packet nature of the second-mode instability waves. In addition, a laser-based radiative heating system was developed to measure the frequency response functions (FRF) of the ALTP sensors used in the wind tunnel test. These measurements were used to assess the stability of the sensor FRFs over time and to correct spectral estimates for any attenuation caused by the finite sensor bandwidth.

Kergerise, Michael A.↗