Core Concepts
Ein neues symmetrisches und direkt energieerhaltenes Verfahren (CIDG-C) wird entwickelt, um die Dynamik geladener Teilchen in elektromagnetischen Feldern numerisch zu lösen.
Abstract
In dieser Arbeit wird ein neuartiges direkt energieerhaltenes Verfahren, das sogenannte CIDG-C-Verfahren, zur numerischen Lösung der Lorentz-Kraft-Gleichung entwickelt. Die Lorentz-Kraft-Gleichung kann als nicht-kanonisches Hamiltonsystem geschrieben werden, sodass geometrische Integrationsverfahren angewendet werden können.
Das CIDG-C-Verfahren basiert auf der Koordinateninkrement-diskreten-Gradienten-Methode und ist symmetrisch sowie direkt energieerhaltend, ohne numerische Quadraturformeln zu verwenden. Es wird gezeigt, dass das CIDG-C-Verfahren im Vergleich zum weit verbreiteten Boris-Verfahren Vorteile bei der Energieerhaltung aufweist.
Numerische Experimente für verschiedene Testfälle, darunter die Bewegung eines geladenen Teilchens in statischen elektromagnetischen Feldern und in der Geometrie eines Tokamak-Reaktors, bestätigen die Leistungsfähigkeit des CIDG-C-Verfahrens. Das neue Verfahren kann die Hamiltonfunktion exakt erhalten, während andere Methoden wie das Boris-Verfahren oder Linienintegral-Verfahren nur eine approximative Energieerhaltung aufweisen.
Stats
Die Hamiltonfunktion des Systems ist H = 1/2 v^T v + U(x).
Die Winkelimpulserhaltung ist gegeben durch p = R^2 d/dt ξ + 1/3 R^3.
Das magnetische Moment μ = v_⊥^2 / (2 R) ist eine adiabatische Invariante.
Quotes
"Das CIDG-C-Verfahren ist symmetrisch und kann die Hamiltonfunktion des Systems exakt erhalten."
"Im Vergleich zum weit verbreiteten Boris-Verfahren weist das CIDG-C-Verfahren Vorteile bei der Energieerhaltung auf."