Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wavelets”

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 181 records · Page 10

Planetary Transmission Diagnostics

This report presents a methodology for detecting and diagnosing gear faults in the planetary stage of a helicopter transmission. This diagnostic technique is based on the constrained adaptive lifting algorithm. The lifting scheme, developed by Wim Sweldens of Bell Labs, is a time domain, prediction-error realization of the wavelet transform that allows for greater flexibility in the construction of wavelet bases. Classic lifting analyzes a given signal using wavelets derived from a single fundamental basis function. A number of researchers have proposed techniques for adding adaptivity to the lifting scheme, allowing the transform to choose from a set of fundamental bases the basis that best fits the signal. This characteristic is desirable for gear diagnostics as it allows the technique to tailor itself to a specific transmission by selecting a set of wavelets that best represent vibration signals obtained while the gearbox is operating under healthy-state conditions. However, constraints on certain basis characteristics are necessary to enhance the detection of local wave-form changes caused by certain types of gear damage. The proposed methodology analyzes individual tooth-mesh waveforms from a healthy-state gearbox vibration signal that was generated using the vibration separation (synchronous signal-averaging) algorithm. Each waveform is separated into analysis domains using zeros of its slope and curvature. The bases selected in each analysis domain are chosen to minimize the prediction error, and constrained to have the same-sign local slope and curvature as the original signal. The resulting set of bases is used to analyze future-state vibration signals and the lifting prediction error is inspected. The constraints allow the transform to effectively adapt to global amplitude changes, yielding small prediction errors. However, local wave-form changes associated with certain types of gear damage are poorly adapted, causing a significant change in the prediction error. The constrained adaptive lifting diagnostic algorithm is validated using data collected from the University of Maryland Transmission Test Rig and the results are discussed.

Lewicki, David G.↗

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this report is to advance the development of methodologies for scattering by airborne composite vehicles. Although the primary focus continues to be the development of a general purpose computer code for analyzing the entire structure as a single unit, a number of other tasks are also being pursued in parallel with this effort. One of these tasks discussed within is on new finite element formulations and mesh termination schemes. The goal here is to decrease computation time while retaining accuracy and geometric adaptability.The second task focuses on the application of wavelets to electromagnetics. Wavelet transformations are shown to be able to reduce a full matrix to a band matrix, thereby reducing the solutions memory requirements. Included within this document are two separate papers on finite element formulations and wavelets.

Volakis, John L.↗

Rough surface scattering from moving ocean surfaces as an indicator of scattering mechanisms

Ever since Crombie used the Doppler spectrum of RF signals scattered from the ocean surface to provide conclusive evidence for Bragg resonant scattering, scattering from the ocean surface has served as a testing ground for the identification of governing scattering mechanisms in rough surface scattering. The moving ocean surface has the advantage for rough surface scattering that, being a dispersive medium, features with different spatial scales the tendency to travel at different speeds. Thus, by looking at the temporal signature of the scattered field, one can try to isolate the scattering mechanism. This feature has been exploited by various authors. To give a recent example, we cite the work of Plant, who examined the Doppler spectrum of the return signal to infer the presence of two-scale scattering in scattering at medium incidence angles. The problem with restricting oneself to the examination of the Doppler spectrum, as has been done traditionally, is that one restricts the phenomena under investigation to temporally homogeneous phenomena: transient phenomena, which are localized in time, are not localized in frequency space. There is increasing evidence that this type of phenomena, due to specular scattering or breaking waves, may also play an important role in determining the scattering mechanism in ocean-like surfaces. To overcome this problem, we introduce the use of the Wavelet Transform to study the frequency-temporal signatures of the scattered field from moving ocean waves. We calculate this signature using various analytic scattering theories and show that the Wavelet Transform provides a useful tool for separating the different scattering mechanisms operating in scattering in ocean-like surfaces. Next we simulate realistic nonlinear moving ocean waves and calculate the temporal scattered field signature by using the method of moments and the stop-start approximation: the surface is assumed stationary during a scattering event, but moves between incident EM pulses. We compute and compare the Doppler frequency spectrum and the Wavelet Transform of the resulting signal for various types of surfaces and show how different scattering mechanisms may be isolated by this technique. Finally, we compare these numerical results against the analytical predictions in order to obtain the domain of applicability of various scattering theories for moving ocean surfaces which may include wave breaking and specular glints.

