Variational Integration for Ideal MHD with Built-in Advection Equations

Technical Report · Tue Aug 05 00:00:00 EDT 2014

DOI:https://doi.org/10.2172/1179782· OSTI ID:1179782

Zhou, Yao ^[1]; Qin, Hong ^[1]; Burby, J. W. ^[1]; Bhattacharjee, A. ^[1]

---

1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)

Newcomb's Lagrangian for ideal MHD in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.

• Research Organization: Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
• Sponsoring Organization: USDOE Office of Science (SC), Fusion Energy Sciences (FES)
• Contributing Organization: Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543,USA, Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
• DOE Contract Number: DE-AC02-09CH11466
• OSTI ID: 1179782
• Report Number(s): PPPL-5053
• Resource Relation: Related Information: Submitted to Physics of Plasmas
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
Computational Physics
Computer Simulation
Magnetohydrodynamics (MHD)