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 217 records · Page 12

Wideband Spectroscopy: The Design and Implementation of a 3 GHz Bandwidth, 8192 Channel, Polyphase Digital Spectrometer

A family of state-of-the-art digital Fourier transform spectrometers has been developed, with a combination of high bandwidth and fine resolution unavailable elsewhere. Analog signals consisting of radiation emitted by constituents in planetary atmospheres or galactic sources are downconverted and subsequently digitized by a pair of interleaved Analog-to-Digital Converters, (ADC). This 6 Gsps (giga-sample per second) digital representation of the analog signal is then processed through an FPGA-based streaming Fast Fourier Transform (FFT), the key development described below. Digital spectrometers have many advantages over previously used analog spectrometers, especially in terms of accuracy and resolution, both of which are particularly important for the type of scientific questions to be addressed with next-generation radiometers. the implementation, results and underlying math for this spectrometer, as well as, potential for future extension to even higher bandwidth, resolution and channel orthogonality, needed to support proposed future advanced atmospheric science and radioastronomy, are discussed.

field programmable gate array↗

Wideband Spectroscopy: The Design and Implementation of a 3 GHz, 2048 Channel Digital Spectrometer

A state-of-the-art digital Fourier Transform spectrometer has been developed, with a combination of high bandwidth and fine resolution unavailable elsewhere. Analog signals consisting of radiation emitted by constituents in planetary atmospheres or galactic sources are downconverted and subsequently digitized by a pair of interleaved Analog-to-Digital Converters (ADC). This 6 Gsps (giga sample per second) digital representation of the analog signal is then processed through an FPGA-based streaming Fast Fourier Transform (FFT), the key development described below. Digital spectrometers have many advantages over previously used analog spectrometers, especially in terms of accuracy and resolution, both of which are particularly important for the type of scientific questions to be addressed with next-generation radiometers. The implementation, results and underlying math for this spectrometer, as well as potential for future extension to even higher bandwidth, resolution and channel orthogonality, needed to support proposed future advanced atmospheric science and radioastronomy, are discussed.

field programmable gate array↗

Fuel Distribution Estimate via Spin Period to Precession Period Ratio for the Advanced Composition Explorer

The spin period to precession period ratio of a non-axisymmetric spin-stabilized spacecraft, the Advanced Composition Explorer (ACE), was used to estimate the remaining mass and distribution of fuel within its propulsion system. This analysis was undertaken once telemetry suggested that two of the four fuel tanks had no propellant remaining, contrary to pre-launch expectations of the propulsion system performance. Numerical integration of possible fuel distributions was used to calculate moments of inertia for the spinning spacecraft. A Fast Fourier Transform (FFT) of output from a dynamics simulation was employed to relate calculated moments of inertia to spin and precession periods. The resulting modeled ratios were compared to the actual spin period to precession period ratio derived from the effect of post-maneuver nutation angle on sun sensor measurements. A Monte Carlo search was performed to tune free parameters using the observed spin period to precession period ratio over the life of the mission. This novel analysis of spin and precession periods indicates that at the time of launch, propellant was distributed unevenly between the two pairs of fuel tanks, with one pair having approximately 20% more propellant than the other pair. Furthermore, it indicates the pair of the tanks with less fuel expelled all of its propellant by 2014 and that approximately 46 kg of propellant remains in the other two tanks, an amount that closely matches the operational fuel accounting estimate. Keywords: Fuel Distribution, Moments of Inertia, Precession, Spin, Nutation

Moment of Inertia↗

Performance of the NASA Beacon Receiver for the Alphasat Aldo Paraboni TDP5 Propagation Experiment

NASA Glenn Research Center (GRC) and the Politecnico di Milano (POLIMI) have initiated a joint propagation campaign within the framework of the Alphasat propagation experiment to characterize rain attenuation, scintillation, and gaseous absorption effects of the atmosphere in the 40 gigahertz band. NASA GRC has developed and installed a K/Q-band (20/40 gigahertz) beacon receiver at the POLIMI campus in Milan, Italy, which receives the 20/40 gigahertz signals broadcast from the Alphasat Aldo Paraboni Technology Demonstration Payload (TDP) no. 5 beacon payload. The primary goal of these measurements is to develop a physical model to improve predictions of communications systems performance within the Q-band. Herein, we describe the design and preliminary performance of the NASA propagation terminal, which has been installed and operating in Milan since June 2014. The receiver is based upon a validated Fast Fourier Transform (FFT) I/Q digital design approach utilized in other operational NASA propagation terminals, but has been modified to employ power measurement via a frequency estimation technique and to coherently track and measure the amplitude of the 20/40 gigahertz beacon signals. The system consists of a 1.2-meter K-band and a 0.6-meter Q-band Cassegrain reflector employing synchronous open-loop tracking to track the inclined orbit of the Alphasat satellite. An 8 hertz sampling rate is implemented to characterize scintillation effects, with a 1-hertz measurement bandwidth dynamic range of 45 decibels. A weather station with an optical disdrometer is also installed to characterize rain drop size distribution for correlation with physical based models.

