Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

Toward Mixed Analog-Digital Quantum Signal Processing: Quantum AD/DA Conversion and the Fourier Transform

Signal processing stands as a pillar of classical computation and modern information technology, applicable to both analog and digital signals. Recently, advancements in quantum information science have suggested that quantum signal processing (QSP) can enable more powerful signal processing capabilities. However, the developments in QSP have primarily leveraged digital quantum resources, such as discrete-variable (DV) systems like qubits, rather than analog quantum resources, such as continuous-variable (CV) systems like quantum oscillators. Consequently, there remains a gap in understanding how signal processing can be performed on hybrid CV-DV quantum computers. Here we address this gap by developing a new paradigm of mixed analog-digital QSP. We demonstrate the utility of this paradigm by showcasing how it naturally enables analog-digital conversion of quantum signals—specifically, the transfer of states between DV and CV quantum systems. We then show that such quantum analog-digital conversion enables new implementations of quantum algorithms on CV-DV hardware. This is exemplified by realizing the quantum Fourier transform of a state encoded on qubits via the free-evolution of a quantum oscillator, albeit with a runtime exponential in the number of qubits due to information theoretic arguments. Collectively, this work marks a significant step forward in hybrid CV-DV quantum computation, providing a foundation for scalable analog-digital signal processing on quantum processors.

42 ENGINEERING↗
DOE OSTI2025

Broadband Fourier-Transform Optical Photothermal Infrared Spectroscopy and Imaging

Infrared (IR) spectroscopy is a powerful method for mapping chemical heterogeneity on the microscale. Synchrotron IR radiation uniquely provides a high brightness and broad bandwidth to further extend the capabilities of IR spectroscopic imaging. However, the diffraction-limited spatial resolution of IR spectroscopy is insufficient for studies requiring submicrometer spatial differentiation. Optical photothermal IR (O-PTIR) microscopy is a powerful, emerging method that overcomes the IR diffraction limit in IR hyperspectral imaging by employing a modulated IR beam and a visible probe laser beam to detect local temperature-induced modulation at the visible diffraction limit. In this work, we extend the spectral range of photothermal infrared measurements by incorporating a synchrotron IR source, demonstrating a combined synchrotron-based O-PTIR modality that enables high spatial resolution far-field chemical imaging spanning the entire mid-IR range. Both optical- and fluorescence-detected photothermal modalities were performed using a step-scan interferometer, demonstrating improved spectral range (541-4000 cm-1) when compared to optical photothermal microscopy with commercial laser sources (800-1800 cm-1 for this particular source) and improved spatial resolution, when compared to synchrotron microspectroscopy measurements. Following these initial validation studies, synchrotron Fourier-transform fluorescence-detected photothermal IR spectroscopy in combination with synchrotron microspectroscopy measurements was used to differentiate cells in mouse brain tissue sections, which requires submicron spatial resolutions beyond those accessible by IR spectroscopy alone.

Razumtcev, Aleksandr↗
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 OSTI2025

Wave-Optics Simulation Framework of Fourier Transform Holography with XMCD

We present a wave-optics simulation framework, implemented within the SRW, for modeling FTH with XMCD contrast. The framework propagates coherent, circularly polarized X-ray wavefronts from an undulator source through the sample and optics to the detector, reconstructing real-space images from the resulting holograms. Three principal extensions are introduced: (i) incorporation of polarization-dependent refractive indices in the sample plane, enabling direct simulation of XMCD contrast, (ii) a memory-efficient wavefront-splitting propagation scheme that treats individual apertures independently and coherently combines their fields at detector, and (iii) an integrated reconstruction module that delivers holographic images within the same framework. Together, these advances establish a versatile tool for quantitative exploration of coherence, aperture geometry, and detector sampling, and for the design and interpretation of XMCD-FTH experiments at synchrotron beamlines.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Quantum Fourier Transform with Leakage Control