Rodriguez, Ernesto↗

Coupling of ELF/ULF energy from lightning and MeV particles to the middle atmosphere, inosphere, and global circuit

In an attempt to explain numerous atmospheric electrical phenomena, the elements of the global electrical circuit are reexamined. In addition to being a 'quasi-static 'DC' generator' and source of radiated energy at VLF and higher, the thunderstorm is found to be a pulse generator, with most of the external energy contained in ELF and ULF pulse currents to the ionosphere (and Earth). The pulse energy is found to deposit largely in the middle atmosphere above the thunderstorm. The VLF and above components are well understood, as are the ULF components due to the conductivity gradient. However, a previously poorly understood ELF component on the millsecond timescale, or 'slow tail,' contains a large fraction of the electrical energy. This component couples strongly to the ionosphere and also launches a unipolar transverse electromagnetic (TEM) wavelet in the radial Earth-ionosphere transmission line. The increase in charge with distance associated with such wavelets, and their ensemble sum at a point, may explain some large mesospheric 'DC' fields but there are still difficulties explaining other than rare occurrences, except for antipodal reconvergence. These millisecond duration unipolar wavelets also coupled to the ionosphere and may trigger other lightning at a distance. A schema is elucidated by which the charge of MeV particles deposited in the middle atmosphere persists for much longer than the local relaxation time. This also gives rise to unipolar waves of global extent which may explain lower-latitude field perturbations associated with solar/geomagnetic events.

Hale, Leslie C.↗

Advanced Signal Processing Methods Applied to Digital Mammography

The work reported here is on the extension of the earlier proposal of the same title, August 1994-June 1996. The report for that work is also being submitted. The work reported there forms the foundation for this work from January 1997 to September 1997. After the earlier work was completed there were a few items that needed to be completed prior to submission of a new and more comprehensive proposal for further research. Those tasks have been completed and two new proposals have been submitted, one to NASA, and one to Health & Human Services WS). The main purpose of this extension was to refine some of the techniques that lead to automatic large scale evaluation of full mammograms. Progress on each of the proposed tasks follows. Task 1: A multiresolution segmentation of background from breast has been developed and tested. The method is based on the different noise characteristics of the two different fields. The breast field has more power in the lower octaves and the off-breast field behaves similar to a wideband process, where more power is in the high frequency octaves. After the two fields are separated by lowpass filtering, a region labeling routine is used to find the largest contiguous region, the breast. Task 2: A wavelet expansion that can decompose the image without zero padding has been developed. The method preserves all properties of the power-of-two wavelet transform and does not add appreciably to computation time or storage. This work is essential for analysis of the full mammogram, as opposed to selecting sections from the full mammogram. Task 3: A clustering method has been developed based on a simple counting mechanism. No ROC analysis has been performed (and was not proposed), so we cannot finally evaluate this work without further support. Task 4: Further testing of the filter reveals that different wavelet bases do yield slightly different qualitative results. We cannot provide quantitative conclusions about this for all possible bases without further support. Task 5: Better modeling does indeed make an improvement in the detection output. After the proposal ended, we came up with some new theoretical explanations that helps in understanding when the D4 filter should be better. This work is currently in the review process. Task 6: N/A. This no longer applies in view of Tasks 4-5. Task 7: Comprehensive plans for further work have been completed. These plans are the subject of two proposals, one to NASA and one to HHS. These proposals represent plans for a complete evaluation of the methods for identifying normal mammograms, augmented with significant further theoretical work.

Stauduhar, Richard P.↗

