Search NASA⌕ Search

SEARCH · Search NASA

Results for “Harmonic numbers”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Supergravity black holes, Love numbers, and harmonic coordinates

A problem which has now come to a head in general relativity, primarily as a result of the many recent detections of gravitational waves from the coalescence of compact binaries, is that there are very few viable theories against which to “test” it. To perform realistic tests of theories of gravity, we need to be able to look beyond general relativity and evaluate the consistency of a parametrized, physically acceptable, family of black hole metric alternatives with observational data from, especially, gravitational wave detections using, for example, an agnostic Bayesian approach. In this paper we further examine properties of one class of such metrics, which in fact arise as solutions of ungauged supergravity. In particular, we examine the massless, neutral, minimally coupled scalar wave equation in a general stationary, axisymmetric background metric such as that of a charged rotating black hole, when the scalar field is either time independent or in the low-frequency, near-zone limit, with a view to calculating the Love numbers of tidal perturbations, and of obtaining harmonic coordinates for the background metric. For a four-parameter family of charged asymptotically flat rotating black hole solutions of ungauged supergravity theory known as STU black holes, which includes Kaluza-Klein black holes and the Kerr-Sen black hole as special cases, we find that all time-independent solutions, and hence the harmonic coordinates of the metrics, are identical to those of the Kerr solution. In the low-frequency limit we find the scalar fields exhibit the same SL(2,R) symmetry as holds in the case of the Kerr solution. We point out extensions of our results to a wider class of metrics, which includes solutions of Einstein-Maxwell-dilaton theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Illustration of bunch merge simulation in RHIC

Following are examples of a 4 to 2 to 1 merge of proton bunches in RHIC obtained by running the simulation code rhic2mrg23. As the code runs, the user is prompted to enter several numbers. For these examples, RF harmonic number 1440.0 is entered. An RF frequency of 112.597 696 626 MHz at this harmonic is entered. This gives a proton energy of 100 GeV. The merge simulation starts with a uniform distribution of unbunched protons. 2.0 eV-s is entered for the longitudinal emittance of the distribution. The code will recommend a voltage to capture the unbunched protons into 4 harmonic 1440 buckets. A capture voltage of 142.390 kV is recommended. 144.0 kV is entered. It is desirable to capture the protons as “adiabatically” as possible. A capture time of 3318 ms is recommended by the code. 3320.0 ms is entered. A 4 to 2 merge time of 1659 ms is recommended by the code. 1660.0 ms is entered. A 2 to 1 merge time twice that of the 4 to 2 merge is recommended by the code. 3320.0 ms is entered. The code requests a “time fraction” FQMRG to be entered. This is a number (from 0.1 to 1.0) that determines the time at which data are recorded during the 4 to 2 merge. The number 0.5 is entered. This encodes the recording of data halfway through the 4 to 2 merge.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Illustration of bunch merge simulation in AGS

Following are examples of a 4 to 2 to 1 merge of proton bunches in AGS obtained by running the simulation code ags2mrg23. As the code runs, the user is prompted to enter several numbers. For these examples, RF harmonic number 12.0 is entered. An RF frequency of 4.453 913 448 73 MHz at this harmonic is entered. This gives proton G γ = 45.5. The merge simulation starts with a uniform distribution of unbunched protons. 2.0 eV-s is entered for the longitudinal emittance of the distribution. The code will recommend a voltage to capture the unbunched protons into 4 harmonic 12 buckets. A capture voltage of 0.0527 kV is recommended. 0.053 kV is entered. It is desirable to capture the protons as “adiabatically” as possible. Because of the long synchrotron period (372 ms), a capture time of 74,707 ms is recommended by the code. A shorter time, 10,000 ms, is entered to reduce the running time of the code. A 4 to 2 merge time of 37,353 ms is recommended by the code. 4000.0 ms is entered. A 2 to 1 merge time twice that of the 4 to 2 merge is recommended by the code. 8000.0 ms is entered. The code requests a “time fraction” FQMRG to be entered. This is a number (from 0.1 to 1.0) that determines the time at which data are recorded during the 4 to 2 merge. The number 0.5 is entered. This encodes the recording of data halfway through the 4 to 2 merge.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Sun's Series via Cyclotomic Multiple Zeta Values

