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

67 Citations (Scopus)

Abstract

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
Volume173
Issue number2
DOIs
Publication statusPublished - 2001
MoE publication typeA1 Journal article-refereed

    Fingerprint

Keywords

  • Monte Carlo simulation
  • plasma
  • transport

Cite this

Heikkinen, J., Kiviniemi, T. P., Kurki-Suonio, T., Peeters, A. G., & Sipilä, S. K. (2001). Particle Simulation of the Neoclassical Plasmas. Journal of Computational Physics, 173(2), 527-548. https://doi.org/10.1006/jcph.2001.6891