Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time-frequency analysis”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Development of a Model Based Technique for Gear Diagnostics using the Wigner-Ville method

Imperfections in gear tooth geometry often result from errors in the manufacturing process or excessive material wear during operation. Such faults in the gear tooth geometry can result in large vibrations in the transmission system, and, in some cases, may lead to early failure of the gear transmission system. This report presents the study of the effects of imperfection in gear tooth geometry on the dynamic characteristics of a gear transmission system. The faults in the gear tooth geometry are modeled numerically as the deviation of the tooth profile from its original involute geometry. The changes in gear mesh stiffness due to various profile and pattern variations are evaluated numerically. The resulting changes in the mesh stiffness are incorporated into a computer code to simulate the dynamics of the gear transmission system. A parametric study is performed to examine the sensitivity of gear tooth geometry imperfections on the vibration of a gear transmission system. The parameters variations in this study consist of the magnitude of the imperfection, the pattern of the profile variation, and the total number of teeth affected. Numerical results from the dynamic simulations are examined in both the time and the frequency domains. A joint time-frequency analysis procedure using the Wigner-Ville Distribution is also introduced to identify the location of the damaged tooth from the vibration signature. Numerical simulations of the system dynamics with gear faults were compared to experimental results. An optimal tracker was introduced to quantify the level of damage in the gear mesh system. Conclusions are drawn from the results of this numerical study.

Choy, F.↗
NASA NTRS2001

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.↗
NASA NTRS2005

Application of the Hilbert-Huang Transform to Financial Data

A paper discusses the application of the Hilbert-Huang transform (HHT) method to time-series financial-market data. The method was described, variously without and with the HHT name, in several prior NASA Tech Briefs articles and supporting documents. To recapitulate: The method is especially suitable for analyzing time-series data that represent nonstationary and nonlinear phenomena including physical phenomena and, in the present case, financial-market processes. The method involves the empirical mode decomposition (EMD), in which a complicated signal is decomposed into a finite number of functions, called "intrinsic mode functions" (IMFs), that admit well-behaved Hilbert transforms. The HHT consists of the combination of EMD and Hilbert spectral analysis. The local energies and the instantaneous frequencies derived from the IMFs through Hilbert transforms can be used to construct an energy-frequency-time distribution, denoted a Hilbert spectrum. The instant paper begins with a discussion of prior approaches to quantification of market volatility, summarizes the HHT method, then describes the application of the method in performing time-frequency analysis of mortgage-market data from the years 1972 through 2000. Filtering by use of the EMD is shown to be useful for quantifying market volatility.

Huang, Norden↗
NASA NTRS2014

Chromospheric Variability: Analysis of 36 years of Time Series from the National Solar Observatory/Sacramento Peak Ca II K-line Monitoring Program

Analysis of more than 36 years of time series of seven parameters measured in the NSO/AFRL/Sac Peak K-line monitoring program elucidates five elucidates five components of the variation: (1) the solar cycle (period approx. 11 years), (2) quasi-periodic variations (periods approx 100 days), (3) a broad band stochastic process (wide range of periods), (4) rotational modulation, and (5) random observational errors. Correlation and power spectrum analyses elucidate periodic and aperiodic variation of the chromospheric parameters. Time-frequency analysis illuminates periodic and quasi periodic signals, details of frequency modulation due to differential rotation, and in particular elucidates the rather complex harmonic structure (1) and (2) at time scales in the range approx 0.1 - 10 years. These results using only full-disk data further suggest that similar analyses will be useful at detecting and characterizing differential rotation in stars from stellar light-curves such as those being produced by NASA's Kepler observatory. Component (3) consists of variations over a range of timescales, in the manner of a 1/f random noise process. A timedependent Wilson-Bappu effect appears to be present in the solar cycle variations (1), but not in the stochastic process (3). Component (4) characterizes differential rotation of the active regions, and (5) is of course not characteristic of solar variability, but the fact that the observational errors are quite small greatly facilitates the analysis of the other components. The recent data suggest that the current cycle is starting late and may be relatively weak. The data analyzed in this paper can be found at the National Solar Observatory web site http://nsosp.nso.edu/cak_mon/, or by file transfer protocol at ftp://ftp.nso.edu/idl/cak.parameters.

Solar physics↗
NASA NTRS2015

Effects of Vehicle Weight and True Versus Indicated Airspeed on BVI Noise During Steady Descending Flight

