Search NASA⌕ Search

SEARCH · Search NASA

Results for “fourier expansion”

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 145 records · Page 8

The Hot Bands of Methane Between 5 and 10(micro)m

Experimental line intensities of 1727 transitions arising from nine hot bands in the pentad-dyad system of methane are fitted to first and second order using the effective dipole moment expansion in the polyad scheme. The observed bands are V3-V2, V3-V4, V1-V2, V1-V4, 2V4-V4, V2+V4-V2, V2+V4-V4, 2V2-V4, and the intensities are obtained from long-path spectra recorded with the Fourier transform spectrometer located at Kitt Peak National Observatory. For the second order model, some of the 27 intensity parameters are not linearly independant, and so two methods (extrapolation and effective parameters) are propopsed to model the intensities of the hot bands.

methane hot bands dipole moment expansion Kitt Pea↗

Mode-expansion method for predicting radar signature above rough ocean surfaces at low-grazing angle

The Mode-Expansion Method (MEM) is introduced to calculate the electromagnetic (EM) waves scattered by 2-D rough water surfaces at low-grazing angles. The Electric Field Integral Equation (EFIE) is used in defining the problem and is simplified by using the Impedance Boundary Condition (IBC). The surface currents are expressed as the sum of modes expanded as the Fourier series with incident wave as the dominant mode. It is shown that, by the MEM and for the geometry with transmitting and receiving waves at low-grazing angles, very few modes are needed in solving the forward scattering field with reasonable accuracy.

electromagnetic (EM)↗

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗

Design of a linear projector for use with the normal modes of the GLAS 4th order GCM

The design of a linear projector for use with the normal modes of a model of atmospheric circulation is discussed. A central element in any normal mode initialization scheme is the process by which a set of data fields - winds, temperatures or geopotentials, and surface pressures - are expressed ("projected') in terms of the coefficients of a model's normal modes. This process is completely analogous to the Fourier decomposition of a single field (indeed a FFT applied in the zonal direction is a part of the process). Complete separability in all three spatial dimensions is assumed. The basis functions for the modal expansion are given. An important feature of the normal modes is their coupling of the structures of different fields, thus a coefficient in a normal mode expansion would contain both mass and momentum information.

Bloom, S. C.↗

Analysis of the vibratory excitation of gear systems. II - Tooth error representations, approximations, and application

The second part of a theory for predicting the vibratory excitation of gear systems from fundamental descriptions of gear tooth elastic properties and deviations of tooth faces from perfect involute surfaces is presented. The first part of the theory provides expressions for the Fourier-series coefficients of the vibratory excitation, and this paper gives expressions for these Fourier-series coefficients in terms of easily interpreted gear tooth metrics that are readily evaluated from tooth-face measurements. Results are given for rectangular tooth-face contact regions using two-dimensional Legendre polynomial expansions of local tooth-pair stiffnesses and stiffness-weighted deviations of tooth faces from perfect involute surfaces. A rigorous transfer function approach is developed that permits separation of the effects of gear tooth errors and gear design parameters; the theory is applicable to helical and spur gears and is illustrated with measurements of tooth-spacing errors and tooth profiles obtained from a pair of spur gears.

Mark, W. D.↗

Time-dependent treatment of scattering. II - Novel integral equation approach to quantum wave packets

The novel wave-packet propagation scheme presented is based on the time-dependent form of the Lippman-Schwinger integral equation and does not require extensive matrix inversions, thereby facilitating application to systems in which some degrees of freedom express the potential in a basis expansion. The matrix to be inverted is a function of the kinetic energy operator, and is accordingly diagonal in a Bessel function basis set. Transition amplitudes for various orbital angular momentum quantum numbers are obtainable via either Fourier transform of the amplitude density from the time to the energy domain, or the direct analysis of the scattered wave packet.

Sharafeddin, Omar A.↗

MOSAIC: An Effective FFT-based Numerical Method to Assess Aging Properties of Concrete

As the nuclear fleet in the United States ages and subsequent license renewal applications grow, the prediction of concrete durability at extended operation becomes more important. To address this issue, a Fast-Fourier Transform (FFT) method is utilized to simulate aging-related degradation of concrete within the Microstructure Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) software. MOSAIC utilizes compositional phase maps to simulate damage from radiation-induced volumetric expansion (RIVE), applied force, creep, and thermal expansion. This compositional detail allows each mineral in the microstructure to be assigned specific material properties, allowing the simulation to be as accurate and representative as possible. The principal goal of MOSAIC is to simulate the effects of nonlinear aging mechanisms occurring in nuclear concrete on the macroscopic mechanical properties, using only the aggregate microstructure compositional information as a starting point. Here, several realistic example simulations are shown to demonstrate the utility and uniqueness of the MOSAIC software.

