Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fast Fourier Transform”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Dark Energy Spectroscopic Instrument: one-dimensional power spectrum from first Ly α forest samples with Fast Fourier Transform

ABSTRACT We present the one-dimensional Ly α forest power spectrum measurement using the first data provided by the Dark Energy Spectroscopic Instrument (DESI). The data sample comprises 26 330 quasar spectra, at redshift z > 2.1, contained in the DESI Early Data Release and the first 2 months of the main survey. We employ a Fast Fourier Transform (FFT) estimator and compare the resulting power spectrum to an alternative likelihood-based method in a companion paper. We investigate methodological and instrumental contaminants associated with the new DESI instrument, applying techniques similar to previous Sloan Digital Sky Survey (SDSS) measurements. We use synthetic data based on lognormal approximation to validate and correct our measurement. We compare our resulting power spectrum with previous SDSS and high-resolution measurements. With relatively small number statistics, we successfully perform the FFT measurement, which is already competitive in terms of the scale range. At the end of the DESI survey, we expect a five times larger Ly α forest sample than SDSS, providing an unprecedented precise one-dimensional power spectrum measurement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DESI DR1 Ly α 1D power spectrum: the Fast Fourier Transform estimator measurement

Here, we present the one-dimensional Lyman-α forest power spectrum measurement derived from the data release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI). The measurement of the Lyman-α forest power spectrum along the line of sight from high-redshift quasar spectra provides information on the shape of the linear matter power spectrum, neutrino masses, and the properties of dark matter. In this work, we use a Fast Fourier Transform (FFT)-based estimator, which is validated on synthetic data in a companion paper. Compared to the FFT measurement performed on the DESI early data release, we improve the noise characterization with a cross-exposure estimator and test the robustness of our measurement using various data splits. We also refine the estimation of the uncertainties and now present an estimator for the covariance matrix of the measurement. Furthermore, we compare our results to previous high-resolution and eBOSS measurements. In another companion paper, we present the same DR1 measurement using the Quadratic Maximum Likelihood Estimator (QMLE). These two measurements are consistent with each other and constitute the most precise one-dimensional power spectrum measurement to date, while being in good agreement with results from the DESI early data release.

Lyman alpha forest↗
NASA NTRS1987

Analysis of eccentric annular incompressible seals. I - A new solution using fast Fourier transforms for determining hydrodynamic force

A new analysis procedure is presented which solves for the flow variables of an incompressible-flow annular pressure seal in which the rotor has a large static displacement from the centered position. The analysis begins with a set of governing equations based on a turbulent bulk-flow model and Moody's friction equation. No simplification of these bulk-flow equations is required for the solution procedure. Perturbation of the flow variables yields a set of zeroth and first-order equations. The zeroth-order equations (which model the large static displacement) are integrated by means of an efficient new method which employs fast Fourier transforms. Further integration of the zeroth-order pressures yields the hydrodynamic reactive force. Predictions for the hydrodynamic forces from this analysis procedure are in excellent agreement with available experimental results.

Nelson, C. C.↗
DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

tomoCAM : fast model-based iterative reconstruction via GPU acceleration and non-uniform fast Fourier transforms

X-ray-based computed tomography is a well established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and discoveries, with improved capabilities such as faster frame rates, larger fields of view, higher resolution and higher dimensionality. These advancements have enabled the material science community to expand the scope of tomographic measurements towards increasingly in situ and in operando measurements. In these new experiments, samples can be rapidly evolving, have complex geometries and restrictions on the field of view, limiting the number of projections that can be collected. In such cases, standard filtered back-projection often results in poor quality reconstructions. Iterative reconstruction algorithms, such as model-based iterative reconstructions (MBIR), have demonstrated considerable success in producing high-quality reconstructions under such restrictions, but typically require high-performance computing resources with hundreds of compute nodes to solve the problem in a reasonable time. Here, tomoCAM , is introduced, a new GPU-accelerated implementation of model-based iterative reconstruction that leverages non-uniform fast Fourier transforms to efficiently compute Radon and back-projection operators and asynchronous memory transfers to maximize the throughput to the GPU memory. The resulting code is significantly faster than traditional MBIR codes and delivers the reconstructive improvement offered by MBIR with affordable computing time and resources. tomoCAM has a Python front-end, allowing access from Jupyter -based frameworks, providing straightforward integration into existing workflows at synchrotron facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

NEUTRON-IRRADIATION–INDUCED DAMAGE ASSESSMENT IN CONCRETE USING COMBINED PHASE CHARACTERIZATION AND NONLINEAR FAST FOURIER TRANSFORM SIMULATION

