Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Discrete 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Primitive quantum gates for an S U ( 2 ) discrete subgroup: Binary tetrahedral

We construct a primitive gate set for the digital quantum simulation of the binary tetrahedral ($\mathbb{BT}$) group on two quantum architectures. This nonabelian discrete group serves as a crude approximation to $SU(2)$ lattice gauge theory while requiring five qubits or one quicosotetrit per gauge link. The necessary basic primitives are the inversion gate, the group multiplication gate, the trace gate, and the $\mathbb{BT}$ Fourier transform over $\mathbb{BT}$. We experimentally benchmark the inversion and trace gates on ibm nairobi, with estimated fidelities between $14-55\%$, depending on the input state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗
DOE OSTI2026

Quantum simulation of QED in Coulomb gauge

A recent work considered quantum simulation of Quantum Electrodynamics on a lattice in the Coulomb gauge with gauge degrees of freedom represented in the occupation basis in momentum space. In this work, we consider the more efficient representation of the gauge degrees of freedom in field basis in position space and develop a quantum algorithm for real-time simulation. We show that the continuum Coulomb gauge Hamiltonian is equivalent to the temporal gauge Hamiltonian when acting on physical states consisting of fermion and transverse gauge fields. The Coulomb gauge Hamiltonian is discretized by using the Green’s function of the discrete Laplacian operator under the Dirichlet boundary conditions. Both the continuum Coulomb gauge Hamiltonian and the discretized one proposed here guarantee that the unphysical longitudinal gauge fields are decoupled and commute with the corresponding Hamiltonian. Thus there is no need to impose any constraint. The local gauge field basis and the canonically conjugate variable basis are swapped efficiently using the quantum Fourier transform. We prove that the qubit cost to represent physical states and the gate count for real-time simulation scale polynomially with the lattice size, energy, time, accuracy, and Hamiltonian parameters in lattice units. The gate cost here for implementing the time evolution of the gauge field is reduced at least by a factor on the order of 10 8 for modest lattice size and accuracy level compared with the previous work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Primitive quantum gates for dihedral gauge theories

We describe the simulation of dihedral gauge theories on digital quantum computers. The nonabelian discrete gauge group $D_N$-- the dihedral group -- serves as an approximation to $U(1)\times\mathbb{Z}_2$ lattice gauge theory. In order to carry out such a lattice simulation, we detail the construction of efficient quantum circuits to realize basic primitives including the nonabelian Fourier transform over $D_N$, the trace operation, and the group multiplication and inversion operations. For each case the required quantum resources scale linearly or as low-degree polynomials in $n=\log N$. We experimentally benchmark our gates on the Rigetti Aspen-9 quantum processor for the case of $D_4$. The fidelity of all $D_4$ gates was found to exceed $80\%$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q 2 ∼ 4.3 GeV 2 and the nucleon up to Q 2 ∼ 1 GeV 2 . The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼ 250 to ∼ 540 MeV , including the unitary point ∼ 340 MeV . We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Primitive quantum gates for an S U ( 2 ) discrete subgroup: Binary octahedral

We construct a primitive gate set for the digital quantum simulation of the 48-element binary octahedral ( BO ) group. This non-Abelian discrete group better approximates S U ( 2 ) lattice gauge theory than previous work on the binary tetrahedral group at the cost of one additional qubit—for a total of six—per gauge link. The necessary primitives are the inversion gate, the group multiplication gate, the trace gate, and the BO Fourier transform. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Core Characterization of Bakken Shale from the Bedwell 33-52-1-1H Well, Sheridan County, Montana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core of the Bedwell 33-52-1-1H well (API 25-091-21920), a wildcat field in Sheridan County, Montana. The primary impetus of this work is to provide a unique dataset to researchers at universities, national laboratories, geological surveys, and other locations for additional analyses. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/bedwell-33-52-1-1h-well. The equipment and techniques used to characterize the full core were non-destructive, enabling future examinations and analyses to be performed on these cores. However, none of the equipment used was suitable for direct visualization of the pore space in fine-grained structures such as those found in shales; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. A selection of samples (Table 1) were cored, sliced, and powdered for additional characterization analyses. These analyses did involve destructive subsampling of the full core at discrete sections, and the results are presented here in the context of the larger core description. Selected sections of the core were removed from the full core for detailed analyses using the following destructive techniques: bulk inductively coupled plasma-optical emission spectroscopy (ICP-OES), isolation of carbon and sulfur, Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM), and higher resolution CT scan with NETL’s micro-CT scanning systems. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Primitive quantum gates for an S U ( 3 ) discrete subgroup: Σ ( 36 × 3 )

