Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FFT”

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 127 records · Page 7

Frequency Analyzed Laser Light Scattering (Falls)

A simple and robust light scattering technique for spray diagnostics was developed. The technique provides a spray's intact length and breakup frequency. Scattered light is collected from a laser sheet that illuminates a portion of the spray. The scattered light is imaged on a fast photodiode whose output is sampled at 20 kHz. Fast Fourier transform (FFT) analysis is performed on the sampled signal. The shape of the FFT produced can be correlated to one of three distinct spray regimes: intact sheet, ligaments, or droplets. Strobe and laser illuminated photographs are also shown for reference.

Rhys, Noah O.↗

Ultrasonic Resonance Spectroscopy of Composite Rims for Flywheel Rotors

Flywheel energy storage devices comprising multilayered composite rotor systems are being studied extensively for utilization in the International Space Station. These composite material systems were investigated with a recently developed ultrasonic resonance spectroscopy technique. The ultrasonic system employs a continuous swept-sine waveform and performs a fast Fourier transform (FFT) on the frequency response spectrum. In addition, the system is capable of equalizing the amount of energy at each frequency. Equalization of the frequency spectrum, along with interpretation of the second FFT, aids in the evaluation of the fundamental frequency. The frequency responses from multilayered material samples, with and without known defects, were analyzed to assess the capabilities and limitations of this nondestructive evaluation technique for material characterization and defect detection. Amplitude and frequency changes were studied from ultrasonic responses of thick composite rings and a multiring composite rim. A composite ring varying in thickness was evaluated to investigate the full thickness resonance. The frequency response characteristics from naturally occurring voids in a composite ring were investigated. Ultrasonic responses were compared from regions with and without machined voids in a composite ring and a multiring composite rim. Finally, ultrasonic responses from the multiring composite rim were compared before and after proof spin testing to 63,000 rpm.

Harmon, Laura M.↗

Image Processing Software

The Ames digital image velocimetry technology has been incorporated in a commercially available image processing software package that allows motion measurement of images on a PC alone. The software, manufactured by Werner Frei Associates, is IMAGELAB FFT. IMAGELAB FFT is a general purpose image processing system with a variety of other applications, among them image enhancement of fingerprints and use by banks and law enforcement agencies for analysis of videos run during robberies.

Source record↗

Long-range correlation properties of coding and noncoding DNA sequences: GenBank analysis

An open question in computational molecular biology is whether long-range correlations are present in both coding and noncoding DNA or only in the latter. To answer this question, we consider all 33301 coding and all 29453 noncoding eukaryotic sequences--each of length larger than 512 base pairs (bp)--in the present release of the GenBank to dtermine whether there is any statistically significant distinction in their long-range correlation properties. Standard fast Fourier transform (FFT) analysis indicates that coding sequences have practically no correlations in the range from 10 bp to 100 bp (spectral exponent beta=0.00 +/- 0.04, where the uncertainty is two standard deviations). In contrast, for noncoding sequences, the average value of the spectral exponent beta is positive (0.16 +/- 0.05) which unambiguously shows the presence of long-range correlations. We also separately analyze the 874 coding and the 1157 noncoding sequences that have more than 4096 bp and find a larger region of power-law behavior. We calculate the probability that these two data sets (coding and noncoding) were drawn from the same distribution and we find that it is less than 10(-10). We obtain independent confirmation of these findings using the method of detrended fluctuation analysis (DFA), which is designed to treat sequences with statistical heterogeneity, such as DNA's known mosaic structure ("patchiness") arising from the nonstationarity of nucleotide concentration. The near-perfect agreement between the two independent analysis methods, FFT and DFA, increases the confidence in the reliability of our conclusion.

Non-NASA Center↗

Pre-Hardware Optimization and Implementation Of Fast Optics Closed Control Loop Algorithms

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). The FFT is particularly useful in two-dimensional (2-D) image processing (FFT2) within optical systems control. However, timing constraints of a fast optics closed control loop would require a supercomputer to run the software implementation of the FFT2 and its inverse, as well as other image processing representative algorithm, such as numerical image folding and fringe feature extraction. A laboratory supercomputer is not always available even for ground operations and is not feasible for a night project. However, the computationally intensive algorithms still warrant alternative implementation using reconfigurable computing technologies (RC) such as Digital Signal Processors (DSP) and Field Programmable Gate Arrays (FPGA), which provide low cost compact super-computing capabilities. We present a new RC hardware implementation and utilization architecture that significantly reduces the computational complexity of a few basic image-processing algorithm, such as FFT2, image folding and phase diversity for the NASA Solar Viewing Interferometer Prototype (SVIP) using a cluster of DSPs and FPGAs. The DSP cluster utilization architecture also assures avoidance of a single point of failure, while using commercially available hardware. This, combined with the control algorithms pre-hardware optimization, or the first time allows construction of image-based 800 Hertz (Hz) optics closed control loops on-board a spacecraft, based on the SVIP ground instrument. That spacecraft is the proposed Earth Atmosphere Solar Occultation Imager (EASI) to study greenhouse gases CO2, C2H, H2O, O3, O2, N2O from Lagrange-2 point in space. This paper provides an advanced insight into a new type of science capabilities for future space exploration missions based on on-board image processing for control and for robotics missions using vision sensors. It presents a top-level description of technologies required for the design and construction of SVIP and EASI and to advance the spatial-spectral imaging and large-scale space interferometry science and engineering.