Understanding irradiation effects in concrete cavities of light water reactors is critical to the long-term operation (60+ years) of US nuclear power plants. Susceptibility of concrete to neutron irradiation greatly varies as a function of its constituents: coarse aggregates, sand, and hardened cement paste (hcp). In particular, higher irradiation–susceptibility of aggregates was found to be a direct function of radiation-induced volumetric expansion (RIVE), or the propensity of aggregates to swell as a function of mineral contents, structures and textures. Irradiation-induced amorphization is accompanied by significant swelling, especially in silicates. For example, the maximum volumetric expansions of quartz and feldspars have been shown to be as large as ≈ 18% and ≈ 8%, respectively, while the change of density in calcite remains rather low (≈ 0.3–0.5%). High swelling is associated with a high number of covalent bonds (Si–O) and a high degree of polymerization of [SiO4] tetrahedrons. Differential RIVEs between rock-forming minerals cause high stresses and may potentially develop cracking within the aggregates at the interface between aggregates and the hcp and within the hcp. Modeling neutron-irradiation effects in concrete requires addressing the problem at a scale that captures both the details of the aggregates’ mineralogy and texture, which is <~ 100 μm, and the representative size of concrete, which is ~ 10 cm. This requires high-resolution, large mapping capabilities. To this end, the US Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program developed an original tool suite called Microstructure Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) in which Nonlinear 2D-fast Fourier transform (FFT) techniques are employed to model mechanical or irradiation-induced cracking of high-resolution complex microstructures derived from combined mapping techniques, including micro x-ray fluorescence (m-XRF), energy-dispersive x-ray spectroscopy (EDS), electron backscatter diffraction (EBSD), ellipsometry, and scanning electron microscopy (SEM). Irradiated mineral properties are derived from the Irradiated Minerals, Aggregate and Concrete (IMAC) database, which includes available relevant literature data. Rigorous coupled damage-creep simulations are obtained through sequential linear analysis. Illustrations of actual accelerated irradiation-damage propagation in a complex microstructures are presented herein.

Le Pape, Yann↗
DOE OSTI2018

Polycrystal Plasticity Models Based on Green’s Functions: Mean-Field Self-Consistent and Full-Field Fast Fourier Transform Formulations

We review two crystal plasticity methodologies based on Green’s functions for the prediction of microstructure-property relations in polycrystalline aggregates. The first, known as the viscoplastic self-consistent (VPSC) formulation, is a mean-field theory. The second, known as the viscoplastic fast Fourier transform-based (VPFFT) formulation, is a full-field method. Finally, the assumptions and main equations of these Green’s function-based crystal plasticity formulations are presented using a unified notation, pointing out their similarities and differences and cross-validating their predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2026

Fast Fourier transform evaluation of the Fresnel integral for gravitational-wave lensing

Gravitational waves (GWs) exhibit wave-optics effects when their wavelength is comparable to the scale of the gravitational lens. This may occur in lensing from galactic subhalos in GWs emitted by binary black-hole mergers and is gaining interest as a novel probe of dark matter. Predictions for observables in these cases ultimately rely on evaluating a Fresnel integral that quantifies the effect of lensing on the amplitude of a GW at a given frequency. However, numerical evaluation of this Fresnel integral is tricky, and several algorithms and publicly available codes that implement it have been developed. Here, we show that the dependence of this integral on the lens position can be written as a two-dimensional Fourier transform. Modern FFT techniques then enable rapid evaluation at all-sky positions simultaneously for general lenses without symmetry. Vectorization of FFT routines allows for derivatives with respect to model parameters to be obtained with only incremental additional computational cost. If the lens is axisymmetric, further speedups can be achieved with recently developed techniques for nonuniform fast Hankel transforms. To demonstrate, we make available Fresnel Integral Optimization with Nonuniform Transforms (fiona), an efficient and accurate code that is significantly faster than current methods for dense source grids, reaching 2 orders of magnitude speedups for ∼10 6 GW-emitting points. As part of FIONA , we developed code that provides vectorized nonuniform fast Hankel transforms that may have other uses (e.g., calculation of cosmological two-point correlation functions) beyond those considered here.

dark matter↗
DOE OSTI2016

Numerical implementation of non-local polycrystal plasticity using fast Fourier transforms

Here, we present the numerical implementation of a non-local polycrystal plasticity theory using the FFT-based formulation of Suquet and co-workers. Gurtin (2002) non-local formulation, with geometry changes neglected, has been incorporated in the EVP-FFT algorithm of Lebensohn et al. (2012). Numerical procedures for the accurate estimation of higher order derivatives of micromechanical fields, required for feedback into single crystal constitutive relations, are identified and applied. A simple case of a periodic laminate made of two fcc crystals with different plastic properties is first used to assess the soundness and numerical stability of the proposed algorithm and to study the influence of different model parameters on the predictions of the non-local model. Different behaviors at grain boundaries are explored, and the one consistent with the micro-clamped condition gives the most pronounced size effect. The formulation is applied next to 3-D fcc polycrystals, illustrating the possibilities offered by the proposed numerical scheme to analyze the mechanical response of polycrystalline aggregates in three dimensions accounting for size dependence arising from plastic strain gradients with reasonable computing times.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Benchmarking a Proof-of-Concept Performance Portable SYCL-based Fast Fourier Transformation Library

