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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Power Solitons in Inverter-based Electric Energy Systems: Observation, Analysis and Implications

The discovery of an unfamiliar class of nonlinear energy waves in electric energy systems with dominating coupled inverter-based energy resources has been reported in a recent research project. Whenever there are slight but persistent differences among the reference signals used by inverters in pulse-width modulation, i.e., a weak group modulation instability condition, the unfamiliar energy waves appear, which contain not only the familiar frequency components present in electric power systems but also self-organized components in the ultra-high frequency region of the energy spectrum. These components are found to be exploding solitary waves with concentrated energy masses, hence we call them power solitons. This paper explains the work of the discovery reported, and the implications of the presence of power solitons for grid integration of inverter-based resources and power system protection.

14 SOLAR ENERGY↗
DOE OSTI2022

Instabilities of a vortex-ring-bright soliton in trapped binary three-dimensional Bose-Einstein condensates

Instabilities of vortex-ring-bright coherent structures in harmonically trapped two-component three-dimensional Bose-Einstein condensates are studied numerically within the coupled Gross-Pitaevskii equations and interpreted analytically. Interestingly, the filled vortex core with a sufficiently large amount of the bright component is observed to reduce the parametric interval of stability of the vortex ring. Here, we have identified the mechanisms of several linear instabilities and one nonlinear parametric instability in this connection. Two of the linear instabilities are qualitatively different from ones reported earlier and are associated with azimuthal modes of m = 0 and m = 1, i.e., deviations of the vortex from the stationary ring shape. Our nonlinear parametric resonance instability occurs between the m = 0 and m = 2 modes and signals the exchange of energy between them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Multi-Port Autonomous Reconfigurable Solar Power Plant

As the penetration level of power electronics increases and remote photovoltaic (PV) generation is integrated into the alternating current (ac) grid, the short-circuit ratio (SCR) at the point of interconnection of a hybrid PV–energy storage system (ESS) plant may be low. Additionally, the inertia of the alternating current (ac) grid may be low. The low SCRs and inertias can lead to reliability challenges in the power grid. These operating conditions require additional reinforcements, such as synchronous condensers, static var compensators, static synchronous compensators, and high-voltage direct current (HVdc) links/grids. HVdc links or grids may also provide the additional capability of connecting the plant to asynchronous grids and/or connecting asynchronous grids, among others. This scenario leads to discrete development of solar inverters, energy storage inverters, HVdc converters, and several transformers. Some of the problems associated with this discrete development and inverter-based generation include increased cost, lower reliability, and reduced efficiency associated with duplication of power electronics (PEs); competing controls of several individual discrete inverter-based generators due to the presence of multiple PEs, which leads to derating of the system; and transient stability problems arising from inverter-based generation, such as voltage/frequency events leading inverter shutdowns, voltage instability and control interactions in the formed weak grid, and harmonics caused by resonances of multiple inverters.

14 SOLAR ENERGY↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

O-Mode Instability for Non-Thermal Continuum Radiation Generated Near the Plasmapause

Analysis of the ordinary mode (O-mode) instability is performed to comprehend the nonthermal continuum (NTC) radiation near the plasmapause, taking into account the relativistic wave-electron resonance effect. The energy source is the anisotropy in the velocity of the minority suprathermal electron population. Numerical solutions demonstrate that the O-mode can be unstable with multiple narrow frequency bands located close to harmonics of the electron cyclotron frequency above the local electron plasma frequency. These waves have narrow beaming angle bands of nearly 90° relative to the ambient magnetic field. Our findings indicate that NTC radiation generated by this wave-electron resonance instability near the plasmapause can propagate nearer to the magnetic equator with multiple harmonics, which is in agreement with a recent statistical study using Van Allen Probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Excitation of toroidally localized harmonics of global Alfvén eigenmodes

