[论文解读] Chasing Hamiltonian structure in gyrokinetic theory

J. W. Burby|arXiv (Cornell University)|Sep 13, 2015

Quantum and Classical Electrodynamics被引用 9

一句话总结

该论文揭示了磁流体力学等离子体理论中的哈密顿结构，提出了一种适用于高性能计算的无碰撞电磁磁流体力学的新哈密顿形式，并推导出首个能量与动量守恒的全F碰撞型磁流体力学模型。通过随机哈密顿方法，表明所得的碰撞算符类似于洛伦兹算符，但由于马尔可夫记忆假设，其并不精确守恒能量，提示可通过放松这些假设来改进能量守恒模型。

ABSTRACT

Hamiltonian structure is pursued and uncovered in collisional and collisionless gyrokinetic theory. A new Hamiltonian formulation of collisionless electromagnetic theory is presented that is ideally suited to implementation on modern supercomputers. The method used to uncover this structure is described in detail and applied to a number of examples, where several well-known plasma models are endowed with a Hamiltonian structure for the first time. The first energy- and momentum-conserving formulation of full-F collisional gyrokinetics is presented. In an effort to understand the theoretical underpinnings of this result at a deeper level, a \emph{stochastic} Hamiltonian modeling approach is presented and applied to pitch angle scattering. Interestingly, the collision operator produced by the Hamiltonian approach is equal to the Lorentz operator plus higher-order terms, but does not exactly conserve energy. Conversely, the classical Lorentz collision operator is provably not Hamiltonian in the stochastic sense.

研究动机与目标

• 揭示并确立无碰撞与碰撞型磁流体力学等离子体理论中的哈密顿结构。
• 开发一种适用于高性能计算的无碰撞电磁磁流体力学的新哈密顿形式。
• 基于哈密顿原理，构建首个能量与动量守恒的全F碰撞型磁流体力学模型。
• 研究传统上被认为非哈密顿的碰撞算符是否可被纳入随机哈密顿框架。
• 通过分析福克-普朗克动力学中马尔可夫假设的作用，解决随机哈密顿模型中能量守恒的不一致性。

提出的方法

• 利用变分原理和泊松括号结构，推导出无碰撞电磁磁流体力学的新哈密顿形式。
• 应用随机哈密顿建模方法处理磁矩散射，构建哈密顿福克-普朗克方程。
• 将所得碰撞算符与经典洛伦兹算符进行比较，识别偏差及能量非守恒现象。
• 通过分析马尔可夫假设在电子动力学中的作用，研究随机哈密顿模型中人为能量增长的根源。
• 提出放松马尔可夫假设，引入对过去离子位置的记忆，以实现辛几何演化并改善能量守恒。

实验结果

研究问题

• RQ1 能否系统性地在无碰撞与碰撞型磁流体力学等离子体模型中揭示哈密顿结构？
• RQ2 如何在哈密顿框架内构建一个完全能量与动量守恒的全F碰撞型磁流体力学模型？
• RQ3 经典洛伦兹碰撞算符是否与随机哈密顿力学相容？若不相容，原因是什么？
• RQ4 随机哈密顿福克-普朗克模型中人为能量增长的原因是什么，能否加以纠正？
• RQ5 能否开发一种保留记忆的随机哈密顿过程，以消除非物理解能量产生，同时保持辛结构？

主要发现

• 推导出一种适用于现代超级计算机实现的无碰撞电磁磁流体力学新哈密顿形式。
• 构建了首个能量与动量守恒的全F碰撞型磁流体力学模型，标志着碰撞型等离子体建模的重大进展。
• 随机哈密顿方法产生的碰撞算符在主导阶上与洛伦兹算符一致，但无法精确守恒能量。
• 随机哈密顿模型中的人为能量增长源于马尔可夫假设，该假设允许电子出现非物理解轨迹，使其在进入势阱后无法有效释放动能。
• 通过引入对过去离子位置的记忆，放松马尔可夫假设，可消除非物理解能量产生，同时保持辛几何演化与随机哈密顿结构。