Blade-vortex interaction noise measurements are analyzed for an AS350B helicopter operated at 7000 ft elevation above sea level. Blade-vortex interaction (BVI) noise from four, 6 degree descent conditions are investigated with descents flown at 80 knot true and indicated airspeed, as well as 4400 and 3915 pound take-off weights. BVI noise is extracted from the acquired acoustic signals by way of a previously developed time-frequency analysis technique. The BVI extraction technique is shown to provide a better localization of BVI noise, compared to the standard Fourier transform integration method. Using this technique, it was discovered that large changes in BVI noise amplitude occurred due to changes in vehicle gross weight. Changes in BVI noise amplitude were too large to be due solely to changes in the vortex strength caused by varying vehicle weight. Instead, it is suggested that vehicle weight modifies the tip-path-plane angle of attack, as well as induced inflow, resulting in large variations of BVI noise. It was also shown that defining flight conditions by true airspeed, rather than indicated airspeed, provides more consistent BVI noise signals.

Stephenson, James H.↗
NASA NTRS2003

Detailed Vibration Analysis of Pinion Gear with Time-Frequency Methods

In this paper, the authors show a detailed analysis of the vibration signal from the destructive testing of a spiral bevel gear and pinion pair containing seeded faults. The vibration signal is analyzed in the time domain, frequency domain and with four time-frequency transforms: the Short Time Frequency Transform (STFT), the Wigner-Ville Distribution with the Choi-Williams kernel (WV-CW), the Continuous Wavelet' Transform (CWT) and the Discrete Wavelet Transform (DWT). Vibration data of bevel gear tooth fatigue cracks, under a variety of operating load levels and damage conditions, are analyzed using these methods. A new metric for automatic anomaly detection is developed and can be produced from any systematic numerical representation of the vibration signals. This new metric reveals indications of gear damage with all of the time-frequency transforms, as well as time and frequency representations, on this data set. Analysis with the CWT detects changes in the signal at low torque levels not found with the other transforms. The WV-CW and CWT use considerably more resources than the STFT and the DWT. More testing of the new metric is needed to determine its value for automatic anomaly detection and to develop fault detection methods for the metric.

Mosher, Marianne↗
NASA NTRS2010

Searching for Turbulent Spectra in Mutliple Spacecraft Observations

The subject of magnetohydrodynamic simulations in the 1980's, magnetic turbulence has reemerged as a favored explanation for the production of fast magnetic reconnection. Recent three-dimensional simulations reveal the formation of secondary islands that serve to shorten the current sheet and increase the accelerating electric field, while both simulations and observations witness electron holes whose collapse energizes electrons. However, little work has been done to investigate the presence of turbulent magnetic and electric field spectra in association with in situ reconnection observations in space. In this study we address this deficiency by analyzing the wavelength spectra of two reconnection events established by multiple spacecraft observations in previous publications. We apply a k-filtering analysis to determine the wave vectors. To optimize the quality of the k-filtering analysis and minimize aliasing, we select multi-spacecraft observations having optimal spacecraft tetrahedral geometry and dimensions that resolve wavelengths within the ion diffusion region. We have chosen one dayside and one nightside reconnection event each observed by the Cluster spacecraft. Spatial correlation between the reconnecting current sheet and a turbulent magnetic spectrum supports an association between turbulence and reconnection in space. The k-filtering method, a wavelet time-frequency spectrum, and a polarization analysis serve to identify wave modes associated with the turbulence. A calculation of the reconnection electric field based on the normal component of the reconnecting magnetic field provides an estimate of the reconnection rate.

Wendel, Deirdre E.↗
NASA NTRS1989

Wavelet analysis introduction and application to radar scattering from water waves

Double-parameter expansions, leading to an instantaneous multiscale analysis, are used to obtain time-frequency graphs from time-series data for interpretation of scatterometer response to sea surface geometry. This method is used to define, as well as possible, one-dimensional signal input (elevation and slope). However, some important fundamental parameters, such as the asymmetric development of short scales in crests and troughs, the long-wave amplitude modulation, the different zero-crossings, and so forth, are easily extracted from the wavelet images.

Chapron, B.↗
NASA NTRS1992

Time-frequency representation of a highly nonstationary signal via the modified Wigner distribution

A new signal analysis technique called the modified Wigner distribution (MWD) is presented. The new signal processing tool has been very successful in determining time frequency representations of highly non-stationary multicomponent signals in both simulations and trials involving actual Space Shuttle Main Engine (SSME) high frequency data. The MWD departs from the classic Wigner distribution (WD) in that it effectively eliminates the cross coupling among positive frequency components in a multiple component signal. This attribute of the MWD, which prevents the generation of 'phantom' spectral peaks, will undoubtedly increase the utility of the WD for real world signal analysis applications which more often than not involve multicomponent signals.

Zoladz, T. F.↗
NASA NTRS2007

Application of the Hilbert-Huang Transform to the Search for Gravitational Waves