We construct the primitive gate set for the digital quantum simulation of the 108-element Σ ( 36 × 3 ) group. This is the first time a non-Abelian crystal-like subgroup of S U ( 3 ) has been constructed for quantum simulation. The gauge link registers and necessary primitives—the inversion gate, the group multiplication gate, the trace gate, and the Σ ( 36 × 3 ) Fourier transform—are presented for both an eight-qubit encoding and a heterogeneous three-qutrit plus two-qubit register. For the latter, a specialized compiler was developed for decomposing arbitrary unitaries onto this architecture. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spectral Radiative Properties of Polydispersed SiO 2 Particle Beds

The focus of this work is on the measurement and analysis of the radiative properties of polycrystalline SiO 2 particle beds with various layer thicknesses. The particles are polydispersed with average diameters of 222, 150, and 40 μm . The spectral, directional–hemispherical reflectance and transmittance of the particle bed are measured at wavelengths from 0.4 to 1.8 μm using a monochromator, and the reflectance measurement is extended to 15 μm using a Fourier-transform infrared spectrometer. Particles are closely packed between two transparent windows for measuring the radiative properties. In the visible and near-infrared region up to 1.8 μm, the inverse adding–doubling method yields the effective absorption and scattering coefficients. The results suggest that short wavelength absorption needs to be included in modeling the behavior of particle beds due to multiple scattering. A discrete-scale Monte Carlo ray-tracing method is developed to model the radiative properties by assuming monodispersed spherical particles, and the simulated results compare well with measurements. The effective absorption and scattering coefficients of the particle beds obtained from the independent scattering theory are compared to those from the inverse method. As a result, the impact of dependent scattering on the packed beds is observed for smaller-sized particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Interfacial Unit-Dependent Catalytic Activity for CO Oxidation over Cerium Oxysulfate Cluster Assemblies

Atomically precise cerium oxo clusters offer a platform to investigate structure–property relationships that are much more complex in the ill-defined bulk material cerium dioxide. We investigated the activity of the MCe 70 torus family (M = Cd, Ce, Co, Cu, Fe, Ni, and Zn), a family of discrete oxysulfate-based Ce 70 rings linked by monomeric cation units, for CO oxidation. CuCe 70 emerged as the best performing MCe 70 catalyst among those tested, prompting our exploration of the role of the interfacial unit on catalytic activity. Temperature-programmed reduction (TPR) studies of the catalysts indicated a lower temperature reduction in CuCe 70 as compared to CeCe 70 . In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) indicated that CuCe 70 exhibited a faster formation of Ce 3+ and contained CO bridging sites absent in CeCe 70 . Isothermal CO adsorption measurements demonstrated a greater uptake of CO by CuCe 70 as compared to CeCe 70 . The calculated energies for the formation of a single oxygen defect in the structure significantly decreased with the presence of Cu at the linkage site as opposed to Ce. Furthermore, this study revealed that atomic-level changes in the interfacial unit can change the reducibility, CO binding/uptake, and oxygen vacancy defect formation energetics in the MCe 70 family to thus tune their catalytic activity.

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

Quantification of active sites in yttrium containing dealuminated Beta zeolites during conversion of ethanol and acetaldehyde to butadiene