propagation↗

High-Speed On-Board Data Processing for Science Instruments: HOPS

The project called High-Speed On-Board Data Processing for Science Instruments (HOPS) has been funded by NASA Earth Science Technology Office (ESTO) Advanced Information Systems Technology (AIST) program during April, 2012 â€" April, 2015. HOPS is an enabler for science missions with extremely high data processing rates. In this three-year effort of HOPS, Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) and 3-D Winds were of interest in particular. As for ASCENDS, HOPS replaces time domain data processing with frequency domain processing while making the real-time on-board data processing possible. As for 3-D Winds, HOPS offers real-time high-resolution wind profiling with 4,096-point fast Fourier transform (FFT). HOPS is adaptable with quick turn-around time. Since HOPS offers reusable user-friendly computational elements, its FPGA IP Core can be modified for a shorter development period if the algorithm changes. The FPGA and memory bandwidth of HOPS is 20 GB/sec while the typical maximum processor-to-SDRAM bandwidth of the commercial radiation tolerant high-end processors is about 130-150 MB/sec. The inter-board communication bandwidth of HOPS is 4 GB/sec while the effective processor-to-cPCI bandwidth of commercial radiation tolerant high-end boards is about 50-75 MB/sec. Also, HOPS offers VHDL cores for the easy and efficient implementation of ASCENDS and 3-D Winds, and other similar algorithms. A general overview of the 3-year development of HOPS is the goal of this presentation.

Beyon, Jeffrey↗

Validation of a Manually Oscillating Chair for In-The-Field Assessment of Dynamic Visual Acuity on Crewmembers Within Hours of Returning from Long-Duration Spaceflight

Due to the deconditioned state of crewmembers in the initial hours after landing, it is safer and more practical to perform a vision test while seated in a chair versus walking on a treadmill. The purpose of this study was to validate the ability of a manually operated oscillating chair to produce the oscillatory frequency and displacement equivalent of walking on a treadmill at a 4 mph pace. A fast Fourier transform (FFT)was performed on the vertical trunk acceleration to compare the peak and spread of the distribution of oscillation frequencies for each oscillating condition. Peak oscillation frequencies achieved with the manual chair were lower and more variable than those of treadmill walking and the automatic chair. This can mostly be attributed to operator fatigue. However, DVA scores across conditions were not significantly different, indicating that the manual chair can provide adequate vertical oscillation frequency and displacement with the added advantage of being portable enough for testing outside a laboratory. Furthermore the automatic chair very closely matches the oscillation frequency of treadmill walking, making it an ideal method for testing DVA in a laboratory setting.

Kreutzberg, G. A.↗

Structure in the 3D Galaxy Distribution: III. Fourier Transforming the Universe: Phase and Power Spectra

We demonstrate the effectiveness of a relatively straightforward analysis of the complex 3D Fourier transform of galaxy coordinates derived from redshift surveys. Numerical demonstrations of this approach are carried out on a volume-limited sample of the Sloan Digital Sky Survey redshift survey. The direct unbinned transform yields a complex 3D data cube quite similar to that from the Fast Fourier Transform (FFT) of finely binned galaxy positions. In both cases deconvolution of the sampling window function yields estimates of the true transform. Simple power spectrum estimates from these transforms are roughly consistent with those using more elaborate methods. The complex Fourier transform characterizes spatial distributional properties beyond the power spectrum in a manner different from (and we argue is more easily interpreted than) the conventional multi-point hierarchy. We identify some threads of modern large scale inference methodology that will presumably yield detections in new wider and deeper surveys.