ABSTRACT In this paper, we present an early version of a SYCL-based FFT library, capable of running on all major vendor hardware, including CPUs and GPUs from AMD, ARM, Intel and NVIDIA. The current limitations of our library is it supports single-dimension FFTs up to 211 in length and base-2 input sequences. Although preliminary, the aim of this work is to seed further developments for a rich set of features for calculating FFTs. The library has the advantage over existing portable FFT libraries in that it is single-source, and there- fore removes the complexities that arise due to abundant use of pre-processor macros and auto-generated kernels to target different architectures. We exercise two SYCL-enabled compilers, Codeplay ComputeCpp and Intel's open-source LLVM project, to evaluate performance portability of our SYCL-based FFT on various hetero- geneous architectures.We provide studies comparing our portable library with highly optimized vendor-specific FFT libraries, and discuss potential sources hindering performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Development of Fast Fourier Transform (FFT) micro-mechanical simulations of concrete specimens characterized by micro-X-ray fluorescence

Concrete in Nuclear Power Plants (NPPs) can be exposed to a wide range of degradation phenomena. In the past years, the Light Water Reactor Sustainability (LWRS) program has investigated Radiation-Induced Volumetric Expansion (RIVE) as a potential degradation mechanism for concrete biological shields [Graves et al., 2014, Rosseel et al., 2016]. RIVE causes swelling and micro-mechanical damage in concrete due to the amorphization of mineral phases contained in the aggregates under neutron irradiation [Hilsdorf et al., 1978, Rosseel et al., 2016]. For long-term operations, it is critical to assess the durability of concrete after 60 or 80 years of exposure to NPP operating conditions against this phenomenon. RIVE is dependent on the composition of the aggregates used in concrete. Quartz-bearing aggregates are more sensitive to RIVE than calcite-bearing aggregates, for example. However, the aggregate composition of a specific plant is generally not explicitly given in the concrete formulation, which makes it nearly impossible to predict the resistance of that concrete to RIVE. Additional characterization is needed to identify the radiation-sensitive mineral phases contained in the aggregates.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Discrete orthogonal function expansions for non-uniform grids using the fast Fourier transform

A technique for applying discrete Fourier series to infinite domains is presented. The technique uses mappings designed to minimize truncation error and can be applied to solve mixed initial boundary value problems among others. The method is alias-free and yields consistent differentiation and integration operators. The mapping-induced truncation error is explicitly expressible and small in nearly all cases of interest. The method is illustrated for three problems involving convection, diffusion, and vortex interaction.

Cain, A. B.↗
NASA NTRS1987

SIR-B ocean-wave enhancement with fast Fourier transform techniques

Shuttle Imaging Radar (SIR-B) imagery is Fourier filtered to remove the estimated system-transfer function, reduce speckle noise, and produce ocean scenes with a gray scale that is proportional to wave height. The SIR-B system response to speckled scenes of uniform surfaces yields an estimate of the stationary wavenumber response of the imaging radar, modeled by the 15 even terms of an eighth-order two-dimensional polynomial. Speckle can also be used to estimate the dynamic wavenumber response of the system due to surface motion during the aperture synthesis period, modeled with a single adaptive parameter describing an exponential correlation along track. A Fourier filter can then be devised to correct for the wavenumber response of the remote sensor and scene correlation, with subsequent subtraction of an estimate of the speckle noise component. A linearized velocity bunching model, combined with a surface tilt and hydrodynamic model, is incorporated in the Fourier filter to derive estimates of wave height from the radar intensities corresponding to individual picture elements.

Tilley, David G.↗
NASA NTRS1991

Ordered fast Fourier transforms on a massively parallel hypercube multiprocessor

The present evaluation of alternative, massively parallel hypercube processor-applicable designs for ordered radix-2 decimation-in-frequency FFT algorithms gives attention to the reduction of computation time-dominating communication. A combination of the order and computational phases of the FFT is accordingly employed, in conjunction with sequence-to-processor maps which reduce communication. Two orderings, 'standard' and 'cyclic', in which the order of the transform is the same as that of the input sequence, can be implemented with ease on the Connection Machine (where orderings are determined by geometries and priorities. A parallel method for trigonometric coefficient computation is presented which does not employ trigonometric functions or interprocessor communication.

Tong, Charles↗