We prove and generalize several recent conjectures of Z.-W. Sun surrounding binomial coefficients and harmonic numbers. We show that Sun's series and their analogs can be represented as cyclotomic multiple zeta values of levels N ϵ {4,8,12,16,24} , namely Goncharov's multiple polylogarithms evaluated at N -th roots of unity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

3D small-gain formula allowing strong focusing and harmonic lasing for a ring-based x-ray free electron laser oscillator

We present a detailed derivation of a formula for the small-gain calculation for an x-ray free electron laser oscillator (XFELO) based on a medium-energy (3–4 GeV) storage ring. We found harmonic lasing and strong focusing are essential for this beam energy range. Taking the small-signal low-gain formula developed by Kim and his colleagues, we modified it in such a way that the gain can be calculated without the “no focusing approximation,” and a strong focusing can be applied, as well as harmonic lasing. In this formula, the gain is represented as a product of two factors with one of them depending only on the harmonic number, undulator period, and gap. Using this factor, we show that it is favorable to use harmonic lasing to achieve hard x-ray FEL working in the small-signal low-gain regime with the medium-energy electron beam. Our formula also allows FEL optimization by varying the vertical gradient of the undulator, the vertical dispersion, and the horizontal and vertical focusing, independently. As an example, we applied this formula to study the feasibility of an XFELO option for the National Synchrotron Light Source II (NSLS-II) upgrade. Since a quite high peak current is required for the FEL, collective effects of beam dynamics in medium-energy synchrotrons significantly affect the electron beam parameters. We carried out a multiparameter optimization taking collective effects into account. Note, even though our example is for a ring-based XFELO at 3 to 4 GeV, the formula and, in particular, the approach developed here may be applied to other types of FELs. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2020

Medium Energy Hard X-ray FEL Gain

We develop a tool for calculation to study the conditions for a hard x-ray FEL oscillator based on electron beam in the medium enegy range from 3 to 4.5GeV. We show that the approach developed by K.J. Kim, et.al. for the small signal low gain formula can be modified so that the gain can be derived wihout taking the “no focusing approximation” adopted in the approach, so that a strong focusing can be applied. We also derive the formula to allow for gain calculation of harmonic lasing. The gain in this formula can be cast in the form of a product of two factors with one of them only depends on harmonic number, undulator period and gap. Thus this factor can be used to show that it is favorable to use harmonic lasing to achieve hard x-ray FEL working in the medium energy range and in the small signal low gain regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First direct observations of gear-changing in a collider

In this work we perform the first demonstration of gear-changing in a real world collider. Gear-changing refers to a collision scheme where each ring of a collider stores a different harmonic number of bunches. These bunches are kept synchronized using different velocities. Such a system has been theorized, but has now been demonstrated using the Double ElectroStatic Ion Ring ExpEriment (DESIREE) in Stockholm Sweden. The experiment was able to demonstrate a gear-changing system, with both four on three bunches and five on four bunches. We determined a measurable parameter that shows a gear-changing system out to $37500$ turns of the slow beam (1 second). We were also able to detect a measurable longitudinal beam-beam interaction. We developed insights into how to control this type of system, opening up new possibilities, and showing that DESIREE can be used for gear-changing research.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Ripple-induced neoclassical toroidal viscous torque in Augmented-First Plasma operation phase in ITER

Abstract A systematic calculation is performed on the ripple-induced neoclassical toroidal viscous (NTV) torque for new ITER scenarios designed for the Augmented-First Plasma (A-FP) operation phase with the full tungsten wall, where the plasma-wall gap is varied in view of mitigating the impact of tungsten wall-plasma interactions. The torque calculation includes drift kinetic response of the plasma thermal and energetic particles to the n = 18 ( n is the toroidal harmonic number) ripple field. For the plasma scenario with ~45 cm plasma-wall gap at the outboard mid-plane and considering the corrected ripple level of 0.17% by the ferritic steel inserts, the computed net NTV torque acting on the plasma column is in the sub-Nm level. However, with decreasing the plasma-wall gap, the computed net NTV torque can reach a level comparable to that produced by the neutral-beam momentum injection in ITER. Ripple correction by ferritic inserts reduces the net torque by a factor of 3.3 for all the three A-FP scenarios considered. The n ω d = l ω b (with ω d and ω b being the toroidal precession and bounce frequencies of trapped particles, respectively, and l an integer number) type of resonance-enhancement of the NTV torque, due to thermal particles, is found to be weak in ITER despite high- n of 18. The same also holds for the ITER 10 MA steady state scenario from the D-T operation phase, where the aforementioned resonance associated with fusion-born alphas is also included. The ripple-induced NTV torque is well below that produced by the resonant magnetic perturbation applied for controlling the type-I edge-localized mode in ITER.

