[论文解读] Enhanced conservation properties of Vlasov codes through coupling with conservative fluid models

Thomas Alexander Trost, S. Lautenbach+1 more|arXiv (Cornell University)|Feb 1, 2017

Gas Dynamics and Kinetic Theory参考文献 21被引用 3

一句话总结

该论文提出了一种新颖的耦合方法，通过基于矩的反馈利用保守流体模型校正动能代码，从而在Vlasov模拟中增强能量守恒。通过应用广义十矩拟合过程，该方法显著减少了数值加热，同时保持了动能特征，其中十矩拟合表现出稳定且有效，而各向同性或五矩方法则不然。

ABSTRACT

Many phenomena in collisionless plasma physics require a kinetic description. The evolution of the phase space density can be modeled by means of the Vlasov equation, which has to be solved numerically in most of the relevant cases. One of the problems that often arise in such simulations is the violation of important physical conservation laws. Numerical diffusion in phase space translates into unphysical heating, which can increase the overall energy significantly, depending on the time scale and the plasma regime. In this paper, a general and straightforward way of improving conservation properties of Vlasov schemes is presented that can potentially be applied to a variety of different codes. The basic idea is to use fluid models with good conservation properties for correcting kinetic models. The higher moments that are missing in the fluid models are provided by the kinetic codes, so that both kinetic and fluid codes compensate the weaknesses of each other in a closed feedback loop.

研究动机与目标

• 解决由于相空间离散化和人工加热引起的Vlasov模拟中数值能量不守恒问题。
• 开发一种通用的、与方案无关的方法，以改善动能等离子体代码的守恒特性。
• 通过相互反馈，利用流体模型的守恒优势和动能代码的高阶矩保真度。
• 在线性和非线性等离子体区域中测试并验证基于矩的校正的有效性。
• 在多区域模拟中实现流体与动能代码之间稳健的耦合，采用保守矩拟合。

提出的方法

• 通过闭环反馈系统将Vlasov求解器与保守流体模型耦合，其中流体矩校正动能相空间分布。
• 校正过程在速度空间中使用线性变换，将动能分布拟合到预设的流体矩，将先前各向同性拟合推广至各向异性情况。
• 拟合矩阵可通过梯度下降精确计算，或近似计算以降低计算成本，两种方法均经过稳定性与准确性测试。
• 该方法在使用CWENO格式处理空间导数、时间分裂法进行积分的有限差分Vlasov求解器中实现。
• 采用十矩拟合以匹配密度、平均速度和能量密度张量，同时将五矩拟合作为对比基线进行测试。
• 模拟在GEM挑战基准上进行验证，比较再连接率、能量演化和分布函数形态。

实验结果

研究问题

• RQ1 保守流体模型能否在不引入虚假效应的情况下有效校正Vlasov模拟中的数值加热？
• RQ2 各向异性十矩拟合与各向同性或五矩拟合相比，在保持动能特征和守恒定律方面表现如何？
• RQ3 拟合矩阵的近似计算与精确计算相比，对能量守恒和模拟稳定性有何影响？
• RQ4 该耦合方法能否推广至多区域等离子体模拟中的流体-动能界面问题？
• RQ5 反馈回路是否能显著减少人工加热，同时在非线性区域保持物理准确性？

主要发现

• 采用近似矩阵计算的十矩拟合使总能量增加量相比未经校正的Vlasov模拟减少高达50%，显著改善了能量守恒。
• 各向同性拟合程序导致数值不稳定性，因此尽管是通用方法的特例，仍不适合实际应用。
• 五矩拟合破坏了关键的动能特征并引发人工不稳定性，使模拟失去物理意义。
• GEM挑战模拟表明，经校正的Vlasov代码即使在高时间分辨率下，也能保持接近参考解的再连接通量和电流分布。
• 该方法在不同分辨率和时间步长下均表现稳健，对CWENO格式中正则化参数ε的选择不敏感。
• 反馈回路通过校正相空间扩散成功稳定了模拟，同时不损害等离子体的动能描述。