Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast Fourier transform (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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Adaptive line enhancers for fast acquisition

Three adaptive line enhancer (ALE) algorithms and architectures - namely, conventional ALE, ALE with double filtering, and ALE with coherent accumulation - are investigated for fast carrier acquisition in the time domain. The advantages of these algorithms are their simplicity, flexibility, robustness, and applicability to general situations including the Earth-to-space uplink carrier acquisition and tracking of the spacecraft. In the acquisition mode, these algorithms act as bandpass filters; hence, the carrier-to-noise ratio (CNR) is improved for fast acquisition. In the tracking mode, these algorithms simply act as lowpass filters to improve signal-to-noise ratio; hence, better tracking performance is obtained. It is not necessary to have a priori knowledge of the received signal parameters, such as CNR, Doppler, and carrier sweeping rate. The implementation of these algorithms is in the time domain (as opposed to the frequency domain, such as the fast Fourier transform (FFT)). The carrier frequency estimation can be updated in real time at each time sample (as opposed to the batch processing of the FFT). The carrier frequency to be acquired can be time varying, and the noise can be non-Gaussian, nonstationary, and colored.

Yeh, H.-G.↗
NASA NTRS1994

A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth-order central differences through fast Fourier transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large scale features, such as the total circulation around the roll-up region, are adequately resolved.

Weinan, E.↗
NASA NTRS1987

Measurement of human pilot dynamic characteristics in flight simulation

Fast Fourier Transform (FFT) and Least Square Error (LSE) estimation techniques were applied to the problem of identifying pilot-vehicle dynamic characteristics in flight simulation. A brief investigation of the effects of noise, input bandwidth and system delay upon the FFT and LSE techniques was undertaken using synthetic data. Data from a piloted simulation conducted at NASA Ames Research Center was then analyzed. The simulation was performed in the NASA Ames Research Center Variable Stability CH-47B helicopter operating in fixed-basis simulator mode. The piloting task consisted of maintaining the simulated vehicle over a moving hover pad whose motion was described by a random-appearing sum of sinusoids. The two test subjects used a head-down, color cathode ray tube (CRT) display for guidance and control information. Test configurations differed in the number of axes being controlled by the pilot (longitudinal only versus longitudinal and lateral), and in the presence or absence of an important display indicator called an 'acceleration ball'. A number of different pilot-vehicle transfer functions were measured, and where appropriate, qualitatively compared with theoretical pilot- vehicle models. Some indirect evidence suggesting pursuit behavior on the part of the test subjects is discussed.

Reedy, James T.↗
NASA NTRS1997

Registration assisted mosaic generation

This paper presents a general strategy for assembling mosaics from numerous individual images where uncertainty exists in the position and orientation of those images. Both of the presented applications relate to remotely operated camera platforms, the first being the Galileo solid state imaging (SSI) camera presently in orbit around Jupiter, and the second being the Imager for Mars Pathfinder (IMP) stereo camera on Mars. A basic strategy in both applications is to determine the correct relative camera pointing followed by direct map projection of the images. It is assumed that approximate camera pointing exists sufficient to locate adjacent images and to place initial tiepoints within reach of the correlator. Spatial correlation is used to fix tiepoints whose initial locations are predicted by the camera pointing. We use either an fast fourier transform (fft) algorithm or a variant of Gruen's scheme permitting limited image rotation and skew. The Gruen correlator has three hierarchical modes: 1) A classical spatial least squares correlation on integral pixel boundaries used when rotation is small. 2) An annealing non-deterministic search used when rotations are unknown. A simplex deterministic search used for the end game. The correlation operation can be performed either interactively or autonomously. The final camera pointing solution relies upon a simplex downhill search in 2n or 3n dimensions where n is the number of images comprising the mosaic and the objective function to be minimized is the disagreement between tiepoint locations predicted from the camera pointing with those observed by the correlator. For Galileo the 3n unknowns are euler angles defining camera pointing in planet coordinates, and for Mars Pathfinder they are 2n unknowns representing commanded azimuth and elevation in the Lander coordinate system.

Lorre, Jean J.↗
NASA NTRS1996

Fast Nitrogen Atoms from Dissociative Excitation of N2 by Electron Impact

The Doppler profiles of one of the fine structure lines of the N I (1200 A) g (sup 4)S(sup 0)-(sup 4)P multiplet and of the N II (1085 A) g (sup 3)p(sup O)-(sup 3)D multiplet have been measured. Excitation of the multiplets is produced by electron impact dissociative excitation of N2. The experimental line profiles are evaluated by fast Fourier transform (FFT) techniques and analysis of the profiles yields the kinetic energy distribution of fragments. The full width at half maximum (FWHM) of N I (1200 A) increases from 27+/-6 mA at 30 eV to 37+/-4 mA at 100 eV as the emission cross section of the dissociative ionization excitation process becomes more important relative to the dissociative excitation process. The FWHM of the N II (1085 A) line is 36+/-4 mA at 100 eV. For each multiplet the kinetic energy distribution function of each of the two fragment N atoms (ions) is much broader than thermal with a mean energy above 1.0 eV. The dissociation process with the largest cross section is predissociation and predominantly produces N atoms with kinetic energy distributions having mean energies above 0.5 eV. Dissociative processes can lead to a substantial escape flux of N I atoms from the satellites, Titan and Triton of the outer planets.

Ajello, Joseph M.↗
NASA NTRS1999

Technical Report Series on Global Modeling and Data Assimilation: Filtering Techniques on a Stretched Grid General Circulation Model - Volume 16

This report documents the techniques used to filter quantities on a stretched grid general circulation model. Standard high-latitude filtering techniques (e.g., using an FFT (Fast Fourier Transformations) to decompose and filter unstable harmonics at selected latitudes) applied on a stretched grid are shown to produce significant distortions of the prognostic state when used to control instabilities near the pole. A new filtering technique is developed which accurately accounts for the non-uniform grid by computing the eigenvectors and eigenfrequencies associated with the stretching. A filter function, constructed to selectively damp those modes whose associated eigenfrequencies exceed some critical value, is used to construct a set of grid-spaced weights which are shown to effectively filter without distortion. Both offline and GCM (General Circulation Model) experiments are shown using the new filtering technique. Finally, a brief examination is also made on the impact of applying the Shapiro filter on the stretched grid.

Takacs, Lawrence L.↗
NASA NTRS1997

Tracking of Ice Edges and Ice Floes by Wavelet Analysis of SAR Images

This paper demonstrates the use of wavelet transforms in the tracking of sequential ice features in the ERS-1 synthetic aperture radar (SAR) imagery, especially in situations where feature correlation techniques fail to yield reasonable results. Examples include the evolution of the St. Lawrence polynya and summer sea ice change in the Beaufort Sea. For the polynya, the evolution of the region of young ice growth surrounding a polynya can be easily tracked by wavelet analysis due to the large backscatter difference between the young and old ice. Also within the polynya, a 2D fast Fourier transform (FFT) is used to identify the extent of the Langmuir circulation region, which is coincident with the wave-agitated frazil ice growth region, where the sea ice experiences its fastest growth. Therefore, the combination of wavelet and FFT analysis of SAR images provides for the large-scale monitoring of different polynya features. For summer ice, previous work shows that this is the most difficult period for ice trackers due to the lack of features on the sea ice cover. The multiscale wavelet analysis shows that this method delineates the detailed floe shapes during this period, so that between consecutive images, the floe translation and rotation can be estimated.

Liu, Antony K.↗
NASA NTRS1999

Measuring Directional Wave Spectra and Wind Speed with a Scanning Radar Altimeter

The geometry for the NASA Scanning Radar Altimeter (SRA) is shown. It transmits a 8-ns duration pulse at Ka-band (8.3 mm) and measures time of flight as it scans a 1 degree (two-way) beam from left to right across the aircraft ground track. The most recent configuration determines the surface elevation at 64 points spaced at uniform angular intervals of about 0.7 across a swath whose width is about 0.8 times the aircraft altitude. The system generates these raster lines of the surface topography beneath the aircraft at about a 10 Hz rate. In postflight processing the SRA wave topographic data are transformed with a two-dimensional Fast Fourier Transformation (FFT) and Doppler corrected to produce directional wave spectra. The SRA is not absolutely calibrated in power, but by measuring the relative fall-off of backscatter with increasing incidence angle, the SRA can also determine the mean square slope (mss) of the sea surface, a surrogate for wind speed. For the slope-dependent specular point model of radar sea surface scattering, an expression approximated by a geometric optics form, for the relative variation with incidence angle of the normalized backscatter radar cross section would be sigma (sup 0) (sub rel) = sec (exp 4) theta exp (-tan squared theta/mss) where theta is the off-nadir incidence angle.

Walsh, E. J.↗
NASA NTRS2000

Advanced Engine Health Management Applications of the SSME Real-Time Vibration Monitoring System

The Real Time Vibration Monitoring System (RTVMS) is a 32-channel high speed vibration data acquisition and processing system developed at Marshall Space Flight Center (MSFC). It Delivers sample rates as high as 51,200 samples/second per channel and performs Fast Fourier Transform (FFT) processing via on-board digital signal processing (DSP) chips in a real-time format. Advanced engine health assessment is achieved by utilizing the vibration spectra to provide accurate sensor validation and enhanced engine vibration redlines. Discrete spectral signatures (such as synchronous) that are indicators of imminent failure can be assessed and utilized to mitigate catastrophic engine failures- a first in rocket engine health assessment. This paper is presented in viewgraph form.

Fiorucci, Tony R.↗
NASA NTRS2001

Thermodynamic Behavior of Nano-sized Gold Clusters on the (001) Surface

We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Paik, Sun M.↗
NASA NTRS1994

Applications Performance Under MPL and MPI on NAS IBM SP2

On July 5, 1994, an IBM Scalable POWER parallel System (IBM SP2) with 64 nodes, was installed at the Numerical Aerodynamic Simulation (NAS) Facility Each node of NAS IBM SP2 is a "wide node" consisting of a RISC 6000/590 workstation module with a clock of 66.5 MHz which can perform four floating point operations per clock with a peak performance of 266 Mflop/s. By the end of 1994, 64 nodes of IBM SP2 will be upgraded to 160 nodes with a peak performance of 42.5 Gflop/s. An overview of the IBM SP2 hardware is presented. The basic understanding of architectural details of RS 6000/590 will help application scientists the porting, optimizing, and tuning of codes from other machines such as the CRAY C90 and the Paragon to the NAS SP2. Optimization techniques such as quad-word loading, effective utilization of two floating point units, and data cache optimization of RS 6000/590 is illustrated, with examples giving performance gains at each optimization step. The conversion of codes using Intel's message passing library NX to codes using native Message Passing Library (MPL) and the Message Passing Interface (NMI) library available on the IBM SP2 is illustrated. In particular, we will present the performance of Fast Fourier Transform (FFT) kernel from NAS Parallel Benchmarks (NPB) under MPL and MPI. We have also optimized some of Fortran BLAS 2 and BLAS 3 routines, e.g., the optimized Fortran DAXPY runs at 175 Mflop/s and optimized Fortran DGEMM runs at 230 Mflop/s per node. The performance of the NPB (Class B) on the IBM SP2 is compared with the CRAY C90, Intel Paragon, TMC CM-5E, and the CRAY T3D.

Saini, Subhash↗
NASA NTRS2002

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.↗
NASA NTRS2002

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.↗
NASA NTRS2003

On the Hilbert-Huang Transform Data Processing System Development

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 Fourier view of nonlinear mechanics that had existed for a long time, and the associated FFT (fairly recent development), carry strong a-priori assumptions about the source data, such as linearity and of being stationary. Natural phenomena measurements are essentially nonlinear and nonstationary. A very recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT) proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using the Empirical Mode Decomposition (EMD) followed by the Hilbert Transform of the empirical decomposition data (HT), the HHT allows spectrum analysis of nonlinear and nonstationary data by using an engineering a-posteriori data processing, based on the EMD algorithm. This results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF) that can be further analyzed for spectrum interpretation by the classical Hilbert Transform. This paper describes phase one of the development of a new engineering tool, the HHT Data Processing System (HHTDPS). The HHTDPS allows applying the "T to a data vector in a fashion similar to the heritage FFT. It is a generic, low cost, high performance personal computer (PC) based system that implements the HHT computational algorithms in a user friendly, file driven environment. This paper also presents a quantitative analysis for a complex waveform data sample, a summary of technology commercialization efforts and the lessons learned from this new technology development.

