Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fourier analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

An optimized finite-difference scheme for wave propagation problems

Two fully-discrete finite-difference schemes for wave propagation problems are presented, a maximum-order scheme and an optimized (or spectral-like) scheme. Both combine a seven-point spatial operator and an explicit six-stage time-march method. The maximum-order operator is fifth-order in space and is sixth-order in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a second-order and a fourth-order method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reduction in global error compared to the maximum-order scheme with no additional computational expense.

Zingg, D. W.↗

Direct calculations of waves in fluid flows using a high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.

Yu, Sheng-Tao↗

Normal incidence X-ray telescope power spectra of X-ray emission from solar active regions. I - Observations. II - Theory

Fourier analysis is applied to very high resolution image of coronal active regions obtained by the Normal Incidence X-Ray Telescope is used to find a broad isotropic power-law spectrum of the spatial distribution of soft X-ray intensities. Magnetic structures of all sizes are present down to the resolution limit of the instrument. Power spectra for the X-ray intensities of a sample of topologically different active regions are found which fall off with increasing wavenumber as 1/k-cubed. A model is presented that relates the basic features of coronal magnetic fluctuations to the subphotospheric hydrodynamic turbulence that generates them. The model is used to find a theoretical power spectrum for the X-ray intensity which falls off with increasing wavenumber as 1/k-cubed. The implications of a turbulent regime in active regions are discussed.

Gomez, Daniel O.↗

Position and morphology of the compact non-thermal radio source at the Galactic Center

We have determined with VLBI the position of the compact nonthermal radio source at the Galactic Center, commonly referred to as SgrA*, in the J2000.0 reference frame of extragalactic radio sources. We have also determined the size of SgrA* at 1.3, 3.6, and 13 cm wavelengths and found that the apparent size of the source increases proportionally to the observing wavelength squared, as expected from source size broadening by interstellar scattering and as reported previously by other authors. We have also established an upper limit of about 8 mJy at 3.6 cm wavelength for any ultracompact component. The actual size of the source is less than 15 AU. Fourier analysis of our very sensitive 3.6 cm observations of this source shows no significant variations of correlated flux density on time scales from 12 to 700 s.

Marcaide, J. M.↗

Dynamic Response of Turbine-blade Temperature to Exhaust-gas Temperature for Gas-turbine Engines

The frequency of blade temperature to exhaust gas temperature is presented for two locations in the blade and at several operating conditions. The frequency response was determined by Fourier analysis of transient data. Two analytical methods are presented, and results are compared with experimental data. Dynamic response of turbine-blade temperature to exhaust-gas temperature exhibited the form of an approximate first-order lag. The results are used to predict blade temperature for a typical controlled and uncontrolled gas-turbine engine.

HEAT TRANSFER THEORY AND EXPERIMENT↗

On the statistical theory of turbulence

A study is made of the spectrum of isotropic turbulence with the aid of the customary method of Fourier analysis. The spectrum of the turbulent motion is derived to the smallest wave lengths, that is, into the laminar region, and correlation functions and pressure fluctuations are calculated. A comparison with experimental results is included. Finally, an attempt is made to derive the numerical value of a constant characteristic of the energy dissipation in isotropic turbulence.

Heisenberg, W↗

Effects of window size and shape on accuracy of subpixel centroid estimation of target images

A new algorithm is presented for increasing the accuracy of subpixel centroid estimation of (nearly) point target images in cases where the signal-to-noise ratio is low and the signal amplitude and shape vary from frame to frame. In the algorithm, the centroid is calculated over a data window that is matched in width to the image distribution. Fourier analysis is used to explain the dependency of the centroid estimate on the size of the data window, and simulation and experimental results are presented which demonstrate the effects of window size for two different noise models. The effects of window shape were also investigated for uniform and Gaussian-shaped windows. The new algorithm was developed to improve the dynamic range of a close-range photogrammetric tracking system that provides feedback for control of a large gap magnetic suspension system (LGMSS).

Welch, Sharon S.↗

Multigrid schemes for viscous hypersonic flows

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm employs upwind spatial discretization with explicit multistage time stepping. Two-level versions of the various multigrid algorithms are applied to the two-dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving two different hypersonic flow problems. Some new multigrid schemes, based on semicoarsening strategies, are shown to be quite effective in relieving the stiffness caused by the high-aspect-ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6).

Swanson, R. C.↗

Detector noise statistics in the non-linear regime

The statistical behavior of an idealized linear detector in the presence of threshold and saturation levels is examined. It is assumed that the noise is governed by the statistical fluctuations in the number of photons emitted by the source during an exposure. Since physical detectors cannot have infinite dynamic range, our model illustrates that all devices have non-linear regimes, particularly at high count rates. The primary effect is a decrease in the statistical variance about the mean signal due to a portion of the expected noise distribution being removed via clipping. Higher order statistical moments are also examined, in particular, skewness and kurtosis. In principle, the expected distortion in the detector noise characteristics can be calibrated using flatfield observations with count rates matched to the observations. For this purpose, some basic statistical methods that utilize Fourier analysis techniques are described.

Shopbell, P. L.↗

The complex lightcurve of 1992 NA

Amor asteroid 1992 NA was monitored during three nights at a large phase angle of -65 deg. The lightcurves obtained did not reveal a repeatable curve with two maxima and two minima. However, some features suggested a periodicity with three maxima and three minima. A satisfactory composite lightcurve of this form was obtained by means of an 'eyeball' fit and by Fourier analysis. Individual and composite lightcurves are presented. The observed colors are consistent with the C class.

Wisniewski, Wieslaw Z.↗

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

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

Weng, Hengyi↗

Large-scale variations in ozone from the first two years of UARS MLS data

