Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic instability”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

MARS-Q modeling of low- n resistive kink-peeling modes and edge harmonic oscillations in DIII-D

Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low-n (n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low-n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

EHO↗
DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Absence of High-Pressure Ground-State Reentrant Ferroelectricity in PbTiO 3

We study ferroelectricity in the classic perovskite ferroelectric PbTiO 3 to high pressures with density functional theory (DFT) and experimental diamond-anvil techniques. We use second harmonic generation spectroscopy to detect lack of inversion symmetry. Consistent with early understanding and experiments, we find that ferroelectricity disappears at moderate pressures. However, DFT computations show that the disappearance arises from the overtaking of zone boundary instabilities, and not from the squeezing out of the off-centering ferroelectric displacements with pressure, as previously thought. Moreover, at high pressures the distorted perovskite phases are metastable with respect to a new dense centrosymmetric postperovskite phase with P⁢2 1 /m symmetry and 8-coordinated Ti, which becomes stable at about 70 GPa.

36 MATERIALS SCIENCE↗
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 OSTI2024

Ion cyclotron emission in Maxwellian plasmas

Destabilization of ion cyclotron waves—waves with frequencies close to ion cyclotron harmonics—in inhomogeneous plasmas with Maxwell velocity distribution is considered. A new mechanism of destabilization of these waves is found, whereas the known interaction between cyclotron waves and drift waves is shown to be hardly able to lead to instabilities in realistic tokamak plasmas. Our finding is that the resonance wave–particle interaction in the presence of temperature gradient can change the diamagnetic drift frequency in such a way that the destabilizing influence of the diamagnetic drift exceeds Landau damping. This occurs when the energy of resonance particles is sufficiently high, which is always the case due to the infinite tail of the Maxwell distribution. Particles with higher energies provide a larger ratio of drive to damping, but the number of resonance particles and, thus, instability growth rates are exponentially small when particle energy is very high. Therefore, only moderately suprathermal particles can lead to observable instabilities. A condition for instabilities driven by these particles is obtained. Destabilization of electrostatic cyclotron waves and ordinary cyclotron waves is studied.

Kolesnichenko, Ya. I. (ORCID:0000000323219681)↗
DOE OSTI2024

Extreme-value statistics in nonlinear optics

We show that, although nonlinear optics may give rise to a vast multitude of statistics, all these statistics converge, in their extreme-value limit, to one of a few universal extreme-value statistics. Specifically, in the class of polynomial nonlinearities, such as those found in the Kerr effect, weak-field harmonic generation, and multiphoton ionization, the statistics of the nonlinear-optical output converges, in the extreme-value limit, to the exponentially tailed, Gumbel distribution. Exponentially growing nonlinear signals, on the other hand, such as those induced by parametric instabilities and stimulated scattering, are shown to reach their extreme-value limits in the class of the Fréchet statistics, giving rise to extreme-value distributions (EVDs) with heavy, manifestly nonexponential tails, thus favoring extreme-event outcomes and rogue-wave buildup.

Zheltikov, Aleksei M. (ORCID:0000000291380576)↗
DOE OSTI2021

Particle resonances in toroidal fusion devices

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through avalanche. The existence of a mode particle resonance depends on properties of the equilibrium, particle trajectories, and perturbation mode harmonic content. Resonance location is a function of particle energy and equilibrium field line helicity. Different methods for finding resonance location and energy dependence are developed. The properties of resonances are discussed using examples from magnetic fusion devices. Here, we show that if mode resonances exist at low particle energy, they very likely also exist at high energy, thus modifying high energy beam particles and fusion products. It is possible for a resonance to appear due to mode induced orbit helicity modification when it is forbidden at low mode amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unlocking New Regimes in Fractional Quantum Hall Effect with Quaternions

We demonstrate that formulating the composite-fermion theory of the fractional quantum Hall (FQH) effect in terms of quaternions greatly expands its reach and opens the door into many interesting issues that were previously not amenable to quantitative theoretical investigation. As an illustration, we explore the possibility of a nematic or a charge-density wave instability of the composite-fermion Fermi sea at half-filled Landau level and of the nearby FQH states by looking for a gap closing instability of the neutral magneto-roton excitation. As a result, our quaternion formulation of the FQH effect has been inspired by mathematical developments in the theoretical analyses of gravitational wave modes and cosmic microwave background radiation, where an important role is played by spin-weighted spherical harmonics that are nothing but monopole harmonics appearing in the spherical geometry for the FQH effect.

Composite fermions↗
DOE OSTI2023

Species mix & distribution function dependence of energetic-ion driven magnetospheric instabilities on DIII-D (Final Report)

Energetic ions drive a rich variety of instabilities in the complex multi-species plasma of Earth’s magnetosphere. This experiment exploited the highly flexible, high power particle beam injection capability and sophisticated diagnostic suite of the DIII-D tokamak to advance understanding of two particular instabilities of significant interest to the magnetospheric research community that have closely-related counterparts in beam-heated tokamak plasmas: electromagnetic ion cyclotron waves (EMIC) [1], observed at frequencies somewhat below the ion cyclotron frequencies of bulk ion species, and equatorial noise [2], observed at energetic-ion cyclotron harmonics. In the magnetic fusion energy domain, such waves are called ICE (Ion Cyclotron Emission). Recent substantial advances in diagnostic capability in the ion cyclotron range of frequencies in DIII-D, including measurement of wave polarization, spa)al structure and amplitude created an exciting opportunity to challenge and drive advances in theoretical understanding of these phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Light-induced Weyl semiconductor-to-metal transition mediated by Peierls instability

