Automation of the guiding center expansion
- Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543 (United States)
We report on the use of the recently developed Mathematica package VEST (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation's numerous non-unique forms or representations. (2) The procedure proceeds entirely in cartesian position and velocity coordinates, thereby producing manifestly gyrogauge invariant results; the commonly used perpendicular unit vector fields e{sub 1},e{sub 2} are never even introduced. (3) It is easy to apply in the derivation of higher-order contributions to the guiding center transformation without fear of human error. Our code therefore stands as a useful tool for exploring subtle issues related to the physics of toroidal momentum conservation in tokamaks.
- OSTI ID:
- 22227924
- Journal Information:
- Physics of Plasmas, Vol. 20, Issue 7; Other Information: (c) 2013 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 1070-664X
- Country of Publication:
- United States
- Language:
- English
Similar Records
VEST: Abstract Vector Calculus Simplification in Mathematica
From liquid crystal models to the guiding-center theory of magnetized plasmas