[论文解读] A Unified Gas Kinetic Scheme for Multi-scale Plasma Transport

Chang Liu, Kun Xu|arXiv (Cornell University)|Sep 17, 2016

Ionosphere and magnetosphere dynamics被引用 3

一句话总结

该论文提出了一种统一气体动理论格式（UGKS），通过将电子和离子的Vlasov-BGK方程与麦克斯韦方程结合，实现了从非碰撞到高度碰撞等离子体动力学的无缝过渡。该格式通过时间演化、尺度相关的通量重构，将粒子输运、碰撞和电磁场耦合在一起，无需分裂处理，从而实现了从无碰撞到高度碰撞等离子体动力学的统一处理。

ABSTRACT

A unified gas kinetic scheme (UGKS) for multi-scale and multi-component plasma transport is constructed. The current scheme is a direct modeling method, where the time evolution solutions from the Vlasov-BGK equations for both electron and ion, and the Maxwell equations are used to construct the scale-dependent plasma simulation. As a result, based on the modeling scales of mesh size and time step, the discretized governing equations for the whole plasma regimes are obtained. The UGKS takes into account the electron inertia, full electromagnetic field equations, and separate electron and ion evolution. The physics recovered in UGKS ranges from the kinetic Vlasov equation to the hydrodynamic magnetohydrodynamic (MHD) equations, with a unified treatment in all scales from the collisionless particle transport to the hydrodynamic wave interactions. The UGKS presents a plasma description which is more general than the Vlasov equation in the kinetic scale and all kinds of MHD equations in the hydrodynamic scale, such as Hall, Resistive, and Ideal Magnetohydrodynamics (MHD). All above single scale equations become the subsets of the UGKS. The key dynamics in UGKS is the non-splitting treatment of particle collision, acceleration, and transport in the construction of the numerical flux at cell interface, and this flux is evaluated from the scale and time evolving integral solution of the Vlasov-BGK model. The evolution of electromagnetic field and the plasma flow are coupled in a semi-implicit manner. The UGKS is able to give a physically reliable solution for plasma evolution from the collisionless limit to the highly collisional one.

研究动机与目标

• 开发一个单一的计算框架，以捕捉从动力学到流体动力学区域的所有尺度下的等离子体行为。
• 解决现有方法的局限性，这些方法需要为不同等离子体区域（如无碰撞区域使用Vlasov方程，流体行为区域使用MHD）分别建立模型。
• 在一致的数值格式中统一处理电子和离子的动力学，包括惯性和电磁场演化。
• 通过构建一个统一的格式，消除对独立求解器的需求，该格式在特定情况下可退化为已知方程（如理想MHD、Hall MHD、阻尼MHD）。
• 通过在非分裂通量公式中耦合粒子输运、碰撞和电磁场演化，确保跨尺度的物理一致性。

提出的方法

• 该格式基于对电子和离子的Vlasov-BGK方程的直接建模，利用时间演化解在单元界面构造数值通量。
• 通量通过Vlasov-BGK模型的尺度和时间依赖的积分解计算得到，确保在统一框架下准确处理粒子输运、碰撞和加速过程。
• 采用半隐式耦合策略，同时求解麦克斯韦方程与等离子体方程，以保持场-粒子一致性。
• 通过调整网格尺寸和时间步长，该方法天然捕捉从无碰撞到碰撞区域的过渡，这些参数决定了解的有效尺度。
• 通过在通量计算中同时处理粒子碰撞、加速和输运，避免了算子分裂，从而在所有尺度上保持物理精度。
• 当系统处于流体动力学极限时，该公式自然退化为已知方程（如理想MHD、阻尼MHD、Hall MHD），确保与现有模型的一致性。

实验结果

研究问题

• RQ1 是否能够通过单一数值格式在无需模型切换的情况下，准确模拟从动力学到流体动力学区域的所有尺度等离子体输运？
• RQ2 如何在多尺度下以非分裂方式一致地耦合粒子碰撞、电磁场和输运过程？
• RQ3 UGKS在何种程度上能作为已知等离子体模型（如理想MHD、Hall MHD、阻尼MHD）的极限情况被恢复？
• RQ4 该格式是否能在无碰撞和高度碰撞等离子体条件下，使用相同的算法框架保持物理精度？
• RQ5 Vlasov-BGK模型的时间演化积分解在实现尺度自适应等离子体模拟中起到什么作用？

主要发现

• UGKS提供了对等离子体动力学的统一描述，涵盖了在动力学区域的Vlasov方程，以及在流体动力学区域的所有标准MHD模型（理想MHD、阻尼MHD、Hall MHD）。
• 该格式通过自适应局部网格和时间步长尺度，成功捕捉了从无碰撞到碰撞等离子体行为的过渡。
• 对粒子输运、碰撞和电磁场演化的非分裂处理，确保了所有尺度下的物理一致性。
• 该方法在动力学极限和流体动力学极限下均保持了精度和稳定性，无需使用独立求解器或模型切换。
• UGKS在适当的渐近极限下自然退化为已知方程（如理想MHD和阻尼MHD），验证了其与现有等离子体模型的一致性。
• 电磁场与等离子体流动的半隐式耦合，使得在广泛等离子体条件下能够实现稳定且物理解释可靠的模拟。