Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic oscillator”

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.

249 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phonon-enhanced nonlinearities in hexagonal boron nitride

Abstract Polar crystals can be driven into collective oscillations by optical fields tuned to precise resonance frequencies. As the amplitude of the excited phonon modes increases, novel processes scaling non-linearly with the applied fields begin to contribute to the dynamics of the atomic system. Here we show two such optical nonlinearities that are induced and enhanced by the strong phonon resonance in the van der Waals crystal hexagonal boron nitride (hBN). We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies beyond the third.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline-based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive-sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

current measurement↗
DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis: Preprint

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline- based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive- sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

grid following inverter↗
DOE OSTI2022

Dark Energy Survey Year 3 results: A 2.7% measurement of baryon acoustic oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of baryon acoustic oscillations (BAO) in an optimized sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg 2 with 0.6 < z photo < 1.1 and a typical redshift uncertainty of 0.03 ( 1 + z ) . The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift, and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive preunblinding protocol. We used two different statistics; angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance D M at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain D M ( z eff = 0.835 ) / r d = 18.92 ± 0.51 , which is consistent with, but smaller than, the Planck prediction assuming flat Λ CDM , at the level of 2.3 σ . The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at z > 0.75 . When combined with DES 3 x 2 pt + SNIa , they lead to improvements in H 0 and Ω m constraints by ∼ 20 % .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tidally Excited Modes and δ Scuti Pulsations in the Eclipsing Triple Star IM Persei

IM Per is a triple star system whose eclipsing pair masses and radii are accurate to within 1%. The Transiting Exoplanet Survey Satellite (TESS) light curve of the program target exhibits partial eclipses and multiple oscillations with millimagnitude-level amplitudes. It is found that the oscillations affect eclipse timing measurements. Binary modeling of the high-quality TESS data indicates that the eclipsing components of the triple system are twin dwarfs with parameters of M {sub 2}/M {sub 1} = 0.995, R {sub 2}/R {sub 1} = 0.901, and Δ(T{sub eff,1}–T{sub eff,2}) = 12 K in an eccentric (e = 0.049), detached configuration. The third light of l {sub 3} = 0.054 may mostly come from a G-type tertiary companion. Our predicted parallax of 1.52 ± 0.09 mas is concurrent with the Gaia measurement of 1.52 ± 0.05 mas. Multifrequency analysis of the outside-eclipse residuals reveals 22 significant pulsation signals: 4 in the gravity-mode region (0.03−2.22 day{sup −1}) and 18 in the pressure-mode region (9.19−25.12 day{sup −1}). Of the low frequencies, f {sub 11} and f {sub 14} are orbital harmonics that can be identified as tidally excited modes. The pulsation periods and constants for the high frequencies, and the position in the Cepheid instability strip demonstrate that the pulsating component of IM Per is a δ Sct variable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ultrafast Nanoimaging of Spin-Mediated Shear Waves in an Acoustic Cavity

Strong spin–lattice coupling in van der Waals (vdW) magnets shows potential for innovative magneto-mechanical applications. Here, nanoscale and picosecond imaging by ultrafast electron microscopy reveal heterogeneous spin-mediated coherent acoustic phonon dynamics in a thin-film cavity of the vdW antiferromagnet FePS3. The harmonics of the interlayer shear acoustic modes are observed, in which the even and odd harmonics exhibit distinct nanoscopic dynamics. Corroborated by acoustic wave simulation, the role of defects in forming even harmonics is elucidated. Above the Néel temperature (T N ), the interlayer shear acoustic harmonics are suppressed, while the in-plane traveling wave is predominantly excited. The dominant acoustic dynamics shifts from the out-of-plane shear to the in-plane traveling wave across T N , demonstrating that magnetic properties can influence phonon scattering pathways. In conclusion, the spatiotemporally resolved structural characterization provides valuable nanoscopic insights for interlayer-shear-mode-based acoustic cavities, opening up possibilities for magneto-mechanical applications of vdW magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗
DOE OSTI2021

TESS cycle 1 observations of roAp stars with 2-min cadence data

ABSTRACT We present the results of a systematic search for new rapidly oscillating Ap (roAp) stars using the 2-min cadence data collected by the Transiting Exoplanet Survey Satellite (TESS) during its Cycle 1 observations. We identify 12 new roAp stars. Amongst these stars we discover the roAp star with the longest pulsation period, another with the shortest rotation period, and six with multiperiodic variability. In addition to these new roAp stars, we present an analysis of 44 known roAp stars observed by TESS during Cycle 1, providing the first high-precision and homogeneous sample of a significant fraction of the known roAp stars. The TESS observations have shown that almost 60 per cent (33) of our sample of stars are multiperiodic, providing excellent cases to test models of roAp pulsations, and from which the most rewarding asteroseismic results can be gleaned. We report four cases of the occurrence of rotationally split frequency multiplets that imply different mode geometries for the same degree modes in the same star. This provides a conundrum in applying the oblique pulsator model to the roAp stars. Finally, we report the discovery of non-linear mode interactions in α Cir (TIC 402546736, HD 128898) around the harmonic of the principal mode – this is only the second case of such a phenomenon.