Data Mining and Optimization Tools for Developing Engine Parameters Tools

This project was awarded for understanding the problem and developing a plan for Data Mining tools for use in designing and implementing an Engine Condition Monitoring System. From the total budget of $5,000, Tricia and I studied the problem domain for developing ail Engine Condition Monitoring system using the sparse and non-standardized datasets to be available through a consortium at NASA Lewis Research Center. We visited NASA three times to discuss additional issues related to dataset which was not made available to us. We discussed and developed a general framework of data mining and optimization tools to extract useful information from sparse and non-standard datasets. These discussions lead to the training of Tricia Erhardt to develop Genetic Algorithm based search programs which were written in C++ and used to demonstrate the capability of GA algorithm in searching an optimal solution in noisy datasets. From the study and discussion with NASA LERC personnel, we then prepared a proposal, which is being submitted to NASA for future work for the development of data mining algorithms for engine conditional monitoring. The proposed set of algorithm uses wavelet processing for creating multi-resolution pyramid of the data for GA based multi-resolution optimal search. Wavelet processing is proposed to create a coarse resolution representation of data providing two advantages in GA based search: 1. We will have less data to begin with to make search sub-spaces. 2. It will have robustness against the noise because at every level of wavelet based decomposition, we will be decomposing the signal into low pass and high pass filters.

Dhawan, Atam P.↗

Abrupt Changes at the Permian/Triassic Boundary: Tempo of Events from High-Resolution Cyclostratigraphy

The Permian/Triassic (P/Tr) boundary (251.4 +/- 3 Myr) is marked by the most severe mass extinction in the geologic record. Recently, precise absolute dating has bracketed the marine extinctions and associated carbon-isotope anomaly within less than 1 Myr. We improve this resolution through high-resolution stratigraphy across the P/Tr boundary in the 331-m Gartnerkofel-1 core and nearby Reppwand outcrop section (Carnic Alps, Austria) utilizing FFT and wavelet timeseries analyses of cyclic components in down-hole core logs of density and natural gamma-ray intensity, and carbon-isotopic ratios of bulk samples. The wavelet analysis indicates continuity of deposition across the P/Tr boundary interval, and the timeseries analyses show evidence for persistent cycles in the ratio of approximately 40: 10: 4.7: 2.3 meters, correlated with Milankovitch-band orbital cycles of approximately 412: 100: 40: 20 kyr (eccentricity 1 and 2, obliquity, and precession), and giving a consistent average sedimentation rate of approximately 10 cm/1,000 yr. Milankovitch periods in delta C-13 and density in these shallow-water carbonates were most likely the result of climatically induced oscillations of sea level and climate, coupled with changes in ocean circulation and productivity, that affected sedimentation. Fluctuations in gamma radiation reflect varying input of clay minerals and the presence of shaly interbeds. Throughout the P/Tr boundary interval in the core, the 100,000-year eccentricity cycle seems to be dominant. Weaker obliquity and precession cycles are in line with the location of the Austrian section in the latest Permian, close to the Equator in the western bight of the Tethys, where obliquity and precessional effects on seasonal contrast might be subdued. Using the improved resolution provided by cycle analysis in the GK-1 core, we find that the dramatic change in the faunal record that marks the P/Tr boundary takes place over less than 6m, or less than 60,000 years. In the nearby Reppwand outcrop section, the same faunal changes occurs over only 0.8 m or about 8,000 years, close to the limit of time-resolution induced by bioturbation and reworking in these sediments. The sharp negative global carbon-isotope shift took place within less than or equal to 40,000 yr, and the isotope excursions persisted for approximately 480,000 yr into the Early Triassic. The results indicate that the severe marine faunal event that marks the P/Tr boundary was very sudden, perhaps less than the resolution window in the GK-1 core, and suggest a catastrophic cause. The wavelet-analysis approach to high-resolution cyclostratigraphy can be applied to other P/Tr boundary sections, and when combined with precise absolute dating and magnetostratigraphic methods promises a significant increase in resolution in determining the correlation and tempo of the end-Permian extinctions and related events worldwide.