Kizhner, Semion↗

Use of Modal Acoustic Emission to Monitor Damage Progression in Carbon Fiber/Epoxy Tows and Implications for Composite Structures

This slide presentation reviews the use of Modal Acoustic Emission to monitor damage progression to carbon fiber/epoxy tows. There is a risk for catastrophic failure of composite overwrapped pressure vessels (COPVs) due to burst-before-leak (BBL) stress rupture (SR) failure of carbon-epoxy (C/Ep) COPVs. A lack of quantitative nondestructive evaluation (NDE) is causing problems in current and future spacecraft designs. It is therefore important to develop and demonstrate critical NDE that can be implemented during stages of the design process since the observed rupture can occur with little of no advanced warning. Therefore a program was required to develop quantitative acoustic emission (AE) procedures specific to C/Ep overwraps, but which also have utility for monitoring damage accumulation in composite structure in general, and to lay the groundwork for establishing critical thresholds for accumulated damage in composite structures, such as COPVs, so that precautionary or preemptive engineering steps can be implemented to minimize of obviate the risk of catastrophic failure. A computed Felicity Ratio (FR) coupled with fast Fourier Transform (FFT) frequency analysis shows promise as an analytical pass/fail criterion. The FR analysis and waveform and FFT analysis are reviewed

Waller, Jess M.↗

Use of Modal Acoustic Emission to Monitor Damage Progression in Carbon Fiber/Epoxy and Implications for Composite Structures

Broad-band modal acoustic emission (AE) data was used to characterize micromechanical damage progression in uniaxial IM7 and T1000 carbon fiber-epoxy tows and an IM7 composite overwrapped pressure vessel (COPV) subjected to an intermittent load hold tensile stress profile known to activate the Felicity ratio (FR). Damage progression was followed by inspecting the Fast Fourier Transforms (FFTs) associated with acoustic emission events. FFT analysis revealed the occurrence of cooperative micromechanical damage events in a frequency range between 100 kHz and 1 MHz. Evidence was found for the existence of a universal damage parameter, referred to here as the critical Felicity ratio, or Felicity ratio at rupture (FR*), which had a value close to 0.96 for the tows and the COPV tested. The implications of using FR* to predict failure in carbon/epoxy composite materials and related composite components such as COPVs are discussed. Trends in the FFT data are also discussed; namely, the difference between the low and high energy events, the difference between early and late-life events, comparison of IM7 and T1000 damage progression, and lastly, the similarity of events occurring at the onset of significant acoustic emission used to calculate the FR.

Waller, J. M.↗

Optimal Padding for the Two-Dimensional Fast Fourier Transform

One-dimensional Fast Fourier Transform (FFT) operations work fastest on grids whose size is divisible by a power of two. Because of this, padding grids (that are not already sized to a power of two) so that their size is the next highest power of two can speed up operations. While this works well for one-dimensional grids, it does not work well for two-dimensional grids. For a two-dimensional grid, there are certain pad sizes that work better than others. Therefore, the need exists to generalize a strategy for determining optimal pad sizes. There are three steps in the FFT algorithm. The first is to perform a one-dimensional transform on each row in the grid. The second step is to transpose the resulting matrix. The third step is to perform a one-dimensional transform on each row in the resulting grid. Steps one and three both benefit from padding the row to the next highest power of two, but the second step needs a novel approach. An algorithm was developed that struck a balance between optimizing the grid pad size with prime factors that are small (which are optimal for one-dimensional operations), and with prime factors that are large (which are optimal for two-dimensional operations). This algorithm optimizes based on average run times, and is not fine-tuned for any specific application. It increases the amount of times that processor-requested data is found in the set-associative processor cache. Cache retrievals are 4-10 times faster than conventional memory retrievals. The tested implementation of the algorithm resulted in faster execution times on all platforms tested, but with varying sized grids. This is because various computer architectures process commands differently. The test grid was 512 512. Using a 540 540 grid on a Pentium V processor, the code ran 30 percent faster. On a PowerPC, a 256x256 grid worked best. A Core2Duo computer preferred either a 1040x1040 (15 percent faster) or a 1008x1008 (30 percent faster) grid. There are many industries that can benefit from this algorithm, including optics, image-processing, signal-processing, and engineering applications.

