[论文解读] Electrostatic gyrokinetic simulation of global tokamak boundary plasma and the generation of nonlinear intermittent turbulence

Seung Hoe Ku, R.M. Churchill+6 more|arXiv (Cornell University)|Jan 20, 2017

Magnetic confinement fusion research被引用 3

一句话总结

本论文提出了一种基于XGC1代码的新型静电陀螺动力学模拟框架，用于模拟托卡马克低β边界等离子体，结合了陀螺动力学离子与漂移动力学电子，且无需尺度分离。结果表明，非线性间歇性湍流（包括类似团块的结构）在磁分离器附近及边层区域由动力学微湍流自然产生，验证了该模型在边缘等离子体物理中的保真度。

ABSTRACT

Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-driven, low-beta electrostatic boundary plasma is reported. Gyrokinetic ions and drift-kinetic electrons are utilized without scale-separation between the neoclassical and turbulence dynamics. It is found that the nonlinear intermittent turbulence is a natural gyrokinetic phenomenon in the boundary plasma in the vicinity of the magnetic separatrix surface and in the scrape-off layer.

研究动机与目标

• 开发一种高保真度的托卡马克边界等离子体陀螺动力学模拟方法，该方法因尺度不可分、非局域、多物理场现象而具有挑战性。
• 通过在磁分离器及边层区域实现真实边界条件下的模拟，克服现有以核心为重点的陀螺动力学代码的局限性。
• 在中性流与湍流动力学之间无尺度分离的前提下，对边缘等离子体进行离子与电子的动能建模，捕捉非麦克斯韦分布行为。
• 研究在边界区域自发生成的非线性间歇性湍流（包括团块结构）的机制。
• 通过展示实验观测到的湍流特征（如团块与喷流）在边缘等离子体中的自然涌现，验证模拟框架的可靠性。

提出的方法

• 采用XGC1代码求解离子的5D静电陀螺动力学玻尔兹曼方程，同时使用漂移动力学电子以降低计算成本。
• 采用混合拉格朗日总-δf格式，结合粒子标记与相空间网格，以精确解析非麦克斯韦分布函数。
• 利用拉格朗日运动方程计算平均电场与离子陀螺中心运动，包含磁曲率与平行加速度效应。
• 通过求解准中性条件与平均电势泊松方程，确定静电势，离子密度则由陀螺中心分布推导得出。
• 模拟区域包含X点几何结构、磁分离器，并延伸至离子轴与壁面，实现全局边界条件。
• 采用亚循环动能电子方法，以在保持边层区域稳定性和精度的同时解析快速电子动力学。

实验结果

研究问题

• RQ1 是否能够通过陀螺动力学模拟框架在无需中性流与湍流动力学尺度分离的前提下，准确模拟包含磁分离器与边层的完整边界等离子体？
• RQ2 在低β、通量驱动条件下，动能电子与离子如何相互作用以生成边缘等离子体中的非线性间歇性湍流？
• RQ3 由于分离器附近微湍流的影响，类似团块的结构在模拟中在多大程度上自然涌现？
• RQ4 极向ExB流在驱动边缘区域团块运动与输运中起何种作用？
• RQ5 湍流强度与结构在边缘等离子体中如何变化，特别是在分离器附近与边层区域？

主要发现

• 非线性间歇性湍流（包括类似团块的结构）在磁分离器附近及边层区域作为包含动能电子的陀螺动力学微湍流的自然结果自发产生。
• 湍流强度在分离器处达到峰值，从ψn = 0.96开始迅速增加，此时密度梯度变得显著。
• 团块结构的极向尺寸约为7 cm，径向宽度为1.5 cm，密度扰动可达平衡密度的10–20%。
• 团块以约15 km/s的速度进行极向运动，与当地极向ExB流高度一致，表明在此构型中ExB流主导其运动。
• 团块的径向运动远小于极向运动，表明在高ExB流区域，输运主要由极向流驱动。
• 模拟在1.7 ms的模拟时间后达到边缘区域的稳态，边界层中的稳态热输运亦已实现。