Rampino, M. R.↗

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗

Fusing Image Data for Calculating Position of an Object

A computer program has been written for use in maintaining the calibration, with respect to the positions of imaged objects, of a stereoscopic pair of cameras on each of the Mars Explorer Rovers Spirit and Opportunity. The program identifies and locates a known object in the images. The object in question is part of a Moessbauer spectrometer located at the tip of a robot arm, the kinematics of which are known. In the program, the images are processed through a module that extracts edges, combines the edges into line segments, and then derives ellipse centroids from the line segments. The images are also processed by a feature-extraction algorithm that performs a wavelet analysis, then performs a pattern-recognition operation in the wavelet-coefficient space to determine matches to a texture feature measure derived from the horizontal, vertical, and diagonal coefficients. The centroids from the ellipse finder and the wavelet feature matcher are then fused to determine co-location. In the event that a match is found, the centroid (or centroids if multiple matches are present) is reported. If no match is found, the process reports the results of the analyses for further examination by human experts.

Huntsberger, Terrance↗

Estimation of Time-Varying Pilot Model Parameters

Human control behavior is rarely completely stationary over time due to fatigue or loss of attention. In addition, there are many control tasks for which human operators need to adapt their control strategy to vehicle dynamics that vary in time. In previous studies on the identification of time-varying pilot control behavior wavelets were used to estimate the time-varying frequency response functions. However, the estimation of time-varying pilot model parameters was not considered. Estimating these parameters can be a valuable tool for the quantification of different aspects of human time-varying manual control. This paper presents two methods for the estimation of time-varying pilot model parameters, a two-step method using wavelets and a windowed maximum likelihood estimation method. The methods are evaluated using simulations of a closed-loop control task with time-varying pilot equalization and vehicle dynamics. Simulations are performed with and without remnant. Both methods give accurate results when no pilot remnant is present. The wavelet transform is very sensitive to measurement noise, resulting in inaccurate parameter estimates when considerable pilot remnant is present. Maximum likelihood estimation is less sensitive to pilot remnant, but cannot detect fast changes in pilot control behavior.

Zaal, Peter M. T.↗

Shearlet Features for Registration of Remotely Sensed Multitemporal Images

We investigate the role of anisotropic feature extraction methods for automatic image registration of remotely sensed multitemporal images. Building on the classical use of wavelets in image registration, we develop an algorithm based on shearlets, a mathematical generalization of wavelets that offers increased directional sensitivity. Initial experimental results on LANDSAT images are presented, which indicate superior performance of the shearlet algorithm when compared to classical wavelet algorithms.

Science data processing↗

Shearlet Features for Registration of Remotely Sensed Multitemporal Images

We investigate the role of anisotropic feature extraction methods for automatic image registration of remotely sensed multitemporal images. Building on the classical use of wavelets in image registration, we develop an algorithm based on shearlets, a mathematical generalization of wavelets that offers increased directional sensitivity. Experimental results on multitemporal Landsat images are presented, which indicate superior performance of the shearlet algorithm when compared to classical wavelet algorithms.

Science Data Processing↗

Shearlet Features for Registration of Remotely Sensed Multitemporal Images

We investigate the role of anisotropic feature extraction methods for automatic image registration of remotely sensed multitemporal images. Building on the classical use of wavelets in image registration, we develop an algorithm based on shearlets, a mathematical generalization of wavelets that offers increased directional sensitivity. Experimental results on multitemporal Landsat images are presented, which indicate superior performance of the shearlet algorithm when compared to classical wavelet algorithms.

Science Data Processing↗

Absolute value optimization to estimate phase properties of stochastic time series

Most existing deconvolution techniques are incapable of determining phase properties of wavelets from time series data; to assure a unique solution, minimum phase is usually assumed. It is demonstrated, for moving average processes of order one, that deconvolution filtering using the absolute value norm provides an estimate of the wavelet shape that has the correct phase character when the random driving process is nonnormal. Numerical tests show that this result probably applies to more general processes.