Dean, Bruce H.↗

Technique for Solving Electrically Small to Large Structures for Broadband Applications

Fast iterative algorithms are often used for solving Method of Moments (MoM) systems, having a large number of unknowns, to determine current distribution and other parameters. The most commonly used fast methods include the fast multipole method (FMM), the precorrected fast Fourier transform (PFFT), and low-rank QR compression methods. These methods reduce the O(N) memory and time requirements to O(N log N) by compressing the dense MoM system so as to exploit the physics of Green s Function interactions. FFT-based techniques for solving such problems are efficient for spacefilling and uniform structures, but their performance substantially degrades for non-uniformly distributed structures due to the inherent need to employ a uniform global grid. FMM or QR techniques are better suited than FFT techniques; however, neither the FMM nor the QR technique can be used at all frequencies. This method has been developed to efficiently solve for a desired parameter of a system or device that can include both electrically large FMM elements, and electrically small QR elements. The system or device is set up as an oct-tree structure that can include regions of both the FMM type and the QR type. The system is enclosed with a cube at a 0- th level, splitting the cube at the 0-th level into eight child cubes. This forms cubes at a 1st level, recursively repeating the splitting process for cubes at successive levels until a desired number of levels is created. For each cube that is thus formed, neighbor lists and interaction lists are maintained. An iterative solver is then used to determine a first matrix vector product for any electrically large elements as well as a second matrix vector product for any electrically small elements that are included in the structure. These matrix vector products for the electrically large and small elements are combined, and a net delta for a combination of the matrix vector products is determined. The iteration continues until a net delta is obtained that is within the predefined limits. The matrix vector products that were last obtained are used to solve for the desired parameter. The solution for the desired parameter is then presented to a user in a tangible form; for example, on a display.

Jandhyala, Vikram↗

Method for Pre-Conditioning a Measured Surface Height Map for Model Validation

This software allows one to up-sample or down-sample a measured surface map for model validation, not only without introducing any re-sampling errors, but also eliminating the existing measurement noise and measurement errors. Because the re-sampling of a surface map is accomplished based on the analytical expressions of Zernike-polynomials and a power spectral density model, such re-sampling does not introduce any aliasing and interpolation errors as is done by the conventional interpolation and FFT-based (fast-Fourier-transform-based) spatial-filtering method. Also, this new method automatically eliminates the measurement noise and other measurement errors such as artificial discontinuity. The developmental cycle of an optical system, such as a space telescope, includes, but is not limited to, the following two steps: (1) deriving requirements or specs on the optical quality of individual optics before they are fabricated through optical modeling and simulations, and (2) validating the optical model using the measured surface height maps after all optics are fabricated. There are a number of computational issues related to model validation, one of which is the "pre-conditioning" or pre-processing of the measured surface maps before using them in a model validation software tool. This software addresses the following issues: (1) up- or down-sampling a measured surface map to match it with the gridded data format of a model validation tool, and (2) eliminating the surface measurement noise or measurement errors such that the resulted surface height map is continuous or smoothly-varying. So far, the preferred method used for re-sampling a surface map is two-dimensional interpolation. The main problem of this method is that the same pixel can take different values when the method of interpolation is changed among the different methods such as the "nearest," "linear," "cubic," and "spline" fitting in Matlab. The conventional, FFT-based spatial filtering method used to eliminate the surface measurement noise or measurement errors can also suffer from aliasing effects. During re-sampling of a surface map, this software preserves the low spatial-frequency characteristic of a given surface map through the use of Zernike-polynomial fit coefficients, and maintains mid- and high-spatial-frequency characteristics of the given surface map by the use of a PSD model derived from the two-dimensional PSD data of the mid- and high-spatial-frequency components of the original surface map. Because this new method creates the new surface map in the desired sampling format from analytical expressions only, it does not encounter any aliasing effects and does not cause any discontinuity in the resultant surface map.

Sidick, Erkin↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗

Comparison of Type 2 versus Type 3 Carrier Tracking Loops under High Dynamic Signal Conditions

