[论文解读] Continuum Kinetic Simulations of Plasma Sheaths and Instabilities

Petr Cagaš|arXiv (Cornell University)|Sep 17, 2018

Plasma Diagnostics and Applications参考文献 60被引用 7

一句话总结

本文提出了一种基于间断伽辽金法的连续介质动力学方法，用于模拟等离子体鞘层与不稳定性现象，实现了等离子体-壁相互作用的自洽建模，并基于物理的电子发射机制。该方法准确再现了朗道阻尼、两束流不稳定性及韦伯尔不稳定性，揭示了在韦伯尔不稳定性饱和阶段，电场能量可与磁场能量相媲美；同时，壁面发射模型显著改变了鞘层结构与电场分布。

ABSTRACT

A careful study of plasma-material interactions is essential to understand and improve the operation of devices where plasma contacts a wall. Key contributions of this work include (i) novel continuum kinetic algorithms with novel boundary conditions that directly discretize the Vlasov/Boltzmann equation using the discontinuous Galerkin method, (ii) fundamental studies of a plasma sheath physics with collisions, ionization, and physics-based wall emission, and (iii) theoretical and numerical studies of the linear growth and nonlinear saturation of the kinetic Weibel instability, including its role in plasma sheaths. The continuum kinetic algorithm has been shown to compare well with theoretical predictions of Landau damping of Langmuir waves and the two-stream instability. Benchmarks are also performed using the electromagnetic Weibel instability and excellent agreement is found between theory and simulation. The role of the electric field is significant during nonlinear saturation of the Weibel instability, which was not noted in previous studies. A focus is put on understanding plasma sheath physics which is essential for studying plasma-material interactions. Initial simulations are performed using a collisionless model to match classical sheath theory and the Bohm criterion. A collision operator and volumetric physics-based source terms are introduced and effects of heat flux are briefly discussed. Novel physics-based boundary conditions are developed and included in a general manner with the continuum kinetic algorithm for bounded plasma simulations. These are the first continuum kinetic simulations using self-consistent, wall emission boundary conditions with broad applicability across a variety of regimes.

研究动机与目标

• 推进聚变、推进与空间系统中至关重要的等离子体-壁相互作用高保真数值模拟。
• 开发一种直接使用间断伽辽金方法求解冯·诺依曼-玻尔兹曼方程的连续介质动力学框架，避免粒子-网格方法的局限性。
• 通过量子力学模型，实现壁面自洽、基于物理的电子发射，提升鞘层模拟的真实性。
• 研究电场在韦伯尔不稳定性非线性饱和阶段的作用，此前研究中未被充分关注。
• 建立高质量基准模拟，涵盖鞘层形成、碰撞效应与边界条件，以实现精确的等离子体-材料相互作用建模。

提出的方法

• 在相空间中对冯·诺依曼-玻尔兹曼方程实施间断伽辽金（DG）空间离散化，实现连续介质动力学模拟。
• 利用Gkeyll 2.0框架，结合机器优化、模态DG内核，实现高性能且可复现的模拟。
• 引入通用反射函数作为边界条件，实现灵活且自洽的壁面建模，包括镜面反射与吸收型反射。
• 应用BGK碰撞算子以模拟体积分层碰撞，并在鞘层区域保持分布函数的完整性。
• 以自洽方式实现基于物理的电子发射模型（Furman-Pivi与Bronold-Fehske模型），用于模拟壁面电子发射。
• 采用扫掠拟合技术从模拟中提取增长率与阻尼率，用于与线性动力学理论进行基准对比。

实验结果

研究问题

• RQ1 连续介质动力学DG方法在多大程度上能准确再现朗道阻尼与两束流不稳定性等基本动力学现象？
• RQ2 在韦伯尔不稳定性非线性饱和阶段，电场的作用是什么？其能量与磁场能量相比如何？
• RQ3 自洽的、基于物理的电子发射模型如何影响等离子体-壁系统中鞘层结构与电场分布？
• RQ4 在真实的等离子体模拟中，碰撞与热流效应在多大程度上改变鞘层动力学？
• RQ5 广义反射函数在减少计算域的同时，能否实现准确且对称的鞘层模拟，同时保持物理保真度？

主要发现

• 该连续介质动力学DG方法在两束流不稳定性中与理论增长率的偏差仅为0.3%，即使在粗略的速度分辨率下也验证了其数值精度。
• 在韦伯尔不稳定性非线性阶段，电场能量可达到与磁场能量相当的量级，这是此前未被报道的效应。
• 与理想吸收壁面相比，引入Bronold与Fehske（2015）电子吸收模型使壁面附近的电子密度提高约两倍。
• 采用Bronold与Fehske模型的边界条件时，鞘层区域电场强度降低约60%，表明其具有显著的物理影响。
• 自洽电子发射建模导致鞘层电势与离子加速特性发生可测量变化，表明在真实的等离子体-壁模拟中不能忽略电子发射。
• 该方法成功再现了无碰撞模拟中的玻姆判据与鞘层形成，离子群体在距离壁面几个德拜长度内即达到玻姆速度。