• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter UniPHY Group iResearch App Facebook

Phys. Plasmas 16, 042510 (2009); http://dx.doi.org/10.1063/1.3099055 (7 pages)

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

Hong Qin, Xiaoyin Guan, and William M. Tang

Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543, USA

View MapView Map

(Received 20 January 2009; accepted 24 February 2009; published online 21 April 2009)

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 modern large-scale simulation studies of fusion plasmas where it is critical to use algorithms with long-term accuracy and fidelity. The variational symplectic integrator is expected to have a wide range of applications.

© 2009 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. VARIATIONAL SYMPLECTIC INTEGRATOR FOR THE GUIDING CENTER MOTION
  3. APPLICATIONS TO BANANA AND TRANSIT ORBITS IN 2D TOKAMAK GEOMETRY
  4. CONCLUSION AND FUTURE WORK

RELATED DATABASES

To view database links for this article, you need to log in.

KEYWORDS and PACS

PACS

  • 52.20.Dq

    Particle orbits

  • 52.55.Fa

    Tokamaks, spherical tokamaks

  • 52.25.Xz

    Magnetized plasmas

  • 52.65.Cc

    Particle orbit and trajectory

  • 02.60.Lj

    Ordinary and partial differential equations; boundary value problems

  • 02.30.Xx

    Calculus of variations

ARTICLE DATA

PUBLICATION DATA

ISSN

1070-664X (print)  
1089-7674 (online)

For access to fully linked references, you need to log in.
    W. W. Lee, Phys. Fluids 26, 556 (1983)PFLDAS000026000002000556000001.

    B. I. Cohen, T. J. Williams, A. M. Dimits, and J. A. Byers, Phys. Fluids B 5, 2967 (1993)PFBPEI000005000008002967000001.

    S. E. Parker, W. W. Lee, and R. A. Santoro, Phys. Rev. Lett. 71, 2042 (1993).

    A. M. Dimits, T. J. Williams, J. A. Byers, and B. I. Cohen, Phys. Rev. Lett. 77, 71 (1996).

    W. Dorland, F. Jenko, M. Kotschenreuther, and B. N. Rogers, Phys. Rev. Lett. 85, 5579 (2000).

    J. Candy and R. E. Waltz, Phys. Rev. Lett. 91, 045001 (2003).

    T. Watanabe and H. Sugama, Phys. Plasmas 9, 3659 (2003)PHPAEN000009000009003659000001.

    H. Qin and X. Guan, Phys. Rev. Lett. 100, 035006 (2008).

    X. Xiao and S. Wang, Phys. Plasmas 15, 122511 (2008)PHPAEN000015000012122511000001.

    H. Qin, R. H. Cohen, W. M. Nevins, and X. Q. Xu, Phys. Plasmas 14, 056110 (2007)PHPAEN000014000005056110000001.

    D. H. E. Dubin and T. M. O'Neil, Phys. Rev. Lett. 60, 1286 (1988).

    D. H. E. Dubin and T. M. O'Neil, Phys. Rev. Lett. 78, 3868 (1997).

    D. H. E. Dubin, Phys. Rev. Lett. 79, 2678 (1997).

    D. H. E. Dubin and T. M. O'Neil, Phys. Plasmas 5, 1688 (1998)PHPAEN000005000005001688000001.


For access to citing articles, you need to log in.


Figures (4)

Access to article objects (figures, tables, multimedia) requires a subscription; log in to view available files.
(Access to supplementary files, where available, is free for this journal.)



Close
Google Calendar
ADVERTISEMENT

close