Virtually all deep-space missions today and in the past employed Type 2 tracking loops in their spacecraft receivers. This was acceptable since the uplink from a single ground antenna makes use of very accurate trajectory information where the radiated frequency is Doppler compensated. However, in the desire to conserve ground resources, the technique of Multiple Uplinks per Antenna (MUPA) is being considered when using a single ground antenna to uplink to a constellation of spacecraft within its beam. This single frequency technique requires the use of non-standard signal acquisition techniques since the constellation may realize different orbits and trajectories around a planetary body leading to varying Doppler offsets and dynamics at each member spacecraft. In such a case, a Type 2 tracking loop may not suffice. An earlier study discussed the various methods for acquiring signals in spacecraft receivers with uncompensated Doppler at the uplink station. In the earlier study, three methods of achieving acquisition and lock were discussed for signals with large Doppler frequency offsets as well as high dynamics; 1) a step-and-sweep signal search algorithm, 2) an FFT signal search algorithm, and, 3) on-board tuning of the signal using uploaded trajectory information or lookup tables. On-board tuning appeared a feasible solution given the trend towards utilizing software-defined radios (SDRs). In the absence of on-board tuning, a search algorithm (e.g., step-and-sweep or FFT) could be utilized in the SDR. One such consideration is in the design of the tracking loop filter, where a Type 3 loop may be more amenable to handling large uncompensated Doppler effects. This report presents an analysis of comparing performance of Type 2 versus Type 3 carrier tracking loops when acquiring and tracking signals with large Doppler dynamics.

Abraham, Douglas↗

Digital processing of radiographic images

Some techniques are presented and the software documentation for the digital enhancement of radiographs. Both image handling and image processing operations are considered. The image handling operations dealt with are: (1) conversion of format of data from packed to unpacked and vice versa; (2) automatic extraction of image data arrays; (3) transposition and 90 deg rotations of large data arrays; (4) translation of data arrays for registration; and (5) reduction of the dimensions of data arrays by integral factors. Both the frequency and the spatial domain approaches are presented for the design and implementation of the image processing operation. It is shown that spatial domain recursive implementation of filters is much faster than nonrecursive implementations using fast fourier transforms (FFT) for the cases of interest in this work. The recursive implementation of a class of matched filters for enhancing image signal to noise ratio is described. Test patterns are used to illustrate the filtering operations. The application of the techniques to radiographic images of metallic structures is demonstrated through several examples.

Bond, A. D.↗

Power spectral density estimation by spline smoothing in the frequency domain

An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying FFT techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.

Defigueiredo, R. J. P.↗

Induced shock pulse testing by transient waveform control

A method of synthesizing an arbitrarily shaped transient time pulse on vibration exciters is described. The transient waveform control technique is based on recent developments in digital time series analysis, the real-time FFT processor. A brief description of the theory, error estimates and hardware/software implementation to the JPL Dynamic Environmental Testing Laboratory is presented.

Kim, B. K.↗

Image analysis techniques associated with automatic data base generation.

This paper considers some basic problems relating to automatic data base generation from imagery, the primary emphasis being on fast and efficient automatic extraction of relevant pictorial information. Among the techniques discussed are recursive implementations of some particular types of filters which are much faster than FFT implementations, a 'sequential similarity detection' technique of implementing matched filters, and sequential linear classification of multispectral imagery. Several applications of the above techniques are presented including enhancement of underwater, aerial and radiographic imagery, detection and reconstruction of particular types of features in images, automatic picture registration and classification of multiband aerial photographs to generate thematic land use maps.

Bond, A. D.↗

Numerical comparisons of a modified Rayleigh approach with other rough surface EM scattering solutions.

A modification of the classical Rayleigh technique utilizing the fast Fourier transform, developed into the so-called Rayleigh-FFT method for electromagnetic scattering, is shown to obtain valid scattering results from arbitrary irregular periodic structures composed of a rough layer over a rough half-space. The considerable usefulness of the Rayleigh approach is thus established, and Beckmann's (1968) contention (that the convergence of the method is reasonably rapid only if the roughness is very slight) is refuted.

Jiracek, G. R.↗

Application of Fourier transforms for microwave radiometric inversions

Existing microwave radiometer technology now provides a suitable method for remote determination of the ocean surface's absolute brightness temperature. To extract the brightness temperature of the water from the antenna temperature, an unstable Fredholm integral equation of the first kind is solved. Fourier transform techniques are used to invert the integral after it is placed into a cross correlation form. Application and verification of the methods to a two-dimensional modeling of a laboratory wave tank system are included. The instability of the ill-posed Fredholm equation is examined and a restoration procedure is included which smooths the resulting oscillations. With the recent availability and advances of fast Fourier transform (FFT) techniques, the method presented becomes very attractive in the evaluation of large quantities of data.

Holmes, J. J.↗