SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1993

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

NASA NTRS1993

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

NASA NTRS1992

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

NASA NTRS1994

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

NASA NTRS1994

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↗

NASA NTRS1994

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

NASA NTRS1994

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

NASA NTRS1994

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

NASA NTRS1995

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

NASA NTRS1995

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

NASA NTRS1994

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↗

NASA NTRS1998

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↗

NASA NTRS1999

HRDI Observations of Inertia-Gravity Waves in the Mesosphere and Lower Thermosphere

Vertical profiles of High-resolution Doppler imager (HRDI) mesospheric winds have small-scale structure (vertical wavelengths between 10 and 20 km) that is virtually always present. Fourier analysis of HRDI zonal and meridional wind profiles have been carried out, and the spectral characteristics are sorted by latitude, month and local time. Power spectral density (PSD) exhibits a universal exp(-km) structure in the 10-20km wavelength regime, with K lying between 2 and 3. The observed PSD for wavelengths between 10 and 20 km is a factor of 3 higher than a null spectrum constructed from HRDI reported error bars multiplied by randomly varying numbers between -1 and +1. Stokes parameters were consolidated by month into Northern and Southern hemisphere middle and high latitudes belts (40-72 degrees), tidal belts (32-16 degrees) and a tropical belt (8S-8N). Vertical waves between 10 and 15 km in wavelength are about 10-15% polarized everywhere. The inferred propagation direction in the middle and high latitude Southern hemisphere is predominantly meridional during solstice, and significantly more zonal during equinoxes. In the tropical belt, the wave orientations are nearly North-South during solstices, with a slightly higher east-west component during equinox. In the tidal belts where the background wind includes a strong meridional tidal wind, the preferred wave orientation has a significant zonal component during equinox. These findings are consistent with the interpretation of wave filtering by the background wind.

Lieberman, Ruth S.↗

NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗

NASA NTRS2002

A Two-Dimensional Linear Bicharacteristic Scheme for Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on one-dimensional electromagnetic wave propagation problems. This memorandum extends the Linear Bicharacteristic Scheme for computational electromagnetics to model lossy dielectric and magnetic materials and perfect electrical conductors in two dimensions. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media and for perfect electrical conductors. Both the Transverse Electric and Transverse Magnetic polarizations are considered. Computational requirements and a Fourier analysis are also discussed. Heterogeneous media are modeled through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for two-dimensional model problems on uniform grids, and the Finite Difference Time Domain (FDTD) algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the two-dimensional explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has less phase velocity error.

Beggs, John H.↗

NASA NTRS2001

An Implicit LU/AF FDTD Method

There has been some recent work to develop two and three-dimensional alternating direction implicit (ADI) FDTD schemes. These ADI schemes are based upon the original ADI concept developed by Peaceman and Rachford and Douglas and Gunn, which is a popular solution method in Computational Fluid Dynamics (CFD). These ADI schemes work well and they require solution of a tridiagonal system of equations. A new approach proposed in this paper applies a LU/AF approximate factorization technique from CFD to Maxwell s equations in flux conservative form for one space dimension. The result is a scheme that will retain its unconditional stability in three space dimensions, but does not require the solution of tridiagonal systems. The theory for this new algorithm is outlined in a one-dimensional context for clarity. An extension to two and threedimensional cases is discussed. Results of Fourier analysis are discussed for both stability and dispersion/damping properties of the algorithm. Results are presented for a one-dimensional model problem, and the explicit FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗

NASA NTRS1993

Effect of upright tilt on ventricular/vascular coupling in chronically instrumented primates

Studies of the hydraulic loading conditions on the heart in humans, especially pulsatile load, have primarily been limited to the supine state. Therefore, we have chosen a nonhuman primate model, the baboon, to assess left ventricular/vascular coupling in both supine and upright positions. Primate subjects were studied by catheterization under sedation and then after surgical implantation of transducers. This allowed the evaluation of postural stress in the chronically instrumented conscious baboon and then after light dissociative doses of ketamine. Basic hemodynamic variables were evaluated for baboons in supine and upright positions. Fourier analysis was applied to aortic pressure and flow to obtain input and characteristic impedance and the ratio of pulsatile (Wp) to total (Wt) left ventricular power (Wp/Wt). The aortic reflected, or backward, pressure was also calculated. Peripheral resistance increased (P = 0.01) and reflected pressure decreased (17.74 +/- 1.50 vs. 15.21 +/- 2 mmHg; P < 0.01) in upright subjects. Characteristic impedance and Wp/Wt were unchanged. Postoperatively, peripheral resistance increased (2,651 +/- 311 vs. 3,667 +/- 276; P < 0.05) and mean power and Wt decreased (P < 0.01) without changes in reflected pressure. All variables were unchanged after light dissociative doses of ketamine. Thus there is no significant change in efficiency of left ventricular/vascular coupling formulated in terms of Wp/Wt or input impedance with postural stress.

NASA Program Space Physiology and Countermeasures↗

NASA NTRS1943

Periodic Heat Transfer at Small Pressure Fluctuations

The effect of cyclic gas pressure variations on the periodic heat transfer at a flat wall is theoretically analyzed and the differential equation describing the process and its solution for relatively. Small pressure fluctuations developed, thus explaining the periodic heat cycle between gas and wall surface. The processes for pure harmonic pressure and temperature oscillations, respectively, in the gas space are described by means of a constant heat transfer coefficient and the equally constant phase angle between the appearance of the maximum values of the pressure and heat flow most conveniently expressed mathematically in the form of a complex heat transfer coefficient. Any cyclic pressure oscillations, can be reduced by Fourier analysis to harmonic oscillations, which result in specific, mutual relationships of heat-transfer coefficients and phase angles for the different harmonics.

Pfriem, H.↗