180 likes | 323 Views
Nonlinear plasma - wave interactions in ion cyclotron range of frequency. N Xiang, C. Y Gan , J. L. Chen, D. Zhou Institute of plasma phsycis , CAS, Hefei J. R. Cary University of Colorado and Tech-X corp. outline. Motivation Benchmark with linear theory Parametric decays of ICRF
E N D
Nonlinear plasma-wave interactions in ion cyclotron range of frequency N Xiang, C. Y Gan, J. L. Chen, D. Zhou Institute of plasma phsycis, CAS, Hefei J. R. Cary University of Colorado and Tech-X corp
outline • Motivation • Benchmark with linear theory • Parametric decays of ICRF • Induced particle transport • Summary
(1) ICRF wave heating has been widely used and proved to be an efficient way to heat plasma, especially for EAST tokamak, key for long pulse H-mode (up to 6MW ICRF heating power). (2) Experiments observed many nonlinear processes of ICRF, especially parametric decay, affect wave absorption. This is beyond the capability of mainstream codes. (3) Particle momentum transport induced by RF waves. Hot topic but theoretical work is very limited . Motivation
Linear theory of mode transform (conversion) of ICRF • An electron plasma wave (EPW) can mode-transformed to an ion Bernstein wave (IBW) (Ono,93) • Near lower hybrid resonance (LHR) (B ~0), Ex reaches its maximum • Nonlinear effects could be strongest near LHR.
Simulation model • Full-PIC simulation in the framework of VORPAL(Chet & Cary 2004). • Fully kinetic ions and electrons. • Density gradient in x (radial), external magnetic field in y (toroidal). Simulation model. Antenna is on the left side. Plasma density is along x-direction.
Simulations confirm theoretical predictions • Mode-transform confirmed • Simulated IBWs satisfy the linear dispersion relation
Parametric decays induced by mode conversion • Decay channels of an IBW, w0 = w1+ w2 1) IBW (w0) -> IBW (w1) + ICQM (w2~ mWi) (non-resonant) 2) IBW (w0) -> IBW (w1) + EQM (w2 ~ k//Vte) 3) IBW (w0) -> IBW (w1) + IBW(w2) (resonant) • Experimental results on HT-7 (J.Li et al 01) (1) and (2) were observed. • Particles can be heated via nonlinear Landau damping (Pokolab 72, Xiang 11)
PDI for H-D plasma with w = 1.2WH • Possible decay channels • IBW (w) -> IBW (w-WD~0.56w)+ ICQM(WD~0.44w) • IBW (w) -> IBW (w+WD~1.44w)+ ICQM(WD~0.44w) • IBW (w) -> IBW (w-WD~-0.3w)+ ICQM(3WD~1.3w) 4) IBW (w) -> QM (w-WH~ 0.18w)+ ICQM(WH~0.8w) 5) IBW (w) -> IBW (w+WH~ 1.8w)+ ICQM(WH~0.8w) • IBW (w) -> IBW (w-WH~ -0.6w)+ ICQM(2WH~1.6w) 7) IBW (w) -> IBW (w+k//Vte~ 1.06w)+ EQM(k//Vte~0.06w)
Simulations confirm predicted PDI • Spectra of Ex, • Pattern (1),(3) ,(4) and (6) are observed. • Two ion-species plasma (nH/nD = 0.05/0.95)
Consequence of PDI: ion heating • Hydrogen and Deuterium ions are heated via NLD • Channels (1)-(3) for D, (4)-(6) for D.
Decay channels decrease for pure H-plasma • For pure hydrogen plasma, only decay channels (4)-(7) are possible. (5) and (6) observed. • H ions are heated.
Decay channels increase for pure D plasma • For pure deuterium plasma, in addition to decay channels (1)-(6) , decay into ICQM with 4WD is observed. • D ions are heated
Asymmetric distribution in poloidal direction is observed • As ion heated, particle distribution in Vx(radial) direction is nearly symmetric, but asymmetric in Vz(poloidal) direction -> ion poloidal flow .
Ion flow induced by RF waves • RF induced poloidal flow predicted by theories (Wang et al 94, Jaeger et al 99, Myra et al 04) comes from the balance of momentum equation in z direction . Reynolds stress electromagnetic force • Simulations indicate that poloidal flow is mainly produced by the force in x direction when neglecting collision terms. Reynolds stress electromagnetic force
Simulations show ion flow shear is induced near LHR • Time-averaged poloidal velocity
Contributions of Reynolds stress and electromagnetic force • Reynolds stress causes charge separation near LHR which in turn produces a radial electric field. • Contribution of electromagnetic force mainly comes from radial electric field and is comparable to the Reynolds stress in the region around LHR.
Summary • PICsimulations are performed to study nonlinear effects during ICRF heating. • Simulations agree with the linear theory for small input power. • Parametric decays of IBWs are observed and ion are heated via nonlinear Landau damping near LH resonance. • Poloidal ion shear flow is produced due to the ion heating. Electric force mainly comes from the radial electric field due the charge separation and it is comparable to Reynolds stress near the LHR.