Two years of stratospheric measurements of ozone from the Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) are examined in order to characterize large horizontal scale wave variations. The use of Fourier analysis allows the detection of variations from daily through seasonal and interannual timescales. High-latitude winter variations at 10 hPa often have very large amplitudes, but smaller midlatitude variations are more ubiquitous. Some variations have the characteristics of locally generated instabilities. Correlations of wave features with changes in the zonal-mean ozone at 10 hPa suggest the presence of significant horizontal motions during strong wintertime polar warming events. Such correlations are not evident at other levels. Spectral analysis of the large-scale variations show most waves to be slowly propagating. In contrast to some past observations, equatorial regions are shown to lack large amplitude wave events.

Elson, Lee S.↗

Comparison of low-latitude ion and neutral zonal drifts using DE 2 data

We have used data from the ion drift meter and the wind and temperature spectrometer on the DE 2 spacecraft to make statistical comparisons of the zonal ion and neutral drifts at dip latitudes (DLAT) in the +/- 35 deg range over all local times. Fourier analysis indicates that the superrotation and the diurnal components of both flows are strongly peaked at the dip equator, with the superrotation term becoming negative for the absolute value DLAT greater than or equal 20 deg. One interesting feature is the presence of a period (2200-0500 solar local time) in the 300-400 km altitude region near the dip equator where the ion drift is more strongly eastward than the neutral flow. This would seem to indicate the presence of an electric field source of greater strength than the F region dynamo elsewhere along the geomagnetic field line. Model calculations indicate that a possible mechanism for this source lies in the vertical shear in the zonal neutral wind in the 100-200 km altitude region.

Coley, W. R.↗

Can we identify lensed gamma-ray bursts?

A gravitationally lensed gamma-ray burst (GRB) would appear as multiple bursts with identical light curves, separated in time and differing only by the scaling of their amplitudes. In reality, the detected bursts will be noisy, and therefore they may be difficult to identify as lensed images. Furthermore, faint, intrinsically similar, yet distinct light curves may be falsely identified as lensing events. In this paper we develop two statistical tests to distinguish noisy burst light curves. We use Fourier analysis techniques to analyze the signals for both intrinsic variability and variability due to noise. We are able to determine the noise level, and we compare the bursts only at frequency channels that are signal dominated. Utilizing these methods, we are able to make quantitative statements about whether two bursts are distinct. We apply these statistics to scaled versions of two subbursts of GRB 910503 -- subbursts previously investigated by Wambsganss using a different statistical test. We find that our methods are able to distinguish these bursts at slightly smaller amplitudes than those at which Wambsganss's method succeeds. We then apply our techniques to 'candidate' lensing events taken from the Burst and Transient Source Experiment (BATSE) catalog, and we find that nearly all of them, except for the very shortest ones (durations approx. less than 3 sec), are distinguishable. We therefore expect that a majority of bursts will be distinguishable from one another.

Nowak, Michael A.↗

Lightning-channel morphology by return-stroke radiation field waveforms

Simultaneous video and wideband electric field recordings of 32 cloud-to-ground lightning flashes in Florida were analyzed to show the formation of new channels to ground can be detected by examination of the return-stroke radiation fields alone. The return-stroke E and dE/dt waveforms were subjectively classified according to their fine structure. Then the video images were examined field by field to identify each waveform with a visible channel to ground. Fifty-five correlated waveforms and channel images were obtained. Of these, all 34 first-stroke waveforms (multiple jagged E peaks, noisy dE/dt), 8 of which were not radiated by the chronologically first stroke in the flash, came from new channels to ground (not previously seen on video). All 18 subsequent-stroke waveforms (smoothly rounded E and quiet dE/dt after initial peak) were radiated by old channels (illuminated by a previous stroke). Two double-ground waveforms (two distinct first-return-stroke pulses separated by tens of microseconds or less) coincided with video fields showing two new channels. One `anomalous-stroke' waveform (beginning like a first stroke and ending like a subsequent) was produced by a new channel segment to ground branching off an old channel. This waveform classification depends on the presence or absence of high-frequency fine structure. Fourier analysis shows that first-stroke waveforms contain about 18 dB more spectral power in the frequency interval from 500 kHz to at least 7 MHz than subsequent-stroke waveforms for at least 13 microseconds after the main peak.

Willett, J. C.↗

Two-dimensional structure of slow shocks

The two-dimensional (2-D) structure of switch-off slow magnetosonic shocks is investigated using an electromagnetic hybrid (fluid electrons, kinetic ions) code. It is shown that the basic physical processes occurring at 1-D slow shocks are also operative in 2-D. Specifically, the interaction between the upstream ions and those streaming away from the shock results in the excitation of Alfven/ions-cyclotron (AIC) waves. Depending on the plasma parameters, these waves may either stay in the upstream or convect back into the shock resulting in a non-steady shock behavior which prevents the formation of a trailing wave train. Despite this similarity, some slow shocks which are steady in 1-D are found to be non-steady in 2-D. Fourier analysis of the waves downstream of non-steady shocks identifies them as AIC, demonstrating that the waves remain on the same branch as they convect from upstream into the downstream region.

Omidi, N.↗

Direct calculations of waves in fluid flows using high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.

Yu, Sheng-Tao↗

Effects of Helicity on Lagrangian and Eulerian Time Correlations in Turbulence

Taylor series expansions of turbulent time correlation functions are applied to show that helicity influences Eulerian time correlations more strongly than Lagrangian time correlations: to second order in time, the helicity effect on Lagrangian time correlations vanishes, but the helicity effect on Eulerian time correlations is nonzero. Fourier analysis shows that the helicity effect on Eulerian time correlations is confined to the largest inertial range scales. Some implications for sound radiation by swirling flows are discussed.

Rubinstein, Robert↗