We present the application of a novel method of time-series analysis, the Hilbert-Huang Transform, to the search for gravitational waves. This algorithm is adaptive and does not impose a basis set on the data, and thus the time-frequency decomposition it provides is not limited by time-frequency uncertainty spreading. Because of its high time-frequency resolution it has important applications to both signal detection and instrumental characterization. Applications to the data analysis of the ground and space based gravitational wave detectors, LIGO and LISA, are described.

Camp, Jordan B.↗
NASA NTRS2004

Characterizing Oscillatory Bursts in Single-Trial EEG Data

Oscillatory bursts in numerous bands ranging from low (theta) to high frequencies (e.g., gamma) undoubtedly play an important role in cortical dynamics. Largely because of the inadequacy of existing analytic techniques. however, oscillatory bursts and their role in cortical processing remains poorly understood. To study oscillatory bursts effectively one must be able to isolate them and characterize them in the single trial. We describe a series of straightforward analysis techniques that produce useful indices of burst characteristics. First, stimulus-evoked responses are estimated using Differentially Variable Component Analysis (dVCA), and are subtracted from the single-trial. The single-trial characteristics of the evoked responses are stored to identify possible correlations with burst activity. Time-frequency (T-F), or wavelet, analyses are then applied to the single trial residuals. While T-F plots have been used in recent studies to identify and isolate bursts, we go further by fitting each burst in the T-F plot with a two-dimensional Gaussian. This provides a set of burst characteristics, such as, center time. burst duration, center frequency. frequency dispersion. and amplitude, all of which contribute to the accurate characterization of the individual burst. The burst phase can also be estimated. Burst characteristics can be quantified with several standard techniques (e.g.. histogramming and clustering), as well as Bayesian techniques (e.g., blocking) to allow a more parametric description analysis of the characteristics of oscillatory bursts, and the relationships of specific parameters to cortical excitability and stimulus integration.

Knuth, K. H.↗
NASA NTRS1976

High resolution spectral analysis of the Jovian decametric radiation. II - The band-like emissions

Analysis of the bandlike Jovian decametric emission is presented. A model for the active region that accounts for the observed radiation characteristics is described using the measured parameters of the bandlike emission and a model of the Jovian magnetic field. The active region is characterized not only by the fact that an upward-flowing electron stream is caused to radiate in this region, but the stream itself is broken into radiating electron bunches within the active region. Observed undulations of the emission band on the time-frequency plane are interpreted as motions of the active region along a flux tube. The instantaneous location of the active region along the flux tube shows a dependence on the density of the stream entering the active region. The mechanism responsible for density modulation of the stream appears to be common to both the bandlike and simple-S-burst emission types.

Flagg, R. S.↗
NASA NTRS1994

Analysis and prediction of ocean swell using instrumented buoys

During the period 20-23 September 1990, the remnants of Supertyphoon Flo moved into the central North Pacific Ocean with sustained wind speeds of 28 m/s. The strong wind and large fetch area associated with this storm generated long-period swell that propagated to the west coast of North America. National Data Buoy Center moored-buoy stations, located in a network that ranged from the Gulf of Alaska to the California Bight, provided wave spectral estimates of the swell from this storm. The greatest dominant wave periods measured were approximately 20-25 s, and significant wave heights measured ranged from 3 to 8 m. Wave spectra from an array of three nondirectional buoys are used to find the source of the long-period swell. Directional wave spectra from a heave-pitch-roll buoy are also used to make an independent estimate of the source of the swell. The ridge-line method, using time-frequency contour plots of wave spectral energy density, is used to determine the time of swell generation, which is used with the appropriate surface pressure analysis to infer the swell generation area. The diagnosed sources of the swell are also compared with nowcasts from the Global Spectral Ocean Wave Model of the Fleet Numerical Oceanography Center. A simple method of predicting the propagation of ocean swell, by applying a simple kinematic model of wave propagation to the estimated point and time source, is demonstrated.

Mettlach, Theodore↗
NASA NTRS1980

Ultrasonic spectrum analysis using frequency-tracked gated RF pulses

A new method of ultrasonic frequency analysis is introduced which employs frequency-tracked gated RF drive pulses rather than shock-excited broadband spikes to generate the ultrasonic waveform. The new technique, a variation of the sampled-continuous wave technique, eliminates problems associated with finite pulse widths of conventional methods. It is shown to yield correct ultrasonic wave velocity measurements of the sample irrespective of receiver gate width or position provided any portions of two successive echoes are gated simultaneously into the spectrum analyzer. The experimental observations are substantiated by a theoretical model based on the time-frequency domain formulation of ultrasonic frequency analysis.

Cantrell, J. H., Jr.↗
NASA NTRS1994

Wavelets, period-doubling, and time-frequency localization with application to organization of convection over the tropical western Pacific