36 MATERIALS SCIENCE↗

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

An investigation of the diffraction of an acoustic plane wave by a curved surface of finite impedance

Phenomena associated with long range propagation of sound over irregular topography motivated this work, which was to analyze the diffraction effects which would occur near the tops of hills and ridges. The diffraction of a high frequency plane wave due to its grazing of a two-dimensional curved surface of finite impedance was also studied. Laboratory scale models were constructed and measurements were made of the field on, above, and behind either of two curved surfaces possessing distinctly different impedances; that is, one was soft while the other was hard. The experimental technique consisted of simultaneously measuring the pressure at a reference point and at a field point due to a transient pulse generated by an electric spark. The pressure waveforms were digitized and processed. The ratio of the discrete Fourier transforms of the two waveforms provided an estimate of the insertion loss between them. The results of the measurements were compared with the predictions of a theory which was derived by Pierce using the method of Matched Asymptotic Expansions (MAE). The predictions relied upon the experimental evaluation of the impedance of each surface at grazing angles of incidence. This evaluation was achieved by a fairly standard technique involving empirical models of various generic types of surfaces. An example was shown of the important role that the structural intricacies of a surface play in the determination of an appropriate model. The comparison between the measurements and predictions indicated that the theory gives an excellent description of the field anywhere near a curved surface. Further, with a simple modification, the theory was also shown to give nearly as good of a description of the field surrounding a curved surface even at distances far behind the surface yet near the line of sight.

Kearns, James A.↗

Ideal shrinking and expansion of discrete sequences

Ideal methods are described for shrinking or expanding a discrete sequence, image, or image sequence. The methods are ideal in the sense that they preserve the frequency spectrum of the input up to the Nyquist limit of the input or output, whichever is smaller. Fast implementations that make use of the discrete Fourier transform or the discrete Hartley transform are described. The techniques lead to a new multiresolution image pyramid.

Watson, Andrew B.↗

An Overdetermined System for Improved Autocorrelation Based Spectral Moment Estimator Performance

Autocorrelation based spectral moment estimators are typically derived using the Fourier transform relationship between the power spectrum and the autocorrelation function along with using either an assumed form of the autocorrelation function, e.g., Gaussian, or a generic complex form and applying properties of the characteristic function. Passarelli has used a series expansion of the general complex autocorrelation function and has expressed the coefficients in terms of central moments of the power spectrum. A truncation of this series will produce a closed system of equations which can be solved for the central moments of interest. The autocorrelation function at various lags is estimated from samples of the random process under observation. These estimates themselves are random variables and exhibit a bias and variance that is a function of the number of samples used in the estimates and the operational signal-to-noise ratio. This contributes to a degradation in performance of the moment estimators. This dissertation investigates the use autocorrelation function estimates at higher order lags to reduce the bias and standard deviation in spectral moment estimates. In particular, Passarelli's series expansion is cast in terms of an overdetermined system to form a framework under which the application of additional autocorrelation function estimates at higher order lags can be defined and assessed. The solution of the overdetermined system is the least squares solution. Furthermore, an overdetermined system can be solved for any moment or moments of interest and is not tied to a particular form of the power spectrum or corresponding autocorrelation function. As an application of this approach, autocorrelation based variance estimators are defined by a truncation of Passarelli's series expansion and applied to simulated Doppler weather radar returns which are characterized by a Gaussian shaped power spectrum. The performance of the variance estimators determined from a closed system is shown to improve through the application of additional autocorrelation lags in an overdetermined system. This improvement is greater in the narrowband spectrum region where the information is spread over more lags of the autocorrelation function. The number of lags needed in the overdetermined system is a function of the spectral width, the number of terms in the series expansion, the number of samples used in estimating the autocorrelation function, and the signal-to-noise ratio. The overdetermined system provides a robustness to the chosen variance estimator by expanding the region of spectral widths and signal-to-noise ratios over which the estimator can perform as compared to the closed system.

