Particle Simulation of the Neoclassical Plasmas

Jukka Heikkinen, T.P. Kiviniemi, T. Kurki-Suonio, A.G. Peeters, S.K. Sipilä

Research output: Contribution to journalArticleScientificpeer-review

72 Citations (Scopus)


A 5D Monte Carlo particle simulation method for advancing rotating plasmas in tori is presented. The method exploits the neoclassical radial current balance (quasineutrality condition). Including the ion polarization current gives the time rate of change of the radial electric field and related evolution of the rotation velocity components. A special orbit initialization for a quiescent start and an efficient radial flux solving algorithm with reduced numerical noise are developed. Numerical stability of the method with respect to the strength of the perpendicular viscosity and Mach number of the poloidal rotation is investigated. This new approach enables one to separate the nonambipolar transport characteristics from the ambipolar ones. Because nonambipolar transport can support sheared flows, this model can provide a very efficient tool for studying transport barriers and related neoclassical mechanisms in toroidal plasmas.
Original languageEnglish
Pages (from-to)527-548
JournalJournal of Computational Physics
Issue number2
Publication statusPublished - 2001
MoE publication typeA1 Journal article-refereed


  • Monte Carlo simulation
  • plasma
  • transport


Dive into the research topics of 'Particle Simulation of the Neoclassical Plasmas'. Together they form a unique fingerprint.

Cite this