[论文解读] A Singularity-Free Coordinate System for X-point Geometries

Farah Hariri, P. Hill+2 more|arXiv (Cornell University)|Jun 13, 2014

Magnetic confinement fusion research参考文献 14被引用 9

一句话总结

本论文提出了一种磁通坐标无关（FCI）的场对齐坐标系，消除了X点磁几何中的度量奇点，从而实现了托卡马克和 stellarator 中等离子体湍流的稳健数值模拟。该方法通过磁力线追踪与插值避免使用磁通坐标，在具有磁岛的平衡态中表现出快速收敛性和守恒性。

ABSTRACT

A Flux Coordinate Independent (FCI) approach for anisotropic systems, not based on magnetic flux coordinates has been introduced in [F. Hariri and M. Ottaviani, Comput. Phys. Commun., 184,2419 (2013)]. In this paper, we show that the approach can tackle magnetic configurations including X-points. Using the code FENICIA, an equilibrium with a magnetic island has been used to show the robustness of the FCI approach to cases in which a magnetic separatrix is present in the system, either by design or as a consequence of instabilities. Numerical results are in good agreement with the analytic solutions of the sound-wave propagation problem. Conservation properties are verified. Finally, the critical gain of the FCI approach in situations including the magnetic separatrix with an X-point is demonstrated by a fast convergence of the code with the numerical resolution in the direction of symmetry. The results highlighted in this paper show that the FCI approach should be able to address turbulent transport problems in X-point geometries.

研究动机与目标

• 解决磁通坐标系在处理磁构型中 X 点奇点时的根本局限性。
• 实现在包括磁岛和 X 点在内的复杂磁几何中的等离子体湍流数值模拟。
• 在分离子线和 O 点奇点区域验证 FCI 方法的鲁棒性和收敛性。
• 通过声波传播的解析解和 X 点区域的数值收敛性测试验证该方法。

提出的方法

• 开发一种磁通坐标无关（FCI）坐标系，避免在极向平面中使用磁通坐标，从而消除 X 点处的奇点。
• 通过在垂直平面之间进行磁力线追踪与插值来构建平行导数，保留各向异性特性，且不依赖于基于磁通的度量。
• 在 FENICIA 代码中实现 FCI 方法，该代码是一种专为场对齐坐标中各向异性输运设计的模块化数值方案。
• 使用 FENICIA 的平板版本模拟磁岛平衡态，测试跨分离子线的收敛性和守恒性。
• 利用椭圆积分和坐标变换，为磁岛几何中的声波传播推导解析解。
• 对磁岛的外部和内部区域应用 FCI 方法，对 X 点区域采用有界积分限的特殊处理。

实验结果

研究问题

• RQ1 FCI 方法是否能准确模拟传统磁通坐标因奇点而失效的 X 点几何中的等离子体输运？
• RQ2 FCI 方法在具有 X 点的分离子线处的收敛性和守恒性表现如何？
• RQ3 FCI 方法是否能对场对齐坐标中的 O 点（磁轴）和 X 点实现自然且稳定的处理？
• RQ4 FCI 方法在涉及磁岛的湍流输运模拟中，能在多大程度上提升数值效率和鲁棒性？

主要发现

• FCI 方法在无度量奇点的情况下成功处理了 X 点几何，如具有分离子线的磁岛平衡态所示。
• 声波传播的数值解在磁岛的外部和内部区域均快速收敛至解析解，证实了方法的准确性。
• 在对称方向上分辨率提高时，该方法表现出快速收敛，凸显其计算效率。
• 通过数值验证了系统的守恒性质，证实了 FCI 公式的鲁棒性。
• 利用第一类椭圆积分推导出波问题的解析解，并给出了 $ g(\rho) $ 在 $ \rho > 1 $ 和 $ \rho < 1 $ 区域的显式表达式。
• FCI 方法通过 $ \theta $-类坐标实现了场线标签的连续性和正则性，即使在靠近 X 点时也避免发散，方法是通过限制积分边界。