ITER↗
DOE OSTI2023

Correlations between flow and transverse momentum in Xe + Xe and Pb + Pb collisions at the LHC with the ATLAS detector: A probe of the heavy-ion initial state and nuclear deformation

The correlations between flow harmonics v n for n = 2, 3, and 4 and mean transverse momentum [p T ] in 129 Xe + 129 Xe and 208 Pb + 208 Pb collisions at $\sqrt{s}$ = 5.44 and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are potentially sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from nonflow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and an event-activity selection based on particle production at very forward rapidity. The v n -[p T ] correlations show strong dependencies on centrality, harmonic number n, p T , and pseudorapidity range. Current models qualitatively describe the overall centralityand system-dependent trends but fail to quantitatively reproduce all features of the data. In central collisions, where models generally show good agreement, the v 2 -[p T ] correlations are sensitive to the triaxiality of the quadruple deformation. Comparison of the model with the Pb + Pb and Xe + Xe data confirms that the 129 Xe nucleus is a highly deformed triaxial ellipsoid that has neither a prolate nor oblate shape. This provides strong evidence for a triaxial deformation of the 129 Xe nucleus from high-energy heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Further Measurements of Beam-Beam Interactions in a Gear-Changing System in DESIREE

In this work we detail experiments performed on a gear-changing system using the Double ElectroStatic Ion Ring ExpEriment (DESIREE). A gear-changing system is one where there are different harmonic numbers in each ring. This experiment used carbon and nitrogen beams in a 4 on 3 gear-changing arrangement, with the last bunch of each left off. The bunch length can be measured and synchrotron motion detected. We performed this measurement on three different values of carbon current, and present the differences in the bunch length frequency spectrum here, which correspond to twice the synchrotron frequencies.

Nissen, E. A.↗
DOE OSTI2026

New measurement techniques for gear-changing research using DESIREE

In this work we cover some of the newer techniques developed to measure the effects of a gear changing system maintained in DESIREE at Stockholm University. Gear-changing is a collider synchronization method where two rings with different harmonic numbers in them maintain collisions through different velocities, pathlengths or a combination of the two. This system has been demonstrated using the low energy ion collider DESIREE at Stockholm university. We have not only continued our previous methods of studying the beam using a repeating pattern technique where one bucket in each ring is intentionally left empty, but we now also use recently installed pickups outside of the merger region to study the beams separately while they collide.

Accelerator Physics↗
DOE OSTI2024

Dimensionally Aligned Signal Projection Algorithms Library

Dimensionally aligned signal projection (DASP) algorithms are used to analyze fast Fourier transforms (FFTs) and generate visualizations that help focus on the harmonics for specific signals. At a high level, these algorithms extract the FFT segments around each harmonic frequency center, and then align them in equally sized arrays ordered by increasing distance from the base frequency. This allows for a focused view of the harmonic frequencies, which, among other use cases, can enable machine learning algorithms to more easily identify salient patterns. This work seeks to provide an effective open-source implementation of the DASP algorithms proposed by Vann et al. (2018) as well as functionality to help explore and test how these algorithms work with an interactive dashboard and signal-generation tool. The DASP library is implemented in Python and contains four types of algorithms for implementing these feature engineering techniques: fixed harmonically aligned signal projection (HASP), decimating HASP, interpolating HASP, and frequency aligned signal projection (FASP). Each algorithm returns a numerical array, which can be visualized as an image. The HASP algorithms are variations of the algorithms originally presented by Vann et al. (2018). For consistency, FASP, which is the terminology used for the short-time Fourier transform (STFT), has been implemented as part of the library to provide a similar interface to the STFT of the raw signal. Additionally, the library contains an algorithm to generate artificial signals with basic customizations such as the base frequency, sample rate, duration, number of harmonics, noise, and number of signals. Finally, the library provides multiple interactive visualizations, each of which is implemented using IPyWidgets and works in a Jupyter environment. A dashboard-style visualization is provided, which contains some common signal-processing visual components (signal, FFT, spectogram) updating in unison with the HASP functions (see Figure 1 below). Separate from the dashboard, an independent visualization is provided for each of the DASP algorithms as well as for the artifical signal generator. These visualizations are included in the library to aid in developing an intuitive understanding how the algorithms are affected by different input signals and parameter selections.

