1 / 55

Energetic Particle Physics in Tokamak Plasmas

Energetic Particle Physics in Tokamak Plasmas. Guoyong Fu Princeton Plasma Physics Laboratory Princeton University Princeton, NJ 08543, USA. 6 th Workshop on Nonlinear Plasma Science and IFTS Workshop Suzh o u, China, S ept. 24. Example of EPM: fishbone instability.

brinda
Download Presentation

Energetic Particle Physics in Tokamak Plasmas

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Energetic Particle Physics in Tokamak Plasmas Guoyong Fu Princeton Plasma Physics Laboratory Princeton University Princeton, NJ 08543, USA 6th Workshop on Nonlinear Plasma Science and IFTS Workshop Suzhou, China, Sept. 24

  2. Example of EPM: fishbone instability Mode structure is of (m,n)=(1,1) internal kink; Mode is destabilized by energetic trapped particles; Mode frequency is comparable to trapped particles’ precessional drift frequency K. McGuire, R. Goldston, M. Bell, et al. 1983, Phys. Rev. Lett.50, 891 L. Chen, R.B. White and M.N. Rosenbluth 1984, Phys. Rev. Lett.52, 1122

  3. First observation of TAE in TFTR . K.L. Wong, R.J. Fonck, S.F. Paul, et al. 1991, Phys. Rev. Lett.66, 1874

  4. RSAE (Alfven cascades) were observed in JET plasmas

  5. NSTX observes that multi-mode TAE bursts can lead to larger fast-ion losses than single-mode bursts 1% neutron rate decrease: 5% neutron rate decrease: • TAE avalanches cause enhanced fast-ion losses. • Potential to model island overlap condition with full diagnostic set. E. Fredrickson, Phys. Plasmas 13, 056109 (2006)

  6. Outline • Roles of energetic particles in fusion plasmas • Single particle confinement • MHD limit: Shear Alfven continuum and Alfven eigenmodes • Linear Kinetic Stability: discrete AE and EPM • Nonlinear Physics: saturation mechanisms, frequency chirping, multi-mode coupling • Prospect for ITER and DEMO • Important problems for future

  7. Roles of energetic particles in fusion plasmas • Heat plasmas via Coulomb collision, drive plasma rotation and plasma current (e.g., NBI injection) • Stabilize MHD modes (e.g., internal kink, RWM) • Destabilize shear Alfven waves via wave-particle resonance • Energetic particle loss can damage reactor wall • Energetic particle redistribution may affect thermal plasma via plasma heating profile, plasma rotation and current drive. • Energetic particle-driven instability can be beneficial (alpha channeling, diagnostic of q profile etc)

  8. Single Particle Confinement • For an axi-symmetric torus, energetic particles are well confined (conservation of toroidal angular momentum). • Toroidal field ripple can induce stochastic diffusion (important in advanced plasma regime with high q) • Symmetry-breaking MHD modes can also cause energetic particle anomalous transport.

  9. Shear Alfven spectrum and continuum damping • Shear Alfven wave dispersion relation and continuum spectrum

  10. Discrete Alfven Eigenmodes can exist near continuum accumulation point due to small effects such as toroidicity, shaping, magnetic shear, and energetic particle effects. Coupling of m and m+k modes breaks degeneracy of Alfven continuum : K=1 coupling is induced by toroidicity K=2 coupling is induced by elongation K=3 coupling is induced by triangularity. or EPM ----- RSAE

  11. Discrete Shear Alfven Eigenmodes • Cylindrical limit  Global Alfven Eigenmodes • Toroidal coupling  TAE and Reversed shear Alfven eigenmodes • Elongation  EAE and Reversed shear Alfven eigenmodes • Triangularity  NAE • FLR effectsKTAE

  12. Toroidal Alfven Eigenmode (TAE) can exist inside continuum gap TAE mode frequencies are located inside the toroidcity-induced Alfven gaps; TAE modes peak at the gaps with two dominating poloidal harmonics. C.Z. Cheng, L. Chen and M.S. Chance 1985, Ann. Phys. (N.Y.)161, 21

  13. Reversed shear Alfven eigenmode (RSAE) can exist above maximum of Alfven continuum at q=qmin U wA q wRSAE rmin r rmin r rmin r w = (n-m/qmin)/R

  14. RSAE exists due to toroidicity, pressure gradient or energetic particle effects H.L. Berk, D.N. Borba, B.N. Breizman, S.D. Pinches and S.E. Sharapov 2001, Phys. Rev. Lett.87 185002 S.E. Sharapov, et al. 2001, Phys. Lett. A289, 127 B.N. Breizman et al, Phys. Plasmas 10, 3649 (2003) G.Y. Fu and H.L. Berk, Phys. Plasmas 13,052502 (2006)

  15. Linear Stability • Energetic particle destabilization mechanism • Kinetic/MHD hybrid model • TAE stability: energetic particle drive and dampings • EPM stability: fishbone mode

  16. Destabilize shear Alfven waves via wave-particle resonance • Destabilization mechanism (universal drive) Wave particle resonance at For the appropriate phase (n/w >0), particles will lose energy going outward and gain energy going inward. As a result, particles will lose energy to waves. Energetic particle drive Spatial gradient drive Landau damping (or drive) due to velocity space gradient

  17. Kinetic/MHD Hybrid Model

  18. Quadratic form

  19. Drift-kinetic Equation for Energetic Particle Response

  20. Perturbative Calculation of Energetic Particle Drive G.Y.Fu and J.W. Van Dam, Phys. Fluids B1, 1949 (1989) R. Betti et al, Phys. Fluids B4, 1465 (1992).

  21. Dampings of TAE • Ion Landau damping • Electron Landau damping • Continuum damping • Collisional damping • “radiative damping” due to thermal ion gyroradius G.Y.Fu and J.W. Van Dam, Phys. Fluids B1, 1949 (1989) R. Betti et al, Phys. Fluids B4, 1465 (1992). F. Zonca and L. Chen 1992, Phys. Rev. Lett.68, 592 M.N. Rosenbluth, H.L. Berk, J.W. Van Dam and D.M. Lindberg 1992, Phys. Rev. Lett.68, 596 R.R. Mett and S.M. Mahajan 1992, Phys. Fluids B4, 2885

  22. Fishbone dispersion relation as an example of EPM For EPM, the energetic particle effects are non-perturbative and the mode frequency is determined by particle orbit frequency. L. Chen, R.B. White and M.N. Rosenbluth 1984, Phys. Rev. Lett.52, 1122

  23. Discrete Alfven Eigenmodes versus Energetic Particle Modes • Discrete Alfven Eigenmodes (AE): Mode frequencies located outside Alfven continuum (e.g., inside gaps); Modes exist in the MHD limit; Are typical weakly damped due to thermal plasma kinetic effects; energetic particle effects are often perturbative. • Energetic Particle Modes (EPM): Mode frequencies located inside Alfven continuum and determined by energetic particle dynamics; Energetic effects are non-perturbative; Requires strong energetic particle drive to overcome continuum damping.

  24. Nonlinear Physics • Nonlinear dynamics of a single mode Bump-on-tail problem: saturation due to wave particle trapping frequency chirping • Multiple mode effects: mode-mode coupling, resonance overlap

  25. Bump-on-tail problem: definition It can be shown that bump-on-tail problem is nearly equivalent to that of energetic particle-driven instability. H.L. Berk and B.N. Breizman 1990, Phys. Fluids B2, 2235

  26. saturation due to wave particle trapping We first consider case of no source/sink and no damping. The instability then saturates at The instability saturates when the distribution is flattened at the resonance region Width of flattened region is on order of

  27. Bump-on-tail problem: steady state saturation with damping, source and sink Collisions tend to restore the original unstable distribution. Balance of nonlinear flattening and collisional restoration leads to mode saturation. It can be shown that the linear growth rate is reduced by a factor of . Thus, the mode saturates at H.L. Berk and B.N. Breizman 1990, Phys. Fluids B2, 2235

  28. H.L. Berk et al, Phys. Plasmas 2, 3007 (1995).

  29. Transition from steady state saturation to explosive nonlinear regime (near threshold) B.N. Breizman et al, Phys. Plasmas 4, 1559(1997).

  30. Hole-clump creation and frequency chirping • For near stability threshold and small collision frequency, hole-clump will be created due to steepening of distribution function near the boundary of flattening region. • As hole and clump moves up and down in the phase space of distribution function, the mode frequency also moves up and down. H.L. Berk et al., Phys. Plasma 6, 3102 (1999).

  31. Experimental observation of frequency chirping M.P. Gryaznevich et al, Plasma Phys. Control. Fusion 46 S15, 2004.

  32. Saturation due to mode-mode coupling • Fluid nonlinearity induces n=0 perturbations which lead to equilibrium modification, narrowing of continuum gaps and enhancement of mode damping. D.A. Spong, B.A. Carreras and C.L. Hedrick 1994, Phys. Plasmas1, 1503 F. Zonca, F. Romanelli, G. Vlad and C. Kar 1995, Phys. Rev. Lett.74, 698 L. Chen, F. Zonca, R.A. Santoro and G. Hu 1998, Plasma Phys. Control. Fusion40, 1823 • At high-n, mode-mode coupling leads to mode cascade to lower frequencies via ion Compton scattering. As a result, modes saturate due to larger effective damping. T.S. Hahm and L. Chen 1995, Phys. Rev. Lett.74, 266

  33. Multiple unstable modes can lead to resonance overlap and stochastic diffusion of energetic particles . H.L. Berk et al, Phys. Plasmas 2, 3007 (1995).

  34. Recent Nonlinear Simulation Results from M3D code • Strong frequency chirping of fishbone due to kinetic nonlinearity alone; • Simulations of multiple TAEs in NSTX show strong mode-mode coupling due to resonance overlap; • Nonlinearly generated modes;

  35. M3D XMHD Model

  36. Recent M3D results: • Alpha particle stabilization • of n=1 kink in ITER; • (2) Nonlinear frequency chirping • of fishbone; • (3) beam-driven TAEs in DIII-D; • (4) beam-driven TAEs in NSTX.

  37. Saturation amplitude scale as square of linear growth rate

  38. Simulation of fishbone shows distribution fattening and strong frequency chirping distribution

  39. Evidence of a nonlinearly driven n=2 mode n=1 n=2

  40. Multi-mode simulations show strong mode-mode interaction. n=2 n=3

  41. Strong interaction between different modes is due to wave-particle resonance overlap L=0.2 0.4 Pj 0.6 0.8 v v

  42. Prospect for ITER/DEMO • Multiple high-n modes are expected to be unstable, especially in DEMO; • Key question: (1) how strong is instability? (2) can multiple modes induce global fast ion transport due to resonance overlap and/or avalanche? (3) can fast ion instabilities affect thermal plasma significantly?

  43. Linear Stability of TAE • Alpha particle drive • Plasma dampings: ion Landau damping electron collisional damping “radiative damping” due to FLR  stability is sensitive to plasma parameters and profiles !

  44. Alpha particle drive is maximized at kqra ~ 1 G.Y. Fu et al, Phys. Fluids B4, 3722 (1992)

  45. Parameters of Fusion Reactors

  46. Critical alpha parameters of TFTR/JET, ITER and DEMO DEMO_Japan ba/bc DEMO_EU ARIES-AT ARIES-ST ITER JET TFTR a/ra

  47. Multiple high-n TAEs are expected be excited in ITER from NOVA-K Instability is maximized at kqra ~ 1 N.N. Gorelenkov, Nucl. Fusion 2003

More Related