Abstract The spherical tokamak NSTX and the upgraded, higher toroidal field, NSTX-U, routinely create plasmas with large un-thermalized populations of super-Alfvénic fast ions. This population of non-thermal energetic ions excites a broad spectrum of Alfvénic waves, from the lower frequency (≈30 kHz–≈200 kHz) toroidal Alfvén eigenmodes (Cheng and Chance 1986 Phys. Fluids 29 3695), to the higher frequency (400 kHz–3 MHz) global Alfvén eigenmodes (GAEs) and compressional Alfvén eigenmodes (Goedbloed 1975 Phys. Fluids 18 1258). In this paper we present evidence that the GAE non-linearly excite modes, presumably GAE, at frequencies consistent with non-linear or 3-wave coupling. The observation of the excitation of 2 nd harmonic GAE through the intrinsic non-linearity of Global Alfvén modes demonstrates that the non-linear terms can act as an exciter-antenna inside the plasma, broadcasting at harmonics of the mode frequency and with concomitantly shorter wavelengths. As with experiments using external antenna to excite otherwise weakly stable Alfvénic modes, this data can provide information about the stability of modes at harmonics of the GAE. It may also provide information on the nature of the non-linearities in the wave dispersion equation or a direct measure of the mode amplitude (Smith et al 2006 Phys. Plasmas 13 042504). We also report that the short wavelength 2 nd harmonic GAE can be strongly toroidally localized. The observation that the shorter wavelength ( n ≈ 20) 2 nd harmonic modes can be toroidally localized potentially has implications for the impact of fast-ion driven instabilities on fast-ion confinement in ITER and future fusion reactors (Gorelenkov et al 2014 Nucl. Fusion 54 125001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Booster cavity damper redesign for PIP-II

A new Higher Order Mode (HOM) damper was designed and is undergoing testing for the Booster accelerator cavity at Fermilab. In anticipation of the PIP-II upgrade, it was discovered that the higher intensity of PIP-II may cause beam instability due to an excited mode at 106 MHz. This unfortunately corresponds with the cavity’s 2nd order harmonic mode, which will sweep from 86-105.7 MHz. The new damper is a modification of an existing damper that was designed to reduce an existing static HOM at 83 MHz, with the new design intending to cover the 2nd order HOM as well. The existing damper uses an inductive coupling loop to extract RF energy from the cavity which then goes through a filter in order to reflect the fundamental frequency back into the cavity while passing HOMs to a dump load. The new damper intends to replace the filter portion of the system with a wider band variant while also changing the topology from a coaxial cable loop filter to a componentized PCB-based design. Primary design challenges include bandwidth coverage, impedance matching of the various modes, long term thermal and mechanical stability, radiation hardness, and high voltage handling. Initial designs achieved the desired damping but were found to quickly succumb to destructive arcing due to the voltages present. More finalized designs intend to address this problem through circuit design modifications as well as the use of hardier components.

Pieper, Dustin [Fermilab]↗
DOE OSTI2025

Booster Cavity Damper Redesign for PIP-II

A new Higher Order Mode (HOM) damper was designed and tested for the Booster accelerator cavity at Fermilab. In anticipation of the PIP-II upgrade, it was discovered that the higher beam intensity of PIP-II may cause beam instability due to an excited mode at 106 MHz. This unfortunately corresponds with the cavity s 2nd order harmonic mode, which will sweep from 86-105.7 MHz. The new damper is a modification of an existing damper that was designed to reduce an existing static HOM at 83 MHz, with the new design intending to cover the 2nd order HOM as well. The existing damper uses an inductive coupling loop to extract RF energy from the cavity which then goes through a filter in order to reflect the fundamental frequency back into the cavity while passing HOMs to a dump load. The new damper intends to replace the filter portion of the system with a wider band variant while also changing the topology from a coaxial cable loop filter to a componentized PCB-based design. Primary design challenges include bandwidth coverage, impedance matching of the various modes, long term thermal and mechanical stability, radiation hardness, and high voltage handling. Initial designs achieved the desired damping but were found to quickly succumb to destructive arcing due to the voltages present. More finalized designs intend to address this problem through circuit design modifications as well as the use of hardier components.

Pieper, Dustin [Fermilab]↗
DOE OSTI2025

Reliable Determination of Pulses and Pulse-Shape Instability in Ultrashort Laser Pulse Trains Using Polarization-Gating and Transient-Grating Frequency-Resolved Optical Gating Using the RANA Approach

Devices that measure the presence of instability in the pulse shapes in trains of ultrashort laser pulses do not exist, so this task necessarily falls to pulse-measurement devices, like Frequency-Resolved Optical Gating (FROG) and its variations, which have proven to be a highly reliable class of techniques for measuring stable trains of ultrashort laser pulses. Fortunately, multi-shot versions of FROG have also been shown to sensitively distinguish trains of stable from those of unstable pulse shapes by displaying readily visible systematic discrepancies between the measured and retrieved traces in the presence of unstable pulse trains. However, the effects of pulse-shape instability and algorithm stagnation can be indistinguishable, so a never-stagnating algorithm—even when instability is present—is required and is generally important. In previous work, we demonstrated that our recently introduced Retrieved-Amplitude N-grid Algorithmic (RANA) approach produces highly reliable (100%) pulse-retrieval in the second-harmonic-generation (SHG) version of FROG for thousands of sample trains of pulses with stable pulse shapes. Further, it does so even for trains of unstable pulse shapes and thus both reliably distinguishes between the two cases and provides a rough measure of the degree of instability as well as a reasonable estimate of most typical pulse parameters. Here, we perform the analogous study for the polarization-gating (PG) and transient-grating (TG) versions of FROG, which are often used for higher-energy pulse trains. We conclude that PG and TG FROG, coupled with the RANA approach, also provide reliable indicators of pulse-shape instability. In addition, for PG and TG FROG, the RANA approach provides an even better estimate of a typical pulse in an unstable pulse train than SHG FROG does, even in cases of significant pulse-shape instability.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

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