In this paper, preliminary results in using orthogonal and continuous wavelet transform (WT) to identify period doubling and time-frequency localization in both synthetic and real data are presented. First, the Haar WT is applied to synthetic time series derived from a simple nonlinear dynamical system- a first-order quadratic difference equation. Second, the complex Morlet WT is used to study the time-frequency localization of tropical convection based on a high-resolution Japanese Geostationary Meteorological Satellite infrared (IR) radiance dataset. The Haar WT of the synthetic time series indicates the presence and distinct separation of multiple frequencies in a period-doubling sequence. The period-doubling process generates a multiplicity of intermediate frequencies, which are manifested in the nonuniformity in time with respect to the phase of oscillations in the lower frequencies. Wavelet transform also enables the detection of extremely weak signals in high-order subharmonics resulting from the period-doubling bifurcations. These signals are either undetected or considered statistically insignificant by traditional Fourier analysis. The Morlet WT of the IR radiance dataset indicates the presence of multiple timescales, which are localized in both frequency and time. There are two regimes in the variation of IR radiance, corresponding to the wet and dry periods. Multiple timescales, ranging from semidiurnal, diurnal, synoptic, to intraseasonal with embedding structures, are active in the wet regime. In particular, synoptic variability is more prominent during the wet phase of an intensive intraseasonal cycle. These are not only consistent with, but also show more details than, previous findings by using other techniques. The phase-locking relationships among the oscillations with different time-scales suggest that both synoptic and intraseasonal variations may be mixed oscillations due to the interaction of self-excited oscillations in the tropical atmosphere and external forcings such as annual and diurnal solar radiation variations. Both examples show that WT is a powerful tool for analysis of phenomena involving multiscale interactions that exhibit localization in both frequency and time. A discussion on the caveats in the use of WT in geophysical data analysis is also presented.

Weng, Hengyi↗
NASA NTRS1997

Estimation of Modal Parameters Using a Wavelet-Based Approach

Modal stability parameters are extracted directly from aeroservoelastic flight test data by decomposition of accelerometer response signals into time-frequency atoms. Logarithmic sweeps and sinusoidal pulses are used to generate DAST closed loop excitation data. Novel wavelets constructed to extract modal damping and frequency explicitly from the data are introduced. The so-called Haley and Laplace wavelets are used to track time-varying modal damping and frequency in a matching pursuit algorithm. Estimation of the trend to aeroservoelastic instability is demonstrated successfully from analysis of the DAST data.

Lind, Rick↗
NASA NTRS2001

The Application of Time-Frequency Methods to HUMS

This paper reports the study of four time-frequency transforms applied to vibration signals and presents a new metric for comparing them for fault detection. The four methods to be described and compared are the Short Time Frequency Transform (STFT), the Choi-Williams Distribution (WV-CW), the Continuous Wavelet Transform (CWT) and the Discrete Wavelet Transform (DWT). Vibration data of bevel gear tooth fatigue cracks, under a variety of operating load levels, are analyzed using these methods. The new metric for automatic fault detection is developed and can be produced from any systematic numerical representation of the vibration signals. This new metric reveals indications of gear damage with all of the methods on this data set. Analysis with the CWT detects mechanical problems with the test rig not found with the other transforms. The WV-CW and CWT use considerably more resources than the STFT and the DWT. More testing of the new metric is needed to determine its value for automatic fault detection and to develop methods of setting the threshold for the metric.

Pryor, Anna H.↗
NASA NTRS2004

An on-board processor for a spaceborne Doppler precipitation radar : requirements and preliminary design

Use of Doppler velocity measurement in spaceborne precipitation radar is highly desirable, since it can allow more accurate retrieval of atmospheric latent heating, which depends on both cloud and rain microphysical processes and on dynamical processes, namely, vertical updrafts and downdrafts. However, if the rain within the antenna beam is very non-uniform, biases can result with conventional pulse-pair processing, and more sophisticated processing is required. The solution is a time-frequency approach, which uses the full Doppler spectrum at each point in time. Since precipitation radars normally operate continuously and collect large amounts of data, real-time, on-board processing of data is needed to reduce the data rate. The objective of this work is to develop an on-board data processor for spaceborne Doppler precipitation radar, using field-programmable gate array (FPGA) technology. This paper describes the simulations and analysis that have been completed during the first year of the project. This includes both analytical calculations and bit-true simulations that allow the effects of finite word length to be explored. Tradeoffs have also been considered between speed and space within the chosen FPGA part. The results of these analyses have been used to develop a complete set of requirements for the processor. Initial design work has begun, based on these requirements. The design is being implemented in Verilog, and a description of this early design effort is also be given.

on-board processing↗