Search NASA⌕ Search

SEARCH · Search NASA

Results for “fourier expansion”

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 163 records · Page 9

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N -body mock challenge for the quasar sample

ABSTRACT The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock–Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of halo occupation distribution mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover fσ8 to within 3 per cent and α∥ and α⊥ to within 1 per cent. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on fσ8, α∥, and α⊥ at the level of $\sigma _{f\sigma _8}=0.013$, $\sigma _{\alpha _\parallel }=0.012$, and $\sigma _{\alpha _\bot }=0.008$. The systematic error on the combined consensus is $\sigma _{f\sigma _8}=0.011$, $\sigma _{\alpha _\parallel }=0.008$, and $\sigma _{\alpha _\bot }=0.005$, which is used in the final DR16 analysis. For baryon acoustic oscillation fits in configuration and Fourier space, we take conservative systematic errors of $\sigma _{\alpha _\parallel }=0.010$ and $\sigma _{\alpha _\bot }=0.007$.

79 ASTRONOMY AND ASTROPHYSICS↗

Reducing heat load density with asymmetric and inclined double-crystal monochromators: principles and requirements revisited

Asymmetric double-crystal monochromators (aDCMs) and inclined DCMs (iDCMs) can significantly expand the X-ray beam footprint and consequently reduce the heat load density and gradient. Based on rigorous dynamical theory calculations, the major principles and properties of aDCMs and iDCMs are presented to guide their design and development, particularly for fourth-generation synchrotrons. In addition to the large beam footprint, aDCMs have very large bandwidths (up to ∼10 eV) and angular acceptance, but the narrow angular acceptance of the second crystal requires precise control of the relative orientations and strains. Based on Fourier coupled-wave diffraction theory calculations, it is rigorously proved that the iDCM has almost the same properties as the conventional symmetric DCM, including the efficiency, angular acceptance, bandwidth, tuning energy range and sensitivity to misalignment. The exception is that, for the extremely inclined geometry that can achieve very large footprint expansion, the iDCM has (beneficially) a larger bandwidth and wider angular acceptance. Inclined diffraction has the `rho-kick effect' that can be cancelled by the second reflection of the iDCM (even with misalignment), except that inhomogeneous strains may cause non-uniform rho-kick angles. At present, fabrication/mounting-induced strains pose low risk since they can be controlled to <0.5 µrad over large areas. The only uncertain challenge is the thermally induced strains, yet it is estimated that these strains are naturally lowered by the large footprint and may be further mitigated by optimized cryogenic cooling to the 1–2 µrad level. Overall, aDCMs and iDCMs have more stringent requirements than normal DCMs, but they are feasible schemes in practice.

asymmetric monochromator↗

Decomposition of multi-particle azimuthal correlations in Q-cumulant analysis*

Abstract The method of Q -cumulants is a powerful tool for studying the fine details of azimuthal anisotropies in high energy nuclear collisions. This paper presents a new method, based on mathematical induction, to evaluate the analytical form of high-order Q -cumulants. The capability of this method is demonstrated via a toy model that uses the elliptic power distribution to simulate the anisotropic emission of particles, quantified in terms of Fourier flow harmonics . The method can help in studying the large amount of event statistics that can be collected in the future and allow measurements of the very high central moments of the distribution. This can, in turn, facilitate progress in understanding the initial geometry, the input to the hydrodynamic calculations of medium expansion in high energy nuclear collisions, and the constraints on it.

Physics↗

Identifying approximate linear models for simple nonlinear systems

This paper addresses the identification (realization) of approximate linear models from response data for certain nonlinear dynamic systems. Response characteristics for several typical nonlinear joints are analyzed mathematically and represented by series expansions. The parameters of the series expansion are then compared with the modal parameters of a linear model identified by the Eigensystem Realization Algorithm. The agreement of the identified model and the analytically derived representation is excellent for the cases studied. Also laboratory data from a model which exhibited stiffening behavior was analyzed using the Eigensystem Realization algorithm and Fast Fourier Transform. The laboratory experiment demonstrated the ability of the technique to recover the model characteristics using real data.

Horta, L. G.↗

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of load introduction on graphite epoxy compression specimens

Compression testing of modern composite materials is affected by the manner in which the compressive load is introduced. Two such effects are investigated: (1) the constrained edge effect which prevents transverse expansion and is common to all compression testing in which the specimen is gripped in the fixture; and (2) nonuniform gripping which induces bending into the specimen. An analytical model capable of quantifying these foregoing effects was developed which is based upon the principle of minimum complementary energy. For pure compression, the stresses are approximated by Fourier series. For pure bending, the stresses are approximated by Legendre polynomials.

Reiss, R.↗

The direction of heat flow in an expanding current free stellar atmosphere

The energy budgets of the corona are substantially influenced by losses attributed to thermal conduction losses down the transition region, towards the chromosphere. Although Fourier law for such a plasma regime is now universally criticized as inappropriate, the direction of this heat flow seems to be certain. The purpose of this paper is to point out, based on conservation laws, the direction in which the heat must flow in this regime, provided it is known that there is a time independent expansion devoid of parallel currents. It will be argued that the direction of the heat flux is from the base of the transition region up into the corona, rather than the opposite. Under the circumstances mentioned above, the direction of the heat flow should be opposite to the gradient of the electron pressure gradient. Although the transition region is modeled as isobaric, this is only approximately true. Such a likely scenario makes the heat flow almost certainly to flow from the lower temperature base of the TR up into the incipient corona if a steady state current free expansion is to be maintained.

Scudder, J. D.↗

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

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

Hypersonics↗

NO PLIF Rotational Thermometry Characterization of a Hypersonic Boundary Layer

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

Hypersonics↗

Raman spectroscopy of gases with a Fourier transform spectrometer - The spectrum of D2

Fourier transform spectrometry (FTS) is presently used to record the spontaneous incoherent laser Raman spectra of gases. The high resolution, sensitivity, calibration accuracy and spectral coverage achieved demonstrate the viability of FTS for Raman spectroscopy. Attention is given to the coefficients obtained by fitting measurements obtained from the spectrum of D2, containing both v = 0-0 and 1-0 transitions, to the Dunham (1932) expansion of the vibration-rotation energy levels.

Jennings, D. E.↗

Kinetic theory of turbulent compressible flows and comparison with classical theory

A complete set of turbulent correlations is given in terms of expansion coefficients in a double series of Hermite polynomials of a two-particle correlation in phase space. Only two of these coefficients, corresponding to Reynolds stress and turbulent heat flux, are shown to appear in gasdynamic equations of turbulence. This is accomplished by incorporating the phase-space correlation in the collision integral of Boltzmann's equation and by deriving generalized Navier-Stokes and Fourier transport relations. In contrast to the classical formalism, the simple expression remains invariant in form whether the flow is compressible or not. Nevertheless, a bilinear transformation of fluctuating quantities shows that the two formalisms are identical with regard to terms of double correlations. The kinetic theory justifies nonexistence of the higher-order correlations characteristic of compressible-turbulence equations in the classical regime.

Tsuge, S.↗

Broadband Rotational Spectroscopy in Uniform Supersonic Flows: Chirped Pulse/Uniform Flow for Reaction Dynamics and Low Temperature Kinetics

ConspectusThe study of gas-phase chemical reactions at very low temperatures first became possible with the development and implementation of the CRESU (French acronym for Reaction Kinetics in Uniform Supersonic Flows) technique. CRESU relies on a uniform supersonic flow produced by expansion of a gas through a Laval (convergent-divergent) nozzle to produce a wall-less reactor at temperatures from 10 to 200 K and densities of 1016-1018 cm-3 for the study of low temperature kinetics, with particular application to astrochemistry. In recent years, we have combined uniform flows with revolutionary advances in broadband rotational spectroscopy to yield an instrument that affords near-universal detection for novel applications in photodissociation, reaction dynamics, and kinetics. This combination of uniform supersonic flows with chirped-pulse Fourier-transform microwave spectroscopy (Chirped-Pulse/Uniform Flow, CPUF) permits detection of any species with a modest dipole moment, thermalized to the uniform temperature of the gas flow, with isomer, conformer, and vibrational state specificity. In addition, the use of broadband, high-resolution, and time-dependent (microsecond time scale) micro- and mm-wave spectroscopy makes it an ideal tool for characterizing both transient and stable molecules, as well as studying their spectroscopy and dynamics.In this Account, we review recent advances made using the CPUF technique, including studies of photodissociation, radical-radical reaction dynamics, and low temperature kinetics. These studies highlight both the strength of universal and multiplexed detection and the challenges of coupling it to a high-density collisional environment. Product branching and product evolution as a function of time have been measured for astrochemically relevant systems, relying on the detailed characterization of these flow conditions via experiments and fluid dynamics simulations. In the photodissociation of isoxazole, an unusual heterocyclic molecule with a very low-energy conical intersection, we have identified 7 products in 5 reaction channels and determined the product branching, pointing to both direct and indirect pathways. We have also approached the same system from separated NO and C3H3 reactants to explore a broader range of the potential energy surface, demonstrating the power of multichannel branching measurements for complex radical-radical reactions. We determined the product branching in the C3H2 isomers in the photodissociation of the propargyl radical and identified the importance of a hydrogen atom catalyzed isomerization to the lowest energy cyclic form. This then motivated a study of direct D-H exchange reaction in radicals, in which we demonstrate that it is an important and overlooked pathway for deuterium fractionation in astrochemical environments. Recently, we have shown the measurement of low temperature kinetics inside an extended Laval nozzle, after which a shock-free secondary expansion to low temperature and density affords an ideal environment for detection by rotational spectroscopy. These results highlight the power and potential of the CPUF approach, and future prospects will also be discussed in light of these developments.

Dias, Nureshan↗

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗

High-resolution absorption spectroscopy of room-temperature and jet-cooled ammonia between 59,000 and 93,000 cm -1

Here, we present new high-resolution photoabsorption spectra of ammonia spanning the region between 59,000 cm -1 and 93,000 cm -1 that were recorded by using the Fourier Transform Spectrometer at the Synchrotron SOLEIL. This region extends from just above the Franck-Condon envelope for the $\tilde{A}$ 1 $A$ 2 " ← $\tilde{X}$ 1 $A$ 1 ' transition to well above the NH 3 + $\tilde{X}$ + 2 $A$ 2 " ionization threshold. The spectra were recorded at a measured resolution of 0.23 cm -1 in both a room-temperature cell (293 K) and in a slit-jet supersonic expansion (~70 K). The absolute photoabsorption cross section with an uncertainty of ±5% is also reported for the room-temperature spectrum. The present resolution is a factor of 10 – 100 times higher than in other recently reported broad band spectra of ammonia, and many of the observed bands show partially resolved rotational structure. We have attempted to assign this structure for a number of these bands. The oscillator strengths extracted from the data are in good agreement with previous measurements but, in the case of structured bands, the present higher resolution measurements return higher peak absorption cross sections, that increase further when the sample is cooled. The present higher resolution spectra suggest that a number of previous vibronic band assignments that were based on quantum defect considerations may require some revision. Finally, we discuss the substantial differences between the photoabsorption and photoionization data just above the first ionization threshold.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Rapidly-convergent flux-surface shape parameterization

We propose a novel flux-surface parameterization suitable for local MHD equilibrium calculations with strongly-shaped flux surfaces. The method is based on a systematic expansion in a small number of intuitive shape parameters, and reduces to the well-known Miller D-shaped parameterization in the limit where some of the coefficients are set to zero. Here, the new parameterization is valid for up-down asymmetric plasmas and provides an improvement to the Miller form. Simultaneously, the method is rapidly convergent and requires only about half the number of shape parameters as a general Fourier representation in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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.↗

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