Monte Carlo calculation of the neoclassical transport matrix by the Einstein-Helfand relation in nonaxisymmetric toroidal plasmas

Journal Article · Mon Mar 15 00:00:00 EDT 2010 · Physics of Plasmas

DOI:https://doi.org/10.1063/1.3299366· OSTI ID:21347180

Matsuyama, A ^[1]; Hanatani, K ^[2]

---

1. Graduate School of Energy Science, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan)
2. Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan)

A Monte Carlo method that calculates the neoclassical transport matrix with the pitch-angle scattering approximation is described for nonaxisymmetric toroidal plasmas. The numerical scheme is based on the Einstein-Helfand relation, which generalizes the calculation of the neoclassical diffusion coefficient with the mean-square displacement. By calculating all the elements of the transport matrix simultaneously, not only the neoclassical diffusion, but also the parallel, poloidal, and toroidal viscosity coefficients can be evaluated for three-dimensional magnetic-field configurations. The newly developed Monte Carlo code is benchmarked favorably for an l=2 single-helicity model against the results of the variational principle.

• OSTI ID: 21347180
• Journal Information: Physics of Plasmas, Vol. 17, Issue 3; Other Information: DOI: 10.1063/1.3299366; (c) 2010 American Institute of Physics; ISSN 1070-664X
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
DIFFUSION
MAGNETOHYDRODYNAMICS
MATRICES
MONTE CARLO METHOD
NEOCLASSICAL TRANSPORT THEORY
PLASMA
PLASMA SIMULATION
VARIATIONAL METHODS
VISCOSITY
CALCULATION METHODS
CHARGED-PARTICLE TRANSPORT THEORY
FLUID MECHANICS
HYDRODYNAMICS
MECHANICS
SIMULATION
TRANSPORT THEORY