Galactic clusters↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper describes the Dynamic Mode Decomposition (DMD) of the pressures on the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the Ascent Transient Aerodynamics Tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD) of the data matrix. For the matrix of the uPSP measurements of the SLS ATAT, the number of rows is equal to the number of nodes in the grid of the scale model, and the number of columns is equal to the number of frames in the videos taken with 4 Phantom high-speed cameras. In this paper, it is verified that, for the time series with zero mean value, the DMD is equivalent to the decomposition with the Discrete Fourier Transform (DFT). Considering the uPSP is mainly used in the assessment of the unsteady, aerodynamic phenomena, the DMD of the uPSP measurements can be implemented in two steps: (1) subtract the mean value from the uPSP measurement on each of the grid nodes; (2) apply the Fast Fourier Transform (FFT) on the resulting zero-mean time series. The DMD of the uPSP measurements with FFT has two advantages: (1) the computational complexity of FFT is O(N*logN), where N is the length of the time series; (2) compared to the SVD-based DMD algorithm, the DMD with FFT can be easily implemented in parallel processing. A sample matrix of uPSP measurements, at the size of 341 grid nodes and 128 frames, is generated. Figures 1 and 2 show the eigenvalues and the ratios of the eigenvectors, respectively, of the sample matrix, without and with the mean value removed on each of the grid nodes, computed with the SVD-based DMD and the FFT. The figures demonstrate the equivalence of the SVD-based DMD and the decomposition with DFT/FFT for the time series with zero mean value. The results of DMD of the uPSP measurements of the SLS ATAT in September 2019 are presented in the paper. The DMD modes at different frequencies are shown, the aerodynamic phenomena (e.g. shockwave and vortex shedding) are demonstrated and the correlation of the DMD modes with the test configuration parameter (e.g., the Mach Number) is discussed. Figure 3 shows a software tool to visualize the DMD modes. The code to implement the algorithm described in this paper was written in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, and executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper discusses the Dynamic Mode Decomposition (DMD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected with four Phantom high-speed cameras at a constant sample frequency in the Ascent Transient Aerodynamics Test (ATAT) of the Space Launch System (SLS) Block 1 cargo vehicle with the Unitary Plan Wind Tunnel (UPWT) 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD). For the data with zero mean, the DMD is equivalent to the Discrete Fourier Transform (DFT). Since the uPSP is mainly used to determine the unsteady property of the aerodynamic flow, the DMD of the uPSP measurements is implemented in two steps: (1) subtract the mean value from the uPSP measurement; (2) apply the Fast Fourier Transform (FFT) on the resulting data with zero mean. The DMD of the uPSP measurements with FFT has two advantages: (1) the FFT algorithm is well known for its computational efficiency, therefore, compared to the SVD-based DMD algorithm, the DMD with FFT reduces the computation time; (2) the DMD with FFT can be easily implemented in parallel processing. The DMD outputs were generated with the execution in parallel of a code in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, on the NASA Pleiades supercomputer. In this paper, the results of DMD of the uPSP measurements in the tests of Mach sweep runs of the SLS ATAT are presented, and the effectiveness of the DMD of the uPSP measurements in the diagnosis of the unsteady, aerodynamic phenomena is demonstrated. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding for this research was provided by the NASA Aeroscience Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

Regional Variations in the Diurnal Cycle of Tropical Precipitation as Represented by IMERG, ERA5, and Spaceborne Ku Radar

The diurnal cycle of precipitation is highly regional and is typically a product of multiple competing effects that can be highly localized. The diurnal cycle in high precipitation regions such as the Amazon and the Maritime Continent are of particular interest, especially due to the complex coastal effects which take place over the Maritime Continent. The high spatial and temporal resolution provided by the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (GPM) mission (IMERG) dataset, is used in this study to examine the fine-scale features of the diurnal cycle in these regions. Using an 18-year (2000 – 2018) record of IMERG precipitation observations, diurnal and semidiurnal phase and amplitude are calculated using a fast Fourier transform (FFT) method on precipitation averaged for each half-hour of the day at 0.1°x0.1° spatial resolution. We first introduce an objective method of identifying locations where the diurnal signals are robust and strong. Clear patterns of precipitation phase propagation with distance from shore are shown over both regions, with the diurnal phase and amplitude exhibiting a strong dependence on the distance from the coastline. Semidiurnal cycles are generally weaker than the diurnal cycle except in some isolated locations. Similar analysis is also conducted on the ERA5 reanalysis data in order to evaluate the model’s representation of the precipitation diurnal cycle. The model captures the broad scale patterns of diurnal variability but does not capture all the fine scale patterns nor the exact timing that is observed by IMERG. Comparisons are also made to a long record Ku radar dataset created by combining Tropical Rainfall Measuring Mission (TRMM) and GPM observations, thus providing an additional point of comparison for the timing of the ERA5 precipitation peak, since the timing precipitation can be different, even in between observational datasets.

L J Hayden↗

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗

Nitric Oxide Planar Laser Induced Fluorescence Rotational Thermometry Characterization of A Hypersonic Boundary Layer

