Numerical simulation of neoclassical currents, parallel viscosity, and radial current balance in tokamak plasmas

No Thumbnail Available
Journal Title
Journal ISSN
Volume Title
Doctoral thesis (article-based)
Checking the digitized thesis and permission for publishing
Instructions for the author
Degree programme
45, [65]
Helsinki University of Technology publications in engineering physics. A, 807
One of the principal problems en route to a fusion reactor is that of insufficient plasma confinement, which has lead to both theoretical and experimental research into transport processes in the parameter range relevant for fusion energy production. The neoclassical theory of tokamak transport is well-established unlike the theory of turbulence driven anomalous transport in which extensive progress has been made during last few years. So far, anomalous transport has been dominant in experiments, but transport may be reduced to the neoclassical level in advanced tokamak scenarios. This thesis reports a numerical study of neoclassical fluxes, parallel viscosity, and neoclassical radial current balance in tokamaks. Neoclassical parallel viscosity and particle fluxes are simulated over a wide range of collisionalities, using the fully kinetic five-dimensional neoclassical orbit-following Monte Carlo code ASCOT. The qualitative behavior of parallel viscosity derived in earlier analytic models is shown to be incorrect for high poloidal Mach numbers. This is because the poloidal dependence of density was neglected. However, in high Mach number regime, it is the convection and compression terms, rather than the parallel viscosity term, that are shown to dominate the momentum balance. For fluxes, a reasonable agreement between numerical and analytical results is found in the collisional parameter regime. Neoclassical particle fluxes are additionally studied in the banana regime using the three-dimensional Fokker-Planck code DEPORA, which solves the drift-kinetic equation with finite differencing. Limitations of the small inverse aspect ratio approximation adopted in the analytic theory are addressed. Assuming that the anomalous transport is ambipolar, the radial electric field and its shear at the tokamak plasma edge can be solved from the neoclassical radial current balance. This is performed both for JET and ASDEX Upgrade tokamaks using the ASCOT code. It is shown that shear high enough for turbulence suppression can be driven at the Low (L) to High (H) transition conditions without taking into account anomalous processes. In agreement with experiments, simulations indicate a higher threshold temperature for the L–H transition in JET than in ASDEX Upgrade. The parametric dependence of the shear on temperature, density, and magnetic field, however, is similar for both devices. In agreement with some theoretical models and experimental observations, the results also suggest that the critical shear for strong turbulence suppression in JET should be lower than in ASDEX Upgrade.
fusion, plasma, tokamak, neoclassical transport, Monte Carlo, L-H transition, Mach
Other note
  • T.P. Kiviniemi and J.A. Heikkinen, Choice of constants of motion coordinates in numerical solving of the three-dimensional Fokker-Planck equation for tokamaks, Computer Physics Communications 107 (1997) pp. 149-154. [article1.pdf] © 1997 Elsevier Science. By permission.
  • T.P. Kiviniemi, J.A. Heikkinen, A.G. Peeters, Test particle simulation of nonambipolar ion diffusion in tokamaks, Nuclear Fusion 40 (2000) pp. 1587-1596. [article2.pdf] © 2000 IAEA. By permission.
  • T.P. Kiviniemi, J.A. Heikkinen, A.G. Peeters, Effect of poloidal density variation on parallel viscosity for large Mach numbers, Physics of Plasmas 7 (2000) pp. 5255-5258. [article3.pdf] © 2000 American Institute of Physics. By permission.
  • J.A. Heikkinen, T.P. Kiviniemi, A.G. Peeters, T. Kurki-Suonio, S.K. Sipilä, W. Herrmann, W. Suttrop, H. Zohm, Ion orbit loss current in ASDEX Upgrade, Plasma Physics and Controlled Fusion 40 (1998) pp. 693-696. [article4.pdf] © 1998 Institute of Physics Publishing. By permission.
  • J.A. Heikkinen, T.P. Kiviniemi, A.G. Peeters, Neoclassical radial current balance in tokamaks and transition to the H-mode, Physical Review Letters 84 (2000) pp. 487-490. [article5.pdf] © 2000 American Physical Society. By permission.
  • T.P. Kiviniemi, J.A. Heikkinen, A.G. Peeters, S.K. Sipilä, Monte Carlo guiding-centre simulations of ExB flow shear in edge transport barrier, accepted for publication in Plasma Physics and Controlled Fusion. [article6.pdf] © 2001 Institute of Physics Publishing. By permission.
  • T.P. Kiviniemi, T. Kurki-Suonio, S.K. Sipilä, J.A. Heikkinen, A.G. Peeters, L-H transport barrier formation: self-consistent simulation and comparison with ASDEX Upgrade experiments, Czechoslovak Journal of Physics 49, No. 12, Supplement S3, (1999) pp. 81-92. [article7.pdf] © 1999 Czechoslovak Journal of Physics. By permission.
Permanent link to this item