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 379 records · Page 21

Application of an improved Nelson-Nguyen analysis to eccentric, arbitrary profile liquid annular seals

An improved dynamic analysis for liquid annular seals with arbitrary profile based on a method, first proposed by Nelson and Nguyen, is presented. An improved first order solution that incorporates a continuous interpolation of perturbed quantities in the circumferential direction, is presented. The original method uses an approximation scheme for circumferential gradients, based on Fast Fourier Transforms (FFT). A simpler scheme based on cubic splines is found to be computationally more efficient with better convergence at higher eccentricities. A new approach of computing dynamic coefficients based on external specified load is introduced. This improved analysis is extended to account for arbitrarily varying seal profile in both axial and circumferential directions. An example case of an elliptical seal with varying degrees of axial curvature is analyzed. A case study based on actual operating clearances of an interstage seal of the Space Shuttle Main Engine High Pressure Oxygen Turbopump is presented.

Padavala, Satyasrinivas↗

A fast numerical solution of scattering by a cylinder: Spectral method for the boundary integral equations

It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.

Hu, Fang Q.↗

In-flight estimation of gyro noise on the Upper Atmosphere Research Satellite (UARS) and Extreme Ultraviolet Explorer (EUVE) missions

This paper characterizes the low-frequency noise response of the Teledyne dry rotor inertial reference unit (DRIRU) gyroscopes on the Upper Atmosphere Research Satellite (UARS) and the Extreme Ultraviolet Explorer (EUVE). The accuracy of spacecraft attitude estimation algorithms that use gyro data for propagating the spacecraft attitude is sensitive to gyro noise. EUVE gyro data were processed to validate a single-axis gyro noise model, which is used onboard various spacecraft. The paper addresses the potential impact of temperature effects on the gyro noise model and the overall impact on attitude determination accuracy. The power spectral density (PSD) of the gyro noise is estimated from UARS in-flight data by Fast Fourier Transform (FFT). The role of actuator dynamics on the PSD function is also discussed.

Lee, M.↗

Parallel computational fluid dynamics - Implementations and results

The present volume on parallel CFD discusses implementations on parallel machines, numerical algorithms for parallel CFD, and performance evaluation and computer science issues. Attention is given to a parallel algorithm for compressible flows through rotor-stator combinations, a massively parallel Euler solver for unstructured grids, a fast scheme to analyze 3D disk airflow on a parallel computer, and a block implicit multigrid solution of the Euler equations. Topics addressed include a 3D ADI algorithm on distributed memory multiprocessors, clustered element-by-element computations for fluid flow, hypercube FFT and the Fourier pseudospectral method, and an investigation of parallel iterative algorithms for CFD. Also discussed are fluid dynamics using interface methods on parallel processors, sorting for particle flow simulation on the connection machine, a large grain mapping method, and efforts toward a Teraflops capability for CFD.

Simon, Horst D.↗

Analysis of the Effects of Surface Pitting and Wear on the Vibrations of a Gear Transmission System

A comprehensive procedure to simulate and analyze the vibrations in a gear transmission system with surface pitting, 'wear' and partial tooth fracture of the gear teeth is presented. An analytical model was developed where the effects of surface pitting and wear of the gear tooth were simulated by phase and magnitude changes in the gear mesh stiffness. Changes in the gear mesh stiffness were incorporated into each gear-shaft model during the global dynamic simulation of the system. The overall dynamics of the system were evaluated by solving for the transient dynamics of each shaft system simultaneously with the vibration of the gearbox structure. In order to reduce the number of degrees-of-freedom in the system, a modal synthesis procedure was used in the global transient dynamic analysis of the overall transmission system. An FFT procedure was used to transform the averaged time signal into the frequency domain for signature analysis. In addition, the Wigner-Ville distribution was also introduced to examine the gear vibration in the joint time frequency domain for vibration pattern recognition. Experimental results obtained from a gear fatigue test rig at NASA Lewis Research Center were used to evaluate the analytical model.

Choy, F. K.↗

An algorithm for extraction of periodic signals from sparse, irregularly sampled data