Here, in this work, yttrium containing dealuminated Beta zeolites (Y/deAlBeta) were synthesized and characterized by various spectroscopic techniques to improve understanding of ethanol upgrading over these materials. Characterization results indicate yttrium atoms partially condense with framework silanol nests formed during dealumination of parent Al-Beta supports. Active sites for conversion of ethanol and acetaldehyde to butadiene were quantified on a series of Y/deAlBeta catalysts (0.1–10 wt% yttrium) via ex situ chemisorption and transmission Fourier transformed infrared (FTIR) spectroscopy measurements by first measuring the integrated molar extinction coefficient (IMEC) for pyridine bound to Lewis acidic yttrium sites. In situ titrations with pyridine demonstrate that the number of sites quantified by ex situ chemisorption IR is quantitatively similar to the number of sites that catalyze butadiene formation, which varies (from 0.05 to 0.35) across the series of catalysts. In situ pyridine titrations impact butadiene site time yields (STY), but not crotonaldehyde STY, indicating that a distribution of yttrium sites is present, and that discrete yttrium site types participate in distinct steps in the pathway from ethanol to butadiene. Apparent kinetic parameters including activation energies and reaction orders were measured, these suggest differences in reactant (or reactant-derived intermediate) surface coverages result in higher STYs (per mol Y or per Lewis acidic Y site) for samples with low Y loadings relative to those with higher Y loadings. Isotopic labeling experiments evince the existence of other kinetically relevant steps in addition to the crotonaldehyde transformation to crotyl alcohol. Together, these findings provide further guidance into the heterogeneities in site structures in yttrium-containing zeolites and their relevance for the various steps in the pathway from ethanol to C 4 products useful for production of sustainable aviation fuel and renewable butadiene.

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

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Enhanced Catalytic Performance of a Ce/V Oxo Cluster through Confinement in Mesoporous SBA-15

To increase catalytic efficiency, mesoporous supports have been widely applied to immobilize well-defined metal oxide clusters due to their ability to stabilize highly dispersed clusters. Here, for this study, a redox-active heterometallic Ce 12 V 6 -oxo cluster (CeV) was first presynthesized and then incorporated into mesoporous silica, SBA-15, via a straightforward impregnation method. Scanning transmission electron microscopy (STEM) and Fourier transform infrared spectroscopy (FTIR), in concert with scanning electron microscopy and energy-dispersive X-ray spectroscopy (SEM-EDS), verified the successful introduction of the CeV cluster inside the pore of SBA-15. The 51 V magic angle spinning solid-state nuclear magnetic resonance ( 51 V MAS NMR) spectroscopy and differential pair distribution function (dPDF) analysis confirmed the structural integrity of the CeV cluster inside the SBA-15. The composite was then benchmarked for liquid-phase oxidation of 2-chloroethyl ethyl sulfide (CEES) under mild conditions and gas-phase oxidative dehydrogenation (ODH) of propane under high temperatures (up to 550 °C). The catalytic reactivity results demonstrated 8- and 14-fold increase in turnover frequency (TOF) values of the composite (CeV@10SBA-2) than the bulk CeV cluster under the same conditions for CEES oxidation and ODH, respectively. These results highlight the improved reactivity of the catalytically active CeV cluster as attributed to the higher dispersion of the discrete cluster upon immobilization within the SBA-15 support.

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

Evaluating the Durability of Balance of Systems Components Using Combined-Accelerated Stress Testing

The degradation of photovoltaic (PV) balance of systems (BoS) components is not well-studied, but the consequences include: offline-modules, -strings, -inverters; system shutdown; arc-faults; and fires. A utility provider experienced a ~30% failure rate in their power transfer chain, attributed to branch connectors. Field-failed specimen assemblies were therefore examined, consisting of cable connector, branch connector and discrete fuse components. Unused field-vintage specimens are presently being examined using combined-accelerated stress testing (C-AST) to clarify the most influential environmental stressors as well as the effect of external mechanical perturbation. A benchtop prototype fixture was used to develop the perturbation capability for C-AST. The benchtop experiments were also used to develop the in-situ data acquisition of specimen: current, voltage, and temperature. A significant increase in operating temperature, ~100 deg C from ~40 deg C, and a different failure mode was observed immediately once periodic mechanical perturbation was applied. The current at failure was decreased from 35 A (with failure occurring in the fuses) to 15 A (failure at the male/female metal pin connection). After initial examination using X-ray computed tomography, the external plastic was machined away from failed specimens to allow failure analysis, including the extraction of the internal convolute springs for morphological examination (optical- and electron-microscopy). Chemical composition analysis included: energy-dispersive X-ray spectroscopy, differential scanning calorimetry, and Fourier transform infrared spectroscopy.

balance of systems↗