Keel, Byron M.↗

Two-dimensional representations of axisymmetric fields for computer calculations

An accurate representation of axisymmetric fields has been devised by extending the method of ideal current loops to off-axis fields. It is assumed that the data to be simulated are available through measurements or through a solution of a boundary value problem. The method provides an algebraic expression for the fields throughout a two-dimensional region of interest and eliminates the need for the axial expansion formula in approximating off-axis fields. The use of Gaussian and other functions as alternatives to the coil function is proposed. Examples of the technique in simulating a periodic permanent magnet (PPM) focusing field are presented and compared with a Fourier analysis of the problem.

Stankiewicz, N.↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

A singularity free analytical solution of artificial satellite motion with drag

The connection between the existing Delaunay-Similar and Poincare-Similar satellite theories in the true anomaly version is outlined for the J(2) perturbation and the new drag approach. An overall description of the concept of the approach is given while the necessary expansions and the procedure to arrive at the computer program for the canonical forces is delineated. The procedure for the analytical integration of these developed equations is described. In addition, some numerical results are given. The computer program for the algebraic multiplication of the Fourier series which creates the FORTRAN coding in an automatic manner is described and documented.

Scheifele, G.↗

Microstructural Simulation of Ion-Irradiated Natural Rocks and Minerals

This report details the contribution of Oak Ridge National Laboratory to the Nuclear Energy University Program (NEUP) project, "Rapid Characterization of Concrete Mineralogy using Multi-Scale Tools" led by the University of Illinois at Urbana-Champaign. To support the subsequent license renewal (SLR) of the US nuclear power plants (NPPs) fleet, the performance of concrete-forming aggregates against neutron irradiation needs to be assessed. Ion irradiation is proposed as a rapid, cost-effective surrogate method to full neutron irradiation testing. To complement the characterization of ion-irradiated rock specimens, ORNL has run finite-elements and fast-Fourier transform (FFT-based) simulations using the codes MARS and Microstructure-Oriented Scientific Analysis of Irradiated Concrete (MOSAIC). The main conclusion of this analysis is that the apparent post–ion-irradiation step-height underestimates the accumulated free radiation-induced volumetric expansion (RIVE) in the ion-implanted depth by about 15% at full amorphization and about 25% toward the beginning of the ion irradiation experiment. The discrepancy is explained by the fact that the step height is proportional to the sum of the RIVE (i.e., one third of the RIVE) and the irradiation-assisted plastic strains in the vertical direction (lower than two thirds of the RIVE). Because of the large lateral compressive stresses caused by the restraining effect of the unirradiated substrate, the stress field in the mineral grains and at the grain boundary (GB) in the ion-implanted layer is much different than the comparable stress field occurring during neutron irradiation. Hence, the mismatch RIVE causing cracks in the rock-forming minerals leads to different cracking patterns. Ion irradiation continues to be used as a rapid technique to assess the RIVE potential of rocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Dynamic Phasor Modeling of Various Multipulse Rectifiers and a VSI Fed by 18-Pulse Asymmetrical Autotransformer Rectifier Unit for Fast Transient Analysis

To study the fast modelling of dynamic characteristics caused by power electronic switching, the dynamic phasor (DP) theory based on the time-varying Fourier decomposition and frequency shifting is firstly applied to the modelling of various multipulse rectifiers. The DP model for symmetrical 12-pulse phase-shifting reactor rectifier unit (PSR-RU) is derived with model order reduction, relations between ac and dc terminals, and Taylor series expansion. The DP model of asymmetrical 18-pulse autotransformer rectifier unit (AT-RU) is proposed based on the switching functions expressed in the DP domain. Meanwhile, the DP model of voltage source inverter (VSI) fed by asymmetrical 18-pulse AT-RU is built with the harmonic state-space (HSS) equations. Under both balanced and unbalanced conditions, the good calculation accuracy and rapid simulation speed of the developed DP models are validated by the detailed time-domain (TD) simulation.

42 ENGINEERING↗

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