Home Mathematics



SPONTANEOUS GENERATION OF HOLES AND CLUMPS IN ENERGETIC PARTICLE DISTRIBUTION BEYOND THE EXPLOSIVE SCENARIO AND FREQUENCYSWEEPING MODESTable of Contents:
It was noted above that the cubic nonlinearity in (8.16) is insufficient for describing the explosive scenario. In a longtime fully nonlinear modelling [8.4], longliving holes and clamps in the energetic ion distribution were observed for the first time. Then, to accommodate all types of the effective collisionality restoring the distribution function, a fully nonlinear ID model ВОТ was developed for the ВОТ instability including the effects of dynamical friction (drag) and velocity space diffusion on the energetic particles driving the wave. The ВОТ results show that in the early nonlinear phase of the instability, the drag facilitates the explosive scenario of the wave evolution. Later, the drag effect leads to the creation of phase space holes and clumps moving away from the original eigenfrequency. The combined effect of drag and diffusion produces a diverse range of nonlinear scenarios, including “hooked” frequency chirping and undulating regimes. Fully Nonlinear 1D BumponTail (ВОТ) Model for Long Time Nonlinear Scenarios beyond the Explosive PhaseThe fully nonlinear system consists of a purely electrostatic wave in a plasma of three species. The first two are the thermal plasma ions and electrons and the third is a low density population of fast electrons that are subject to weak collisions (much less than the background species) and whose distribution function is treated kinetically. In the normalised units of the BerkBreizman theory, the starting set of coupled nonlinear equations has the form [8.6]:
where the collisionality terms are given by (8.8), ^ = kx cot, and u> = e kE jm_{e} is the nonlinear bounce frequency of the electrons trapped in the field of the wave. For describing holes and clumps identified in Ref. [8.4], we represent the frequency in the form iо = a)_{pe} 8+a»(r). As the holes and clumps evolve, the electric field becomes a sum of noninteracting BGK modes with timedependent frequencies, so the envelope takes the form
and the distribution function can be written in a similar manner. To fully explore the nonlinear system (Eq. 8.17 and 8.18), a numerical scheme similar to the one developed in Ref. [8.7] was employed for longterm evolution of the system. More specifically, the Fourier series representation of F in space transformed (8.17) into a set of coupled partial differential equations in t and u. By Fourier transforming in velocity, a set of advection equations was obtained for numerical processing with the ВОТ code. A fixed nearthreshold parameter, l/J/ft^O.9 will be used for ВОТ simulations throughout this section, as Figure 8.8. shows a satisfactory agreement with more demanding case of x_{rf}/yi=0.99. Figure 8.9 shows the result of the ВОТ code for longtime evolution of the nonlinear system (8.17) and (8.18) in the collisionless limit. The holes and clumps move away from the original resonance, as shown schematically in Figure 8.10. This motion is almost adiabatic and preserves the value of the distribution function for FIGURE 8.8 Comparison of the bumpontail simulation with the cubic equation (8.16) in the nearthreshold regime with diffusion for two nearthreshold parameters y_{d} //, . FIGURE 8.9 Spectrogram of the electric field amplitude £, for the collisionless case close to the threshold. The white line is the best l'^{n} fit passing through the upper and lower frequencysweeping structures*. FIGURE 8.10 Cartoon illustrating the motion of holes and clumps and the wake dotted line that steepens the distribution function, creating a favourable environment for instability*. FIGURE 8.11 Spectrogram of the electric field amplitude showing chirping asymmetry for the pure drag case*. Reproduced from [M.K. Lilley et al„ Phys. Plasma 17 (2010) 092305]. with the permission of AIP Publishing. FIGURE 8.12 Spectrogram of the electric field amplitude showing the “hooked” frequency spectrum with drag and diffusion. the particles trapped by the wave. We note that within this simplified theory there is no limit to the extent of the chirp, which is indeed supported by the simulation. The chirping behaviour shows the correct t^{112} scaling [8.4]. The drag introduces a preferred direction of particle flow into the system (from highenergy source to lowenergy sink). Consequently, one can expect the holeclump symmetry observed in Figure 8.9 to be broken when drag is introduced, which is indeed the case as Figure 8.11 shows. In the general case of both drag and diffusion affecting the distribution function at the resonance, more sophisticated patterns of the frequency sweeping are seen. One of this is the “hooked” frequency sweep shown in Figure 8.12. Nonlinear FrequencySweeping Scenarios Observed Experimentally for TAEs Excited with Energetic Ions Resulted from the Drag RelaxationAmong all machines exhibiting the frequencysweeping spectra, spherical tokamaks (START, MAST, NSTX) with superAlfvenic NBI are the best test beds for investigating such phenomena, as discussed in Section 7.3. The exceptional conditions on STs are determined by the low magnetic fields of STs and low values of T„ at the beginning of NBI, so that the following ordering is valid for the beam ions:
where V^t is the beam speed corresponding to the critical energy (3.3), and V_{u} is the initial speed of the injected beam. Under the conditions (8.20), the beam passes through the principal resonance V(, = V_{A} because of the beam slowingdown due to the electron drag. Figure 8.13 schematically shows that a short blip of NBI moves first unidirectionally from the injection highenergy range to low energy due to the drag, without changing its shape because the diffusion effect is small. Only at the low beam energy comparable to £_{clit}, the shape of the beam blip diffuses, and this diffusion is not unidirectional anymore. Because the resonance area E^ comes in this case to the region of the dominant beam relaxation due to the drag, we have a perfect drag scenario similar to that investigated in the ID ВОТ model. Frequencysweeping modes associated with the holeclump generation are inevitable in such a scenario, which is what we observe in STs with superAlfvenic NBI. Figure 8.14 shows magnetic spectrogram from MAST experiment, in which superAlfvenic NBI drives Alfven instability when the resonance Vj_{tKam} = V_{A} is in phase space region dominated by electron drag of the beam ions. It is seen that FS modes dominate the spectrum, with some modes sweeping in frequency to a very long range of 8(i> / ft) = 0.5. Modelling with the HAGIS code [8.8] FIGURE 8.13 Schematic illustration of the beam blip evolution with time in the case of ordering (8.20) typical for STs. FIGURE 8.14 Spectrogram showing FS Alfven modes driven by NBI in MAST discharge #27177. FIGURE 8.15 Nonlinear HAGIS simulation of Alfven instability in MAST #27177. FIGURE 8.16 Top: Mirnov coil data; Middle: magnetic spectrogram; Bottom: Time evolution of tangential energetic neutral spectrum. The NPA viewing angle is set to 0°. and with the beam pure draglike relaxation model installed was performed to compare nonlinear waveparticle dynamics with the experimentally observed in Figure 8.14. Figure 8.15 shows that the HAGIS modelling reproduces pretty well the characteristic spectrum observed in experiments, although the range of the frequency sweeping is not as large as this observed on MAST. Finally, we note that the dominant transport mechanism for nonlinear FS modes is convection of particles trapped in the wave field. Experimentally, the holeclump formation and transport were observed for the first time with an NPA diagnostic on stellarator LHD [8.9]. Figure 8.16 shows how the flux of energetic beam ions sweeps in energy together with the chirping modes. REFERENCES

<<  CONTENTS  >> 

Related topics 