Nitric oxide (NO) planar laser-induced fluorescence (PLIF) was performed to determine rotational temperature profiles and fluctuations within a hypersonic boundary layer above the surface of a 2.75° half-angle wedge. The experiments were performed in the Texas A&M University Actively Controlled Expansion (ACE) hypersonic blow-down wind tunnel at Mach = 5.7 and 𝑅𝑒 = 6 × 10 6 /m. The NO was introduced to the flow in the settling chamber of the ACE tunnel and probed using two laser sheets near 226 nm. The resulting NO PLIF fluorescence signal was acquired using in-house software, which simultaneously recorded tunnel conditions. After a Fast Fourier Transform (FFT) blurring and statistical treatment was performed, the preshock temperature was evaluated to be 58 ± 2 K (3.5%), while the turbulent boundary layer temperature near the wall was found to be 350 K with fluctuations on the order of ±25 K (7%). The relative temperature fluctuations were determined to be 3 − 5% in the freestream and peaked at 33% in the turbulent boundary layer. The advantages and disadvantages of seeding NO in the tunnel settling chamber for thermometric PLIF measurements are discussed.

Hypersonics↗

NO PLIF Rotational Thermometry Characterization of a Hypersonic Boundary Layer

Nitric oxide (NO) planar laser-induced fluorescence (PLIF) was performed to determine rotational temperature profiles and fluctuations within a hypersonic boundary layer above the surface of a 2.75° half-angle wedge. The experiments were performed in the Texas A&M University Actively Controlled Expansion (ACE) hypersonic blow-down wind tunnel at Mach = 5.7 and 𝑅𝑒 = 6 × 10 6 /m. The NO was introduced to the flow in the settling chamber of the ACE tunnel and probed using two laser sheets near 226 nm. The resulting NO PLIF fluorescence signal was acquired using in-house software, which simultaneously recorded tunnel conditions. After a Fast Fourier Transform (FFT) blurring and statistical treatment was performed, the preshock temperature was evaluated to be 58 ± 2 K (3.5%), while the turbulent boundary layer temperature near the wall was found to be 350 K with fluctuations on the order of ±25 K (7%). The relative temperature fluctuations were determined to be 3 − 5% in the freestream and peaked at 33% in the turbulent boundary layer. The advantages and disadvantages of seeding NO in the tunnel settling chamber for thermometric PLIF measurements are discussed.

Hypersonics↗

Studying Performance Portability of LAMMPS across Diverse GPU-based Platforms

The molecular dynamics simulation software, LAMMPS, utilizes the Kokkos acceleration library to port computation to a diverse set of architectures including those based on GPU accelerators. In addition to Kokkos, LAMMPS contains a vast code base that leverages the CUDA application programming interface using library functions such as cuFFT, CUDA’s fast-fourier transform (FFT) library, and, more recently, also support for AMD’s Heterogeneous Interface for Portability (HIP) that is rapidly growing. While preparing LAMMPS tests for the AMD GPU-based test system precursors to Frontier, we investigated several strategies for accelerating LAMMPS on AMD GPUs, using the AMD Instinct MI100 and MI250X. In this work, we integrated the HIP FFT library, hipFFT, into the particle-particle particle-mesh (PPPM) long-range solver, which allowed the porting of PPPM calculations to the GPUs. Kokkos behavior on the MI100 and MI250X was also investigated through the package kokkos command of LAMMPS, targeting com- munication, memory usage, and particle grid decomposition. The Tersoff, Reax, Lennard-Jones (LJ), EAM, Granular, and PPPM potentials were investigated in this effort, and results from these experiments are provided. The selected potentials were run on Spock (AMD Instinct MI100), Crusher (AMD Instinct MI250X), AFW HPC11 (NVIDIA A100) and Summit (NVIDIA V100), for comparison. Operational roofline models were constructed and analyzed for the Tersoff, Reax, and Lennard-Jones potentials on Crusher and Summit.

Hagerty, Nick↗

Carrier and sampling frequency offset estimation for RF communication with crystal-less nodes

When the ultra-low power mm-scale sensor node does not have a crystal oscillator and phase-lock loop, it inevitably exhibits significant carrier frequency offset (CFO) and sampling frequency offset (SFO) with respect to the reference frequencies in the gateway. This disclosure enables efficient real-time calculation of accurate SFO and CFO at the gateway, thus the ultra-low power mm-scale sensor node can be realized without a costly and bulky clock reference crystal and also power-hungry phase lock loop. In the proposed system, the crystal-less sensor starts transmission with repetitive RF pulses with a constant interval, followed by the data payload using pulse-position modulation (PPM). A proposed algorithm uses a two-dimensional (2D) fast Fourier transform (FFT) based process that identifies the SFO and CFO at the same time to establish successful wireless communication between the gateway and crystal-less sensor nodes.

Kim, Hun-Seok↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

Abstract We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗

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↗