Holdsworth, D. L.↗
DOE OSTI2022

State-selective probing of CO 2 autoionizing inner valence Rydberg states with attosecond extreme ultraviolet four-wave-mixing spectroscopy

Nonlinear spectroscopies can disentangle spectra that are congested due to inhomogeneous broadening. Here, in conjunction with theoretical calculations, attosecond extreme ultraviolet (XUV) four-wave-mixing (FWM) spectroscopy is utilized to probe the dynamics of autoionizing inner valence excited Rydberg states of the polyatomic molecule, CO 2 . This tabletop nonlinear technique employs a short attosecond XUV pulse train and two noncollinear, few-cycle near-infrared pulses to generate background-free XUV wave-mixing signals. FWM emission is observed from the n=5-7 states of the Henning sharp ndσ g Rydberg series that converges to the ionic $\widetilde{B}$ 2 Σ$^{+}_{u}$ state. However, these transient emission signals decay with lifetimes of 33 ± 6, 53 ± 2, and 94 ± 2 fs, respectively, which calculations show are consistent with the lifetimes of the short-lived n=6-8 members of the nsσ g character Henning diffuse Rydberg series. The oscillator strengths of transitions between states involved in all possible resonant FWM processes are calculated, verifying that the nonlinear spectra are dominated by pathways described by an initial excitation to the diffuse nsσ g Rydberg series and emission from the sharp ndσ g Rydberg series. The results substantiate not only that attosecond XUV FWM spectroscopy produces rigorous and meaningful measurements of ultrafast dynamics in polyatomic systems, but also that nonlinear spectroscopic techniques are versatile tools to selectively probe dynamics that are otherwise difficult to access.

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

The r -modes of slowly rotating, stratified neutron stars

ABSTRACT The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental l = |m| solutions, where l and m are the indices of the spherical harmonic $Y_l^m$ that describe the mode’s angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the l ≥ |m| perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star’s ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental l = |m| perturbations change character and become generic inertial modes. Restricting the analysis to l = |m| perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. In order to see which r-modes persist and retain their character in realistic neutron stars, these calculations will need to be brought into full general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Electro-optic characterization of synthesized infrared-visible light fields

The measurement and control of light field oscillations enable the study of ultrafast phenomena on sub-cycle time scales. Electro-optic sampling (EOS) is a powerful field characterization approach, in terms of both sensitivity and dynamic range, but it has not reached beyond infrared frequencies. Here, we show the synthesis of a sub-cycle infrared-visible pulse and subsequent complete electric field characterization using EOS. The sampled bandwidth spans from 700 nm to 2700 nm (428 to 110 THz). Tailored electric-field waveforms are generated with a two-channel field synthesizer in the infrared-visible range, with a full-width at half-maximum duration as short as 3.8 fs at a central wavelength of 1.7 µm (176 THz). EOS detection of the complete bandwidth of these waveforms extends it into the visible spectral range. To demonstrate the power of our approach, we use the sub-cycle transients to inject carriers in a thin quartz sample for nonlinear photoconductive field sampling with sub-femtosecond resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Experimental Investigation of Interactions Between Two Closely Spaced Azimuthal Modes in a Multinozzle Can Combustor

Thermoacoustic instabilities in annular or circular combustors are often coupled with azimuthal modes. These modes can be characterized by a spin ratio (SR), which quantifies the dominant mode between two counter-rotating waves, and a phase difference between them, which is directly related to the orientation of the antinodal line. This study investigates the instability amplitudes, SR, and phase difference of two closely spaced (3% of their mean frequency), yet distinct azimuthal modes; one is the first azimuthal (1A) mode, and the other is a combination of first azimuthal and first longitudinal (1A1L) mode. Each mode itself consists of two peaks that are spaced even more closely in frequency (0.8%). Furthermore, distinct harmonics at 2× and 3× of these frequencies, presumably associated with nonlinearities, are also evident in the spectra. Each mode is bandpass filtered in spectrum to analyze them separately. For the 1A mode, the SR and phase difference exhibit a variety of behaviors—including quasi-periodic standing waves, spinning waves, and intermittency—depending on operating conditions such as thermal power and azimuthal fuel staging. Similar trends are observed for the 1A1 L mode. Moreover, there is clear coupling between the 1A and 1A1 L modes, as their SRs are almost synchronized during the quasi-periodic standing wave. This synchronization is observed in phase differences as well, but not in the instability amplitude. For spinning dominant wave conditions, the SRs of each mode have similar average values, but they fluctuate in a seemingly random fashion. For the phase difference, both average and fluctuation are not correlated. In contrast, the instability amplitudes are strongly correlated, with modulation of the 1A mode leading to that of the 1A1 L mode. Furthermore, these results clearly indicate that complex coupling occurs across closely spaced frequencies under instability conditions, coupling that must be understood in order to capture limit cycle dynamics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗
DOE OSTI2023

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