Elemental tellurium is a strongly spin-orbit coupled Peierls-distorted semiconductor whose band structure features topologically protected Weyl nodes. Here, using time-dependent density functional theory calculations, we show that impulsive optical excitation can be used to transiently control the amplitude of the Peierls distortion, realizing a mechanism to switch tellurium between three states: Weyl semiconductor, Weyl metal and non-Weyl metal. Further, we present experimental evidence of this inverse-Peierls distortion using time-resolved optical second harmonic generation measurements. These results provide a pathway to multifunctional ultrafast Weyl devices and introduce Peierls systems as viable hosts of light-induced topological transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Higher Harmonics in Multipactor Induced Plasma Ionization Breakdown near a Dielectric Surface

In this Letter, the novel physics of higher harmonic (HH) generation of the normal electric field near a dielectric surface is reported in multipactor induced plasma ionization breakdown, as determined by kinetic particle-in-cell simulations. The observed HH frequency is around ten times the fundamental rf driving frequency, but lower than the electron plasma frequency. A theory is constructed which indicates that stream plasma interaction-induced instability is the mechanism of HH generation in the collisional regime. The HH frequency and its corresponding growth rate of the HH oscillation amplitude from the theory are in good agreement with kinetic particle-in-cell simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Ion-driven destabilization of a toroidal electron plasma—A 3D3V PIC simulation

Ion-driven destabilization of a toroidal electron plasma in a small aspect ratio axisymmetric toroidal device is reported for A r + ions of different initial density values using a high fidelity 3D3V PIC solver. Stability of a recently discovered quiescent quasi-steady state (QQS) of a toroidal electron plasma obtained from “seed” solution as a result of entropy extremization at zero inertia is addressed in the presence of a small ion population. An initial value (f0) of the ion fraction (f = n i / n e) and the corresponding secondary electrons are “preloaded” into the system after the electron plasma attains a QQS state. This procedure is regarded as a proxy to the conventional production of ions in the experimental devices via impact ionization. The resulting electron plasma exhibits destabilized “center of charge motion” (m = 1) along with higher order harmonics with dominant power in the second harmonic. Gradual loss of ions (and also electrons) is observed resulting in time varying f values. Beyond a certain value of f0 (≥ 0.005), growth in wall probe current is observed, which saturates at later simulation time due to the loss of particles. In conclusion, trajectories of ion particles indicate ion trapping in the potential well, which is qualitatively similar to the ion resonance instability in pure electron plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accelerated Wind-Turbine Wake Recovery Through Actuation of the Tip-Vortex Instability

Advances in wind-plant control have often focused on more effectively balancing power between neighboring turbines. Wake steering is one such method that provides control-based improvements in a quasi-static way, but this does little to fundamentally change the wake recovery process, and thus, it has limited potential. This study investigates use of another control paradigm known as dynamic wake control (DWC) to excite the mutual inductance instability between adjacent tip-vortex structures, thereby accelerating the breakdown of the structures. The current work carries this approach beyond the hypothetical by applying the excitation via turbine control vectors that already exist on all modern wind turbines: blade pitch and rotor speed control. The investigation leverages a free-vortex wake method (FVWM) that allows a thorough exploration of relevant frequencies and amplitudes of harmonic forcing for each control vector (as well as the phase difference between the vectors for a tandem configuration) while still capturing the essential tip-vortex dynamics. The FVWM output feeds into a Fourier stability analysis working to pinpoint candidate DWC strategies suggesting fastest wake recovery. Near-wake length reductions of >80% are demonstrated, although without considering inflow turbulence. Analysis is provided to interpret these predictions considering the presence of turbulence in a real atmospheric inflow.

17 WIND ENERGY↗
DOE OSTI2024

Three-dimensional reconstruction of x-ray emission volumes in magnetized liner inertial fusion from sparse projection data using a learned basis

The ability to visualize x-ray and neutron emission from fusion plasmas in 3D is critical to understand the origin of the complex shapes of the plasmas in experiments. Unfortunately, this remains challenging in experiments that study a fusion concept known as Magnetized Liner Inertial Fusion (MagLIF) due to a small number of available diagnostic views. Here, we present a basis function-expansion approach to reconstruct MagLIF stagnation plasmas from a sparse set of x-ray emission images. A set of natural basis functions is “learned” from training volumes containing quasi-helical structures whose projections are qualitatively similar to those observed in experimental images. Tests on several known volumes demonstrate that the learned basis outperforms both a cylindrical harmonic basis and a simple voxel basis with additional regularization, according to several metrics. Two-view reconstructions with the learned basis can estimate emission volumes to within 11% and those with three views recover morphology to a high degree of accuracy. The technique is applied to experimental data, producing the first 3D reconstruction of a MagLIF stagnation column from multiple views, providing additional indications of liner instabilities imprinting onto the emitting plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