Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates
The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component BΨP of the magnetic field vector with respect to the poloidal flux function ΨP. In an unperturbed plasma, BΨP can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of BΨP is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when BΨP is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without BΨP yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with BΨP retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE; Japan Society for the Promotion of Science (JSPS)
- Grant/Contract Number:
- AC02-09CH11466; JPMXP1020200103
- OSTI ID:
- 1957529
- Alternate ID(s):
- OSTI ID: 1899410
- Journal Information:
- Physics of Plasmas, Vol. 29, Issue 11; ISSN 1070-664X
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Analytical calculation of Boozer magnetic coordinates for axisymmetric magnetohydrodynamic equilibria
Full-field drift Hamiltonian particle orbits in axisymmetric tokamak geometry