harmonics↗
DOE OSTI2024

Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. Further, the method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
DOE OSTI2021

Nonuniform Absorption and Scattered Light in Direct-Drive Implosions Driven by Polarization Smoothing

Laser-direct-drive symmetric implosions on OMEGA illuminate a target with 60 laser beams and are designed to produce spherical implosions. Each beam is smoothed using orthogonal polarizations obtained by passing through distributed polarization rotators (DPR’s). Observations of the scattered light from OMEGA implosions do not show the expected symmetry and have much larger variation than standard predictions. For the first time, we have quantified the scattered-light nonuniformity from individual beams and identified the DPR’s as the source of the enhanced nonuniformity. An instrument was invented that isolated and measured the variation in the intensity and polarization of the light scattered from each OMEGA beam. The asymmetric intensity and polarization measurements are explained when the on-target offsets between the two orthogonal polarizations produced by the DPR’s are modeled using a 3-D cross-beam energy transfer (CBET) code that tracks the polarizations of each beam. The time-integrated nonuniformity in laser absorption and scattered light due to CBET and the DPR polarization offsets during high-performance OMEGA implosions is predicted to be significant and dominated by low spherical harmonic mode numbers. Further, the nonuniformity is predicted to be greatly reduced by replacing the DPR’s with new optics that create smaller offsets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Bounds on spectral gaps of Hyperbolic spin surfaces

We describe a method for constraining Laplacian and Dirac spectra of two dimensional compact orientable hyperbolic spin manifolds and orbifolds. The key ingredient is an infinite family of identities satisfied by the spectra. These spectral identities follow from the consistency between 1) the spectral decomposition of functions on the spin bundle into irreducible representations of SL(2,R) and 2) associativity of pointwise multiplication of functions. Applying semidefinite programming methods to our identities produces rigorous upper bounds on the Laplacian spectral gap as well as on the Dirac spectral gap conditioned on the former. In several examples, our bounds are nearly sharp; a numerical algorithm based on the Selberg trace formula shows that the [0;3,3,5] orbifold, a particular surface with signature [1;3], and the Bolza surface nearly saturate the bounds at genus 0, 1 and 2 respectively. Under additional assumptions on the number of harmonic spinors carried by the spin-surface, we obtain more restrictive bounds on the Laplacian spectral gap. In particular, these bounds apply to hyperelliptic surfaces. We also determine the set of Laplacian spectral gaps attained by all compact orientable two-dimensional hyperbolic spin orbifolds. We show that this set is upper bounded by 12.13798; this bound is nearly saturated by the [0;3,3,5] orbifold, whose first non-zero Laplacian eigenvalue is λ^(0)_1 ≈ 12.13623.

Spectral theory↗
DOE OSTI2020

Limiting the number of quantisation levels of a harmonic lens as a method for improving the quality of the generated image

The characteristics (modulation transfer functions) of harmonic lenses made by direct laser patterning of resist are experimentally studied. The influence of technological errors in the manufacture of lenses on the formation of the point spread function is modelled for harmonic lenses with 32 and 256 quantisation levels. Direct laser writing is used to produce lenses with these numbers of levels. It is experimentally shown that the modulation transfer function of the harmonic lens with 32 quantisation levels is significantly greater at high spatial frequencies (50 mm{sup −1}). (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