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 235 records · Page 13

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗

Collisional Histories of Comets and Trojan Asteroids: Diopside, Magnesite, and Fayalite Impact Studies

Comets and asteroids have weathered dynamic histories, as evidenced by their rough surfaces. The Nice model describes a violent reshuffling of small bodies during the Late Heavy Bombardment, with collisions acting to grind these planetesimals away. This creates an additional source of impact material that can re-work the surfaces of the larger bodies over the lifetime of the solar system. Here, we investigate the possibility that signatures due to impacts (e.g. from micrometeoroids or meteoroids) could be detected in their spectra, and how that can be explained by the physical manifestation of shock in the crystalline structure of minerals. All impact experiments were conducted in the Johnson Space Center Experimental Impact Laboratory using the vertical gun. Impact speeds ranged from approx.2.0 km/s to approx.2.8 km/s. All experiments were conducted at room temperature. Minerals found in comets and asteroids were chosen as targets, including diopside (MgCaSi2O6, monoclinic pyroxene), magnesite (MgCO3, carbonate), and fayalite (FeSiO4, olivine). Impacted samples were analyzed using a Fourier Transform Infrared Spectrometer (FTIR) and a Transmission Electron Microscope (TEM). Absorbance features in the 8-13 m spectral region demonstrate relative amplitude changes as well as wavelength shifts. Corresponding TEM images exhibit planar shock dislocations in the crystalline structure, attributed to deformation at high strain and low temperatures. Elongating or shortening the axes of the crystalline structure of forsterite (Mg2SiO4, olivine) using a discrete dipole approximation model (Lindsay et al., submitted) yields changes in spectral features similar to those observed in our impacted laboratory minerals.

Lederer, S. M.↗

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↗

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↗

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↗

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗

International Symposium on Remote Sensing of Environment, 10th, University of Michigan, Ann Arbor, Mich., October 6-10, 1975, Proceedings. Volumes 1 & 2

Topics treated include the application of a Fourier transform spectrometer to infrared remote sensing, the performance optimization of a satellite-borne thematic mapper, a data handling system to be integrated with a digital airborne multispectral scanner, infrared thermography for micro- and mesometeorological measurements, satellite interrogated data collection platforms for river and flood forecasting and the automatic measurement of sea surface temperature from a GOES satellite. Solar and atmospheric effects on satellite imagery derived from aircraft reflectance measurements, methods for determining haze levels from multispectral scanner data, restoration of Landsat images by discrete two-dimensional deconvolution and the automatic classification of aircraft and satellite multispectral images using mixed integer programming are also discussed. Individual items are announced in this issue.

Source record↗

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↗

Program for the analysis of time series

A digital computer program for the Fourier analysis of discrete time data is described. The program was designed to handle multiple channels of digitized data on general purpose computer systems. It is written, primarily, in a version of FORTRAN 2 currently in use on CDC 6000 series computers. Some small portions are written in CDC COMPASS, an assembler level code. However, functional descriptions of these portions are provided so that the program may be adapted for use on any facility possessing a FORTRAN compiler and random-access capability. Properly formatted digital data are windowed and analyzed by means of a fast Fourier transform algorithm to generate the following functions: (1) auto and/or cross power spectra, (2) autocorrelations and/or cross correlations, (3) Fourier coefficients, (4) coherence functions, (5) transfer functions, and (6) histograms.

Brown, T. J.↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

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↗

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↗

The stress intensity factor for the double cantilever beam

Fourier transforms and the Wiener-Hopf technique are used in conjunction with plane elastostatics to examine the singular crack tip stress field in the double cantilever beam (DCB) specimen. In place of the Dirac delta function, a family of functions which duplicates the important features of the concentrated forces without introducing unmanageable mathematical complexities is used as a loading function. With terms of order h-squared/a-squared retained in the series expansion, the dimensionless stress intensity factor is found to be K (h to the 1/2)/P = 12 to the 1/2 (a/h + 0.6728 + 0.0377 h-squared/a-squared), in which P is the magnitude of the concentrated forces per unit thickness, a is the distance from the crack tip to the points of load application, and h is the height of each cantilever beam. The result is similar to that obtained by Gross and Srawley by fitting a line to discrete results from their boundary collocation analysis.

Fichter, W. B.↗

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↗

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

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↗

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↗