Scargle, J. D.↗

The magnetospheric response to 8-minute period strong-amplitude upstream pressure variations

This paper documents a series of brief, strong (delta p/p = 1), dynamic pressure oscillations that occurred in the region upstream of the earth's bow shock during a period of radial interplanetary magnetic field. The analyzed set of oscillations, which may be either intrinsic solar wind or bow shock-related phenomena, recur approximately every 8-10 min, and their magnetic field signatures occur nearly simultaneously over great distances transverse to the earth-sun line. The pressure oscillations appear to drive tailward-moving magnetopause surface wavelets. In turn, the surface wavelets can be identified as hydromagnetic waves with strong compressional components in the outer magnetosphere and as quasi-periodic variations in electron precipitation and high-latitude ground pulsations. Observations by spacecraft in the outer dayside magnetosphere are used to predict geosynchronous and subsolar magnetic field strengths, the location of the subsolar magnetopause, the solar wind dynamic pressure, and variations in the energetic magnetospheric ion flux.

Sibeck, D. G.↗

Shiftable multiscale transforms

A type of translational invariance, referred to as shiftability, is defined for wavelet transforms. The property of shiftability is first discussed with respect to individual parameters: spatial position, orientation, and scale. The discussion then focuses on transformations that are simultaneously shiftable with respect to subsets of these parameters. It is shown that the critical sampling condition on the wavelet transforms must be relaxed to achieve shiftability. Two example transforms are implemented and applied to several signal and image processing problems.

Simoncelli, Eero P.↗

A new approach to the photon localization problem

Since wavelets form a representation of the Poincare group, it is possible to construct a localized superposition of light waves with different frequencies in a Lorentz-covariant manner. This localized wavelet satisfies a Lorentz-invariant uncertainty relation, and also the Lorentz-invariant Parseval's relation. A quantitative analysis is given for the difference between photons and localized waves. It is then shown that this localized entity corresponds to a relativistic photon with a sharply defined momentum in the non-localization limit. Waves are not particles. It is confirmed that the wave-particle duality is subject to the uncertainty principle.

Han, D.↗

Detection of signature consistent with cosmological time dilation in gamma-ray bursts

If gamma ray bursts are at cosmological distances-as suggested by their isotropic distribution on the sky and by their number-intensity relation-then the burst profiles will be stretched in time, by an amount proportional to the redshift, 1 + Z. We have tested data from the Compton Gamma Ray Observatory's (CGRO's) Burst and Transient Source Experiment (BATSE) for such time dilation. Out of 590 bursts observed by BATSE, 131 bursts were analyzed; bursts with durations shorter than 1.5 s were excluded. We used three tests to compare the timescales of bright and dim bursts, the latter, on average, being more distant than the former. Our measures of timescale are constructed to avoid selection effects arising from intensity differences by rescaling all bursts to fiducial levels of peak intensity and noise bias. (1) We found that the total rescaled count above background for the dim burst ensemble is approximately twice that for the brightest bursts-translating into longer durations for the dim bursts. (2) Wavelet-transform decompositions of the burst profiles confirmed that this dilation operates over a broad range of timescales. (3) Structure on the shortest timescales was examined using a procedure which aligns the highest peaks of profiles from which the noise has been optimally removed using a wavelet threshold technique. In all three tests, the dim bursts are stretched by a factor of approximately 2 relative to the bright ones, over seven octaves of timescale. We calibrated the measurements by dilating synthetic bursts that approximate the temporal characteristics of bright BATSE bursts. Results are consistent with bursts at BATSE's peak-flux completeness limit being at cosmological distances corresponding to Z approximately equal to 1, and thus with independent cosmological interpretations of the BATSE number-intensity relation. Alternative explanations of our results, arising from the nature of physical processes in bursts, are still possible.

Norris, J. P.↗