280 likes | 446 Views
Computer simulations of fast frequency sweeping mode in JT-60U and fishbone instability. Y. Todo (NIFS) Y. Shiozaki (Graduate Univ. Advanced Studies) K. Shinohara, M. Takechi, M. Ishikawa (JAERI). October 6-8, 2003
E N D
Computer simulations of fast frequency sweeping mode in JT-60U and fishbone instability Y. Todo (NIFS) Y. Shiozaki (Graduate Univ. Advanced Studies) K. Shinohara, M. Takechi, M. Ishikawa (JAERI) October 6-8, 2003 8th IAEA Technical Meeting on Energetic Particles in Magnetic Confinement Systems General Atomics, San Diego
Fast Frequency Sweeping Mode observed in the JT-60U plasma with NNB injection K. Shinohara et al., Nucl. Fusion 41, 603 (2001). Frequency sweeping takes place both upward and downward by 10-20kHz in 1-5 ms.
Particle-MHD hybrid simulation[Y. Todo and T. Sato, Phys. Plasmas 5, 1321 (1998)] • Plasma is divided into “energetic ions” + “MHD fluid”. • Electromagnetic fields are given by MHD equations. • Energetic ions are described by the drift-kinetic equation.
Equilibrium • An extended Grad-Shafranov equation for anisotropic energetic ion pressure is found (equivalent to Eq (9) in [Belova et al. Phys. Plasmas 10, 3240 (2003)]). We solve a stream function K in the zero orbit width limit.
Numerical method • MHD: finite difference with 4th order accuracy in space and time. • Number of grid points: (101, 16-48, 101) for (R, j, z) coordinates • Time steps: (2-7)105 • The df particle simulation is applied to energetic ions. • Number of particles N: 5105-4106 • Retaining fluctuations with toroidal mode number n=1 to reduce numerical noise (N=5105). Numerical convergence is checked with a simulation where toroidal mode number n=0-6 is retained and N= 4106.
Condition (E36379, t=4.0s) • R0 = 3.4 m, a = 1.0 m • B = 1.2 T • DD plasma • ne at the plasma center 1.6 1019 m-3 • Alfvén velocity at the plasma center 4.6 106 m/s • Injection energy 346 keV (5.75 106 m/s) • Only parallel velocity is considered. • Maximum velocity of the slowing down distribution is assumed to be 80% of the injection velocity. • Finite viscosity and resistivity (10-5 or 210-5 R0vA) are assumed for nonlinear simulation.
Investigation for different energetic ion pressure profiles • Red curve is a classical distribution based on the OFMC code calculation. • The energetic ion beta value at the plasma center is 1.9% for the red curve. • For the other curves central beta value is reduced. • All the profiles are similar at r/a>0.6.
q-profile and shear Alfvén continuous spectra Toroidal mode number n=1 q(0)=1.6 q=2.5 at r/a=0.8 Mode spatial profile depends on energetic ion pressure profile. They do not peak at the TAE gap at q=2.5. ------> Energetic Particle mode (EPM)
Mode spatial profile (toroidal electric field) depends on energetic ion pressure profile Profile I: bh(0)=1.9% Profile III: bh(0)=0.8%
A new kind of EPM The EPMs are not on the shear Alfvén continuum. They are different from the local (resonant) EPM [L. Chen (1994), C.Z. Cheng et al. (1995)]. Nonlocal EPM is studied for reversed shear plasma and trapped energetic ions [H.L. Berk (2001), F. Zonca (2002)]. The EPMs found in this work is similar to the nonlocal EPM but in monotonic shear plasma and with passing energetic ions. Theoretical exploration is needed!
Investigation on different q profileswith energetic ion pressure profile I • The safety factor q is uniformly reduced by 0.2; q=qexp-0.2. • A gap at q=1.5 (r/a~0.3) appears. • A local (resonant) EPM is destabilized. • Frequency is on the Alfvén continuum.
Nonlocal EPM and local (resonant) EPM Test q profile (q=qexp-0.2) Local (resonant) EPM Spiral structure q profile from experiment Nonlocal EPM
Nonlocal EPM and local (resonant) EPM (II) Cosine part of n=1 electrostatic potential. Sine part of n=1 electrostatic potential. Local (resonant) EPM Nonlocal EPM
Time evolution of the nonlocal EPM(energetic ion pressure profile I: classical) Left: Time evolution of cosine part of the toroidal electric field. dBmax/B~10-2. Right: Frequency shifts only downward by 33% (~17kHz) in 6 102 Alfvén time (~0.5ms).
Consider reduced energetic ion pressure profiles • In the experiments fast FS modes and abrupt large events take place with time intervals much shorter than the slowing down time. • These activities lead to redistribution and loss of energetic ions. • The classical distribution gives an overestimate of energetic ion pressure profile. • We consider reduced pressure profiles whose shape is the same as the classical distribution.
Time evolution of the nonlocal EPM(classical distribution is reduced to 2/5) Left: Time evolution of cosine part of the toroidal electric field. dBmax/B~210-3 Right: Frequency shifts upward by 14% (~7kHz) and downward by 23% (~12kHz) in 103 Alfvén time (~0.8ms). Close to the experiment.
Radial beta profiles Initial: reduced pressure profile. 15% reduction in central beta. Flattened only at r/a<0.2. Initial: classical pressure profile. 40% reduction in central beta. Flattened at r/a<0.4.
Time evolution of mode spatial profile(toroidal electric field) Classical pressure profile At=84 At=252 At=335 Reduced pressure profile Change in mode profile is small. At=252 At=418 At=924
Comparison with hole-clump pair creation[H.L. Berk et al. Phys. Lett. A 234, 213 (1997); 238, 408(E) (1998).] Blue circles: Frequency upshift and downshift. Red triangles: Only downshift. Consistent with the hole-clump pair creation. For gd/wA=0.025, gL/wA=0.05, the theory predicts frequency shift dw=0.44 gL(gdt)1/2=0.11 wA in 103 Alfvén time. This corresponds to 22kHz. Simulation result is 10kHz. In good agreement with a factor of 2.
Summary for fast FS mode A new kind of nonlocal EPM is found near the plasma center with monotonic magnetic shear and passing energetic ions. When a classical distribution is taken for the initial condition, frequency shifts only downward. When a reduced distribution is taken, frequency shifts both upward and downward. The rate of frequency sweeping is close to the fast FS modes in the experiments. Frequency shifts both upward and downward are consistent with the hole-clump pair creation. A great redistribution of energetic ions takes place when the classical distribution is taken for the initial condition.
Computer Simulation of Fishbone oscillation (Y. Shiozaki) Clarify self-consistent nonlinear evolution of precessional fishbone oscillations; spatial profile, frquency evolution, saturation mechanism, energetic ion transport GOAL METHOD Particle + MHD hybrid simulation (current coupling) Bulk plasma : Ideal MHD equation Fast ion : kinetic equation using delta-f method BACKGROUND Fishbone instabilities was first identified during near-perpendicular NBI injection in PDX tokamak.
Initial condition of simulation 1: Major radius 1.43 [m] Minor radius 0.43 [m] 2: Magnetic field at axis 1.5 [T] H+ plasma D+ plasma 3: Bulk plasma component Beam plasma component 4: Density 51019 [m-3] 5: Velocity distribution of NBI fast ion isotropic slowing down distribution, injection energy 50keV 6: Spatial profile of fast ion pressure Gaussian distribution Pressure gradient scale at q=1 surface 0.15 [m]
Frequency and growth rate at various energetic ion beta Case 1 Case 2 When fast ion beta is 1.0%, the frequency is oscillating between 4kHz and –4kHz For relatively high energetic ion beta values fishbone mode is destabilized This result is similar to [W. Park et al., Phys. Plasmas 6, 1796 (1999)].
Mode Structure of kink mode q-profile sin component of radial velocity Case 1 cos component of radial velocity Radial velocity of (1,1) mode is well-known top-hat structure. q Mode Structure of fishbone Case 2 sin component of radial velocity q-profile cos component of radial velocity Around q=1 surface phase of the mode profile changes (spiral structure). Two step structure is slightly seen. Higher resolution is needed to clarify the two step structure. q
Temporal evolution of poloidal magnetic field and frequency of the fishbone mode, bh=3% Frequency shifts downward to zero. Saturation level dB/B~510-3
Time evolution in each toroidal mode number MHD energy ( toroidal mode number N=0 ) MHD energy ( toroidal mode number N=1 ) MHD energy ( toroidal mode number N=2 ) The n=1 mode saturates at 1400 Alfvén time. The n=0 and n=2 modes continues to grow after the n=1 mode saturation. MHD Energy Transferred Energy (from energetic particles to MHD) Total energy Energy transfer from energetic ions to MHD continues even after the saturation of the n=1 mode growth.
Summary for fishbone mode • Linear stability and nonlinear evolution of fishbone mode is investigated. Frequency shifts downward at saturation. Saturation level is dB/B~510-3. • After saturation of the n=1 mode, n=0 and n=2 modes continues to grow. Energy transfer from energetic ions to MHD fluid keeps a constant level. • Future work: 3.1 investigate energetic ion transport. 3.2 investigate generation of nonlinear MHD modes, interaction between energetic ions and the nonlinear MHD modes.