Temporal gaps in discrete sampling sequences produce spurious Fourier components at the intermodulation frequencies of an oscillatory signal and the temporal gaps, thus significantly complicating spectral analysis of such sparsely sampled data. A new fast Fourier transform (FFT)-based algorithm has been developed, suitable for spectral analysis of sparsely sampled data with a relatively small number of oscillatory components buried in background noise. The algorithm's principal idea has its origin in the so-called 'clean' algorithm used to sharpen images of scenes corrupted by atmospheric and sensor aperture effects. It identifies as the signal's 'true' frequency that oscillatory component which, when passed through the same sampling sequence as the original data, produces a Fourier image that is the best match to the original Fourier space. The algorithm has generally met with succession trials with simulated data with a low signal-to-noise ratio, including those of a type similar to hourly residuals for Earth orientation parameters extracted from VLBI data. For eight oscillatory components in the diurnal and semidiurnal bands, all components with an amplitude-noise ratio greater than 0.2 were successfully extracted for all sequences and duty cycles (greater than 0.1) tested; the amplitude-noise ratios of the extracted signals were as low as 0.05 for high duty cycles and long sampling sequences. When, in addition to these high frequencies, strong low-frequency components are present in the data, the low-frequency components are generally eliminated first, by employing a version of the algorithm that searches for non-integer multiples of the discrete FET minimum frequency.

Wilcox, J. Z.↗

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

Extracting Periodic Signals From Irregularly Sampled Data

Successive approximations formed in Fourier space. Algorithm extracts periodic signals from sparse, irregularly sampled sets of measurement data. Pertains to data processed via fast Fourier transforms (FFTs). Data represents signal components with initially unknown frequencies spanning large spectral range and includes frequencies not integer multiples of minimum FFT frequency.

Wilcox, Jaroslava Z.↗

An experimental study of a turbulent boundary layer in the trailing edge region of a circulation-control airfoil

This report describes the development of a new near-wall LDV technique and its use to measure the full turbulent stress tensor in a highly three-dimensional separated boundary layer flow. Laser Doppler velocimetry, with small diameter laser beams allowing small beam-crossing angles (i.e., larger fringe spacings), was used, resulting in measurements close to solid surfaces. Dennis Johnson of NASA/Ames combined the advantages offered by FFT signal processors with a beam-steering probe inserted in the flow to devise a new self-aligning, near-wall LDV approach capable of making direct measurements of the transverse velocity component. This new technique was applied to a 2D, zero-pressure-gradient, flat-plate boundary layer in the High Reynolds Number Pilot Channel at Ames. Accurate and repeatable measurements of U and uu to within 30 mm of the wall were obtained. Next, the near-wall LDV technique was employed in a 3D separated boundary layer flow, consisting of a 2-dimensional wall boundary layer approaching a right-angle wedge. The resulting lateral and adverse streamwise pressure gradients generated significant spanwise velocities and turbulent stresses, and a large reverse-flow region. The near-wall LDV was modified further to facilitate measurements in the 3D Wedge Flow.

Brown, Jeffrey↗

Periodicities of hard x-ray burst during the last solar cycle

By using power spectrum and standard FFT time series analysis, the Hard X-ray burst during solar cycle -22 were studied. This data of Hard X-ray burst spectrometer (HXRBS) on the solar maximum mission from Launch and February 14, 1980, through re-entry on December 2, 1989, by NASA artificial satellite. The results indicate that there are short and intermediate solar periodicities. Also it is found that there is a relation between the short periodicities (few minutes) with similar periodicities in solar radio emissions and in good agreement with the theoretical mode of solar oscillations.

Hady, Ahmed A.↗

Solar Dynamic Power System Fault Diagnosis

The objective of this research is to conduct various fault simulation studies for diagnosing the type and location of faults in the power distribution system. Different types of faults are simulated at different locations within the distribution system and the faulted waveforms are monitored at measurable nodes such as at the output of the DDCU's. These fault signatures are processed using feature extractors such as FFT and wavelet transforms. The extracted features are fed to a clustering based neural network for training and subsequent testing using previously unseen data. Different load models consisting of constant impedance and constant power are used for the loads. Open circuit faults and short circuit faults are studied. It is concluded from present studies that using features extracted from wavelet transforms give better success rates during ANN testing. The trained ANN's are capable of diagnosing fault types and approximate locations in the solar dynamic power distribution system.

Momoh, James A.↗

Coincident 1.3-year Periodicities in the ap Geomagnetic Index and the Solar Wind