Kizhner, Semion↗
NASA NTRS1995

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

Implementation of a Message Passing Interface into a Cloud-Resolving Model for Massively Parallel Computing

The capability for massively parallel programming (MPP) using a message passing interface (MPI) has been implemented into a three-dimensional version of the Goddard Cumulus Ensemble (GCE) model. The design for the MPP with MPI uses the concept of maintaining similar code structure between the whole domain as well as the portions after decomposition. Hence the model follows the same integration for single and multiple tasks (CPUs). Also, it provides for minimal changes to the original code, so it is easily modified and/or managed by the model developers and users who have little knowledge of MPP. The entire model domain could be sliced into one- or two-dimensional decomposition with a halo regime, which is overlaid on partial domains. The halo regime requires that no data be fetched across tasks during the computational stage, but it must be updated before the next computational stage through data exchange via MPI. For reproducible purposes, transposing data among tasks is required for spectral transform (Fast Fourier Transform, FFT), which is used in the anelastic version of the model for solving the pressure equation. The performance of the MPI-implemented codes (i.e., the compressible and anelastic versions) was tested on three different computing platforms. The major results are: 1) both versions have speedups of about 99% up to 256 tasks but not for 512 tasks; 2) the anelastic version has better speedup and efficiency because it requires more computations than that of the compressible version; 3) equal or approximately-equal numbers of slices between the x- and y- directions provide the fastest integration due to fewer data exchanges; and 4) one-dimensional slices in the x-direction result in the slowest integration due to the need for more memory relocation for computation.

