ITER Relevant Simulations of Lower Hybrid and Ion Cyclotron Waves with Self-Consistent Non-Maxwellian Species

The next step toward fusion as a practical energy source is the design and construction of ITER. ITER relies in part on ion-cyclotron radio frequency (ICRF) power to heat the deuterium-tritium plasma to fusion temperatures as well as to provide a portion of the current during flat-top operations. Lower hybrid (LH) RF power is under consideration as an upgrade to the baseline heating and current drive system in order to provide control over the current profile and additional current during start up. We have applied a suite of mature radio frequency (rf) full wave codes using self-consistent particle distributions from a Fokker-Planck code to ITER ICRF scenarios and ITER relevant LH experiments on Alcator C-Mod. We will present three dimensional fullwave simulations showing that the ICRF waves propagate radially inward in ITER with strong central focusing and little toroidal spreading. Fokker-Planck coupled rf simulations show that because of the high plasma density, energetic ion tail formation in ITER is typically weak, with the exception of the minority deuterium heating scheme where strong tails can develop on the minority ion distribution. Absorption by the fast alpha particles can approach five to ten percent of the injected power. Massively parallel full wave simulations in the lower hybrid range of frequencies using 2000 poloidal modes and 1000 radial elements have shown that proper reconstruction of wave fronts in the full wave treatment at caustics and cut-offs, where WKB methods fail, can lead to significant spectral broadening. We demonstrate that this linear mechanism is sufficient to bridge the spectral gap (the difference between the high injected phase velocities and the slower phase velocity at which damping on electrons occurs) and explains the efficient damping of lower hybrid waves. This is seen to affect the amount of broadening in the phase velocity and the current drive location.

[1]  M. Brambilla,et al.  Full-wave simulations of lower hybrid wave propagation in toroidal plasma with nonthermal electron distributions , 2008 .

[2]  M. Brambilla,et al.  Full‐wave Simulations of LH wave propagation in toroidal plasma with non‐Maxwellian electron distributions , 2007 .

[3]  Riccardo Maggiora,et al.  Evolution of nonthermal particle distributions in radio frequency heating of fusion plasmas , 2007 .

[4]  N. Bertelli,et al.  Beam Tracing description of LH waves in tokamaks , 2007 .

[5]  Marco Brambilla,et al.  Full Wave Simulations of Fast Wave Mode Conversion and Lower Hybrid Wave Propagation in Tokamaks , 2004 .

[6]  C. Chapman,et al.  Application of the Maslov Seismogram Method in Three Dimensions , 2002 .

[7]  R. Aymar,et al.  Overview of ITER-FEAT - The future international burning plasma experiment , 2001 .

[8]  Eduardo F. D'Azevedo,et al.  All-orders spectral calculation of radio-frequency heating in two-dimensional toroidal plasmas , 2001 .

[9]  Y Lin,et al.  Plasma curvature effects on microwave reflectometry fluctuation measurements , 2001 .

[10]  R. Nazikian,et al.  LETTER TO THE EDITOR: Plasma curvature effects on microwave reflectometry fluctuation measurements , 2001 .

[11]  Current Drive Chapter 6: Plasma auxiliary heating and current drive , 1999 .

[12]  Marco Brambilla,et al.  Numerical simulation of ion cyclotron waves in tokamak plasmas , 1999 .

[13]  F. G. Rimini,et al.  D-T fusion with ion cyclotron resonance heating in the JET tokamak , 1998 .

[14]  X. Litaudon,et al.  Statistical theory of wave propagation and multipass absorption for current drive in tokamaks , 1993 .

[15]  R. W. Harvey,et al.  The CQL3D Fokker-Planck code , 1992 .

[16]  G. V. Pereverzev,et al.  Use of the multidimensional WKB method to describe propagation of lower hybrid waves in tokamak plasma , 1992 .

[17]  F. Berz On the Theory of Plasma Waves , 1956 .