Recent observations show an approximately 1.3-year period in the speed of the solar wind detected by the IMP 8 and Voyager 2 spacecraft. A similar period is also seen in the north-south (GSE) component of the magnetic field observed by IMP 8. Since both parameters are commonly used as input to models of geomagnetic activity, the 'ap' index (a measure of geomagnetic disturbance) is examined to look for this periodicity. The Lomb-Scargle periodogram method is used on the ap, plasma, and magnetic field data during the 1973-1994 time range. A dynamic FFT periodogram method is also used to analyze the ap data during this time, as well as to look for periods present between 1932 and 1972. A clear 1.3-year periodicity is present in the post-1986 data when the same period is observed in the plasma and field data. The V(2)B(zsm) and V(2)B(s) proxies for geomagnetic activity also show this periodicity. However, the southward (GSM) component of the magnetic field does not have a 1.3-year period, and neither do solar wind or ap data from 1973-1985. This demonstrates that the ap geomagnetic index can act as a proxy for solar wind periodicities at this time scale. Historic ap data are examined, and show that a similar periodicity in ap exists around 1942. Since auroral data show a 1.4-year periodicity, all these similar periods may result from a common underlying solar mechanism.

Paularena, K. I.↗

A High-Order Direct Solver for Helmholtz Equations with Neumann Boundary Conditions

In this study, a compact finite-difference discretization is first developed for Helmholtz equations on rectangular domains. Special treatments are then introduced for Neumann and Neumann-Dirichlet boundary conditions to achieve accuracy and separability. Finally, a Fast Fourier Transform (FFT) based technique is used to yield a fast direct solver. Analytical and experimental results show this newly proposed solver is comparable to the conventional second-order elliptic solver when accuracy is not a primary concern, and is significantly faster than that of the conventional solver if a highly accurate solution is required. In addition, this newly proposed fourth order Helmholtz solver is parallel in nature. It is readily available for parallel and distributed computers. The compact scheme introduced in this study is likely extendible for sixth-order accurate algorithms and for more general elliptic equations.

Sun, Xian-He↗

Three Dimensional Viscous Flow Field in an Axial Flow Turbine Nozzle Passage

The objective of this investigation is experimental and computational study of three dimensional viscous flow field in the nozzle passage of an axial flow turbine stage. The nozzle passage flow field has been measured using a two sensor hot-wire probe at various axial and radial stations. In addition, two component LDV measurements at one axial station (x/c(sum m) = 0.56) were performed to measure the velocity field. Static pressure measurements and flow visualization, using a fluorescent oil technique, were also performed to obtain the location of transition and the endwall limiting streamlines. A three dimensional boundary layer code, with a simple intermittency transition model, was used to predict the viscous layers along the blade and endwall surfaces. The boundary layers on the blade surface were found to be very thin and mostly laminar, except on the suction surface downstream of 70% axial chord. Strong radial pressure gradient, especially close to the suction surface, induces strong cross flow components in the trailing edge regions of the blade. On the end-walls the boundary layers were much thicker, especially near the suction corner of the casing surface, caused by secondary flow. The secondary flow region near the suction-casing surface corner indicates the presence of the passage vortex detached from the blade surface. The corner vortex is found to be very weak. The presence of a closely spaced rotor downstream (20% of the nozzle vane chord) introduces unsteadiness in the blade passage. The measured instantaneous velocity signal was filtered using FFT square window to remove the periodic unsteadiness introduced by the downstream rotor and fans. The filtering decreased the free stream turbulence level from 2.1% to 0.9% but had no influence on the computed turbulence length scale. The computation of the three dimensional boundary layers is found to be accurate on the nozzle passage blade surfaces, away from the end-walls and the secondary flow region. On the nozzle passage endwall surfaces the presence of strong pressure gradients and secondary flow limit the validity of the boundary layer code.

Ristic, D.↗

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

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

A Search for Short-Period Accreting Pulsars with BATSE

Accreting X-ray pulsar systems can be detected and monitored by the BATSE Pulsar Monitor at full sensitivity only above the Nyquist period of 2.048 s, since it uses the DISCLA data type. Systems with shorter periods than 2 s can be observed, but with a reduced sensitivity, since only the aliased power is measured. While persistent sources in the period range below that have been found by earlier spacecraft, transient systems may still be undiscovered, due to limited time coverage of the galactic plane with all-sky instruments. BATSE data with a time resolution of 31 ms (PSR single-sweep mode) has been collected for four energy channels spanning 20-50 keV , at scheduled times for particular detectors and detector combinations which view the galactic plane. The data spans Nov 1993 to present, with a typical exposure of 15000 s per day to approximately 1/8 of the plane. An accelerated FFT-based search is in progress on these data. Data are segmented into approximately 300 second intervals, barycentered for 3 possible source directions spaced by 30 degrees about the detector normal, with the time series accelerated for each of 15 steps, each spaced by 5 x 10(exp -5). To avoid loss of coherence due to system orbital velocities, FFTs are averaged over each 1/4 day interval separately. Results of this search are presented, including examples of detection of Her X-1, Cen X-3, and other pulsars.

Wilson, R. B.↗

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