Variational symplectic algorithm for guiding center dynamics and its application in tokamak geometry

A variational symplectic integrator for the guiding center motion of charged particles in general magnetic fields is developed to enable accurate long-time simulation studies of magnetized plasmas. Instead of discretizing the differential equations of the guiding center motion, the action of the guiding center motion is discretized and minimized to obtain the iteration rules for advancing the dynamics. The variational symplectic integrator conserves exactly a discrete Lagrangian symplectic structure and globally bounds the numerical error in energy by a small number for all simulation time steps. Compared with standard integrators, such as the fourth order Runge–Kutta method, the variational symplectic integrator has superior numerical properties over long integration time. For example, in a two-dimensional tokamak geometry, the variational symplectic integrator is able to guarantee the accuracy for both the trapped and transit particle orbits for arbitrarily long simulation time. This is important for mo...

[1]  T. M. O'Neil,et al.  CROSS-MAGNETIC-FIELD HEAT CONDUCTION IN NON-NEUTRAL PLASMAS , 1997 .

[2]  T. A. Brengle,et al.  An orbit averaged particle code , 1980 .

[3]  D. Dubin.  Collisional transport in non-neutral plasmas , 1998 .

[4]  O'Neil,et al.  Two-dimensional guiding-center transport of a pure electron plasma. , 1988, Physical review letters.

[5]  Robert G. Littlejohn,et al.  Variational principles of guiding centre motion , 1983, Journal of Plasma Physics.

[6]  Hong Qin,et al.  Variational symplectic integrator for long-time simulations of the guiding-center motion of charged particles in general magnetic fields. , 2008, Physical review letters.

[7]  R. Waltz,et al.  Anomalous transport scaling in the DIII-D tokamak matched by supercomputer simulation. , 2003, Physical review letters.

[8]  D. Dubin.  Test Particle Diffusion and the Failure of Integration along Unperturbed Orbits , 1997 .

[9]  Scott E. Parker,et al.  A δ f particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations , 2003 .

[10]  Williams,et al.  Scalings of Ion-Temperature-Gradient-Driven Anomalous Transport in Tokamaks. , 1996, Physical review letters.

[11]  R. Ruth,et al.  Fourth-order symplectic integration , 1990 .

[12]  Wojciech Rozmus,et al.  A symplectic integration algorithm for separable Hamiltonian functions , 1990 .

[13]  J. Marsden,et al.  Discrete mechanics and variational integrators , 2001, Acta Numerica.

[14]  William McCay Nevins,et al.  General Gyrokinetic Equations for Edge Plasmas , 2006 .

[15]  F. Jenko,et al.  Electron temperature gradient turbulence. , 2000, Physical review letters.

[16]  Hideo Sugama,et al.  Kinetic simulation of a quasisteady state in collisionless ion temperature gradient driven turbulence , 2002 .

[17]  William McCay Nevins,et al.  Geometric gyrokinetic theory for edge plasmasa) , 2007 .

[18]  T. Hahm,et al.  Turbulent transport reduction by zonal flows: massively parallel simulations , 1998, Science.

[19]  John M. Dawson,et al.  Fluctuation-induced heat transport results from a large global 3D toroidal particle simulation model , 1996 .

[20]  W. W. Lee,et al.  Gyrokinetic approach in particle simulation , 1981 .

[21]  K. Feng.  Difference schemes for Hamiltonian formalism and symplectic geometry , 1986 .

[22]  Explicit Runge–Kutta integrator with Hamiltonian correction for long-time simulations of guiding-center orbit in tokamak configurations , 2008 .

[23]  Parker,et al.  Gyrokinetic simulation of ion temperature gradient driven turbulence in 3D toroidal geometry. , 1993, Physical review letters.

[24]  C. Scovel,et al.  Symplectic integration of Hamiltonian systems , 1990 .

[25]  Timothy J. Williams,et al.  Gyrokinetic simulations of E×B velocity‐shear effects on ion‐temperature‐gradient modes , 1993 .