In circuit QED (cQED) systems, one can manipulate the cavity Fock states to use quantum processing unit (QPU). Since the coherence times of the Fock states are long, qudit based computation can be possible. One of the most basic gate that can be engineered in this cQED platform is the displacement gate, $D(\alpha)$. However, accessible number of qudit-levels is finite. The truncated displacement gate on the qudit is the source of theoretical error above the logical states and needs to be addressed. In this work, we monitor and minimize the leakage in qudit systems that involves displacement gate by introducing so-called `bumper' states and a modified cost function.

Kurkcuoglu, D. [Fermilab] (ORCID:0000000311097074)↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2026

Jahn–Teller distortion controls electron transfer in photoexcited Cu( I ) donor–acceptor systems

The Jahn–Teller distortion (JTD) is a defining structural response to electronic excitation in Cu( I )-based transition metal complexes, yet its role in photoinduced electron transfer (PET) remains largely unexplored. Here, we demonstrate that the JTD governs charge separation (CS) through a vibronically controlled conical intersection in heteroleptic Cu( I ) bisphenanthroline–naphthalene diimide (CuHETPHEN-NDI) donor–acceptor dyads. Using 20-fs broadband transient absorption spectroscopy combined with coherent vibrational wavepacket (CVWP) analysis and quantum chemical calculations, we directly track nuclear motions that steer the system from the metal-to-ligand charge-transfer ( 1 MLCT) state to the CS state. Steric bulkiness of the pendant groups at the 2,9-positions of the phenanthroline ligand systematically slows both JTD and CS. Short-time Fourier transformation and Fourier filtering analyses identify two key vibrational signatures: a low-frequency breathing mode (∼100 cm −1 ) via a bond distance change between Cu and ligated Ns (Cu–N) that modulates the NDI anion absorption and acts as a vibronic coupling coordinate, and a higher-frequency mode (∼313 cm −1 ) that evolves along the PET trajectory. Normal mode analysis and potential energy surface calculations show that the JTD brings the 1 MLCT and CS states into degeneracy, while the Cu–N breathing motion dynamically modulates donor–acceptor electronic coupling to enable ultrafast nonadiabatic electron transfer. Steric hindrance exerted by the groups at the 2,9 positions of the phenanthroline ligands suppresses this vibronic coupling, leading to faster CVWP decoherence for the 313 cm −1 mode and slower CS. These findings unravel JTD-controlled vibronic coupling at conical intersection as a governing factor for CS and provide insight into designing Cu-based photosensitizers by harnessing structural dynamics to control PET.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗
DOE OSTI2026

Second Order System Study

During my education in mathematics, engineering and physics, I learned transform pairs and their usage mechanics but I never remember seeing the derivations of the solutions to second order ordinary differential equations (ODE) and difference equations. A solution to a question posed in a potential funder meeting put me on a path to solving second order systems using the five principal Fourier based methods: Fourier transform (FT), Z-transform (ZT), discrete time Fourier transform (DTFT), discrete Fourier transform (DFT) and Laplace transform (LT).

42 ENGINEERING↗
DOE OSTI2024

Estimation of Forest Aboveground Biomass from Derivatives of Vegetation-Structure Profiles

Several studies have found that the vertical Fourier transform of lidar, interferometric Synthetic Aperture Radar (SAR), and stereo photogrammetric profiles at empirically-determined spatial frequencies enables high-performance forest aboveground biomass (AGB) estimation. Linear combinations of real and imaginary parts of Fourier transforms of Tomographic (multi-baseline) SAR (TomoSAR) profiles, from Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne data, generate ~20%-precision estimates of AGB in the Saskatchewan area of Canada. We found that this 20% precision can be improved to ~15%, a factor of 30% improvement in root mean square error (RMSE) if, in addition to using Fourier transforms of the profile itself, we use Fourier transforms of the spatial, vertical derivative of the profile. The formulation of this "derivative" algorithm is the subject of this paper.

Treuhaft, Robert↗
DOE OSTI2026

Composition and persistence of soil organic matter along eroding and depositional transects in buried vs. modern soil layers: A case of the Brady paleosol at Wauneta, Nebraska

Paleosols form when soils are buried through deposition by aeolian, colluvial, alluvial or other processes. Burial of former topsoil isolates soil organic matter (SOM) from surface conditions, allowing carbon to accumulate and potentially remain stable for millennia. In this study, SOM composition, distribution, and persistence were analyzed in the Brady Soil of Nebraska, USA to compare SOM spatial variability in modern and buried soils, as well as the impact of erosional exposure on SOM stability. The Brady Soil, formed as a surface soil during the Pleistocene-Holocene transition and now a paleosol buried up to 6 m deep (or more) by loess deposition during the Holocene, was sampled along burial (up to 5.8 m depth) and erosional (up to 1.8 m depth) transects to compare SOM dynamics in different geomorphic settings. Fourier Transform Infrared Spectroscopy (FTIR) and Fourier Transform ion cyclotron resonance mass spectrometry (FTICR-MS) were used to analyze SOM composition, while δ 13 C isotope analyses identified SOM sources and radiocarbon values were used to estimate turnover rates. Results confirmed a vegetation shift from C3 to C4 plants after Brady Soil formation, reflecting warming climatic conditions. Increasing SOM age and decreasing δ 13 C and δ 15 N values with depth indicated slowing of decomposition rate in buried soils. Higher pH in the Brady Soil suggested greater base cation content, supporting SOM stabilization through organo-mineral associations and aggregate formation. However, exposure of the Brady Soil due to surface erosion caused faster SOM turnover. This result suggested susceptibility of buried SOM to losses via decomposition upon erosional exposure, possibly accelerated by priming in response to modern SOM inputs. These findings highlight the potential loss of carbon stocks in buried soils under future climate change, as shifts in soil physicochemical properties may destabilize long-preserved SOM.

Environmental sciences↗
DOE OSTI2026

Extensive analysis of reconstruction algorithms for DESI 2024 baryon acoustic oscillations

Reconstruction of the baryon acoustic oscillation (BAO) signal has been a standard procedure in BAO analyses over the past decade and has helped to improve the BAO parameter precision by a factor of ∼2 on average. The Dark Energy Spectroscopic Instrument (DESI) BAO analysis for the first year (DR1) data uses the “standard” reconstruction framework, in which the displacement field is estimated from the observed density field by solving the linearized continuity equation in redshift space, and galaxy and random positions are shifted in order to partially remove non-linearities. There are several approaches to solving for the displacement field in real survey data, including the multigrid (MG), iterative Fast Fourier Transform (iFFT), and iterative Fast Fourier Transform particle (iFFTP) algorithms. In this work, we analyze these algorithms and compare them with various metrics including two-point statistics and the displacement itself using realistic DESI mocks. We focus on three representative DESI samples, the emission line galaxies (ELG), quasars (QSO), and the bright galaxy sample (BGS), which cover the extreme redshifts and number densities, and potential wide-angle effects. We conclude that the MG and iFFT algorithms agree within 0.4% in post-reconstruction power spectrum on BAO scales with the RecSym convention, which does not remove large-scale redshift space distortions (RSDs), in all three tracers. The RecSym convention appears to be less sensitive to displacement errors than the RecIso convention, which attempts to remove large-scale RSDs. However, iFFTP deviates from the first two; thus, we recommend against using iFFTP without further development. In addition, we provide the optimal settings for reconstruction for five years of DESI observation. The analyses presented in this work pave the way for DESI DR1 analysis as well as future BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Discrete Hankel Transform Approach to Nuclear Data Processing for Fusion Applications

This study introduces advancements to the numerical solutions employed in the processing of nuclear data for fusion applications. It leverages the convolution theorem and Fourier transform techniques to enhance computational efficiency and broaden applicability. Building upon a previously reported discrete Hankel transform approach for Doppler broadening, this work refines the solution of convolution integrals central to these applications. The methodology provides a general and unified framework for evaluating any convolution operation, regardless of whether the underlying problem involves temperature effects in nuclear reactions. The applicability to the nuclear data processing for fusion is demonstrated by deriving the convolution integrals for some of the fusion-related quantities. As before, the convolution operation utilizes a Gaussian-based kernel; however, the discrete Hankel transform of order $𝛼$ = $\frac{1}{2}$ is now applied to the forward Fourier transform of the nonkernel argument, rather than the inverse Fourier transform. This modification eliminates the need for the integration of the nonkernel, cross section–based function, which is a step that posed challenges for certain pointwise cross-section representations. It also removes the requirement for cross-section linearization. Optimized for graphics processing unit architectures, the approach significantly improves computational performance. These advancements are currently under evaluation as the foundation for the next-generation thermonuclear data file processing codes being developed at Lawrence Livermore National Laboratory.

Nuclear science and engineering↗
DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗
DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