Abstract:The present analysis makes use of the Vlasov-Maxwell equations to develop a fully kinetic description of the electrostatic, electron-ion two-stream instability driven by the directed axial motion of a high-intensity ion beam propagating in the z direction with average axial momentum γ_{b}m_{b}β_{b}c through a stationary population of background electrons. The ion beam has characteristic radius r_{b} and is treated as continuous in the z direction, and the applied transverse focusing force on the beam ions is modeled by F_{foc}^{b}=-γ_{b}m_{b}ω_{βb}^{0^{2}}x_{⊥} in the smooth-focusing approximation. Here, ω_{βb}^{0}=const is the effective betatron frequency associated with the applied focusing field, x_{⊥} is the transverse displacement from the beam axis, (γ_{b}-1)m_{b}c^{2} is the ion kinetic energy, and V_{b}=β_{b}c is the average axial velocity, where γ_{b}=(1-β_{b}^{2})^{-1/2}. Furthermore, the ion motion in the beam frame is assumed to be nonrelativistic, and the electron motion in the laboratory frame is assumed to be nonrelativistic. The ion charge and number density are denoted by +Z_{b}e and n_{b}, and the electron charge and number density by -e and n_{e}. For Z_{b}n_{b}>n_{e}, the electrons are electrostatically confined in the transverse direction by the space-charge potential φ produced by the excess ion charge. The equilibrium and stability analysis retains the effects of finite radial geometry transverse to the beam propagation direction, including the presence of a perfectly conducting cylindrical wall located at radius r=r_{w}. In addition, the analysis assumes perturbations with long axial wavelength, k_{z}^{2}r_{b}^{2}≪1, and sufficiently high frequency that |ω/k_{z}|≫v_{Tez} and |ω/k_{z}-V_{b}|≫v_{Tbz}, where v_{Tez} and v_{Tbz} are the characteristic axial thermal speeds of the background electrons and beam ions. In this regime, Landau damping (in axial velocity space v_{z}) by resonant ions and electrons is negligibly small. We introduce the ion plasma frequency squared defined by ω[over ^]_{pb}^{2}=4πn[over ^]_{b}Z_{b}^{2}e^{2}/γ_{b}m_{b}, and the fractional charge neutralization defined by f=n[over ^]_{e}/Z_{b}n[over ^]_{b}, where n[over ^]_{b} and n[over ^]_{e} are the characteristic ion and electron densities. The equilibrium and stability analysis is carried out for arbitrary normalized beam intensity ω[over ^]_{pb}^{2}/ω_{βb}^{0^{2}}, and arbitrary fractional charge neutralization f, consistent with radial confinement of the beam particles. For the moderately high beam intensities envisioned in the proton linacs and storage rings for the Accelerator for Production of Tritium and the Spallation Neutron Source, the normalized beam intensity is typically ω[over ^]_{pb}^{2}/ω_{βb}^{0^{2}}≲ 0.1. For heavy ion fusion applications, however, the transverse beam emittance is very small, and the space-charge-dominated beam intensity is much larger, with ω[over ^]_{pb}^{2}/ω_{βb}^{0^{2}}≲ 2γ_{b}^{2}. The stability analysis shows that the instability growth rate Imω increases with increasing normalized beam intensity ω[over ^]_{pb}^{2}/ω_{βb}^{0^{2}} and increasing fractional charge neutralization f. In addition, the instability is strongest (largest growth rate) for perturbations with azimuthal mode number ℓ=1, corresponding to a simple (dipole) transverse displacement of the beam ions and the background electrons. For the case of overlapping step-function density profiles for the beam ions and background electrons, corresponding to monoenergetic ions and electrons, a key result is that there is no threshold in beam intensity ω[over ^]_{pb}^{2}/ω_{βb}^{0^{2}} or fractional charge neutralization f for the onset of instability. Finally, for the case of continuously varying density profiles with parabolic profile shape, a semiquantitative estimate is made of the effects of the corresponding spread in (depressed) betatron frequency on stability behavior, including an estimate of the instability threshold for the case of weak density nonuniformity.

Sawtooth stability in neutral beam heated plasmas in TEXTOR

Nonlinear Hybrid Simulation of Internal Kink with Beam Ion Effects in DIII-D

Second Stable Regime of Internal Kink Modes Excited by Barely Passing Energetic Ions in Tokamak Plasmas

Neutral beam stabilization of sawtooth oscillations in JET

M3D-K Simulations of Sawteeth and Energetic Particle Transport in Tokamak Plasmas

Control oriented modeling and simulation of the sawtooth instability in nuclear fusion tokamak plasmas

Fast ion confinement and stability in a neutral beam injected reversed field pinch

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Neutral beam effects on Tokamak ballooning mode stability

Counteracting sawtooth crash effects via fluctuation-induced inward transport in HL-2A NBI plasma

A new explanation of the sawtooth phenomena in tokamaks

Sawtooth suppression by flux pumping on HBT-EP

Effects of Sawtooth Heat Pulses on Edge Flows and Turbulence in a Tokamak Plasma

Observation of Off-Axis Sawtooth Oscillations During the Presence of Nonlinear Mode Couplings in HL-2A NBI Heated Plasmas

Ion-temperature-gradient Instability in Toroidal Plasmas

Isotope effects on transport characteristics of edge and core plasmas heated by neutral beam injection (NBI) in an inward shifted configuration at the Large Helical Device

Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma

Current Profile Redistribution Driven by Neutral Beam Injection in A Reversed-Field Pinch

Modelling of sawtooth-induced fast ion transport in positive and negative triangularity in TCV

Kinetic Electromagnetic Instabilities in an ITB Plasma with Weak Magnetic Shear

Kinetic Description of Electron-Proton Instability in High-Intensity Proton Linacs and Storage Rings Based on the Vlasov-Maxwell Equations