Juang, Hann-Ming Henry↗
NASA NTRS1997

Light microscopic image analysis system to quantify immunoreactive terminal area apposed to nerve cells

The present report describes a desktop computer-based method for the quantitative assessment of the area occupied by immunoreactive terminals in close apposition to nerve cells in relation to the perimeter of the cell soma. This method is based on Fast Fourier Transform (FFT) routines incorporated in NIH-Image public domain software. Pyramidal cells of layer V of the somatosensory cortex outlined by GABA immunolabeled terminals were chosen for our analysis. A Leitz Diaplan light microscope was employed for the visualization of the sections. A Sierra Scientific Model 4030 CCD camera was used to capture the images into a Macintosh Centris 650 computer. After preprocessing, filtering was performed on the power spectrum in the frequency domain produced by the FFT operation. An inverse FFT with filter procedure was employed to restore the images to the spatial domain. Pasting of the original image to the transformed one using a Boolean logic operation called 'AND'ing produced an image with the terminals enhanced. This procedure allowed the creation of a binary image using a well-defined threshold of 128. Thus, the terminal area appears in black against a white background. This methodology provides an objective means of measurement of area by counting the total number of pixels occupied by immunoreactive terminals in light microscopic sections in which the difficulties of labeling intensity, size, shape and numerical density of terminals are avoided.

Non-NASA Center↗
NASA NTRS2004

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