[论文解读] A Dispersion Minimized Mimetic Method for Cold Plasma

Vrushali A. Bokil, Vitaliy Gyrya+1 more|arXiv (Cornell University)|Apr 5, 2016

Electromagnetic Simulation and Numerical Methods参考文献 1被引用 3

一句话总结

本文提出了一种用于冷等离子体中麦克斯韦方程组的最小化色散性拟谱有限差分（MFD）方法，结合指数时间差分（ETD）与一种新型的 m-自适应优化技术。通过在空间上应用 m-自适应于广义质量凝聚 MFD 格式，并在时间上采用 ETD，该方法实现了四阶数值色散误差——相比标准二阶格式，误差降低了两个数量级——该结果通过数值与理论分析均得到验证。

ABSTRACT

In this paper we consider the lowest edge-based mimetic finite difference (MFD) discretization in space for Maxwell's equations in cold plasma on rectangular meshes. The method uses a generalized form of mass lumping that, on one hand, eliminates a need for linear solves at every iteration while, on the other hand, retains a set of free parameters of the MFD discretization. We perform an optimization procedure, called m-adaptation, that identified a set of free parameters that lead to the smallest numerical dispersion. The choice of the time stepping proved to be critical for successful optimization. Using exponential time differencing we were able to reduce the numerical dispersion error from second to fourth order of accuracy in mesh size. It was not possible to achieve this order of magnitude reduction in numerical dispersion error using the standard leapfrog time stepping. Numerical simulations independently verify our theoretical findings.

研究动机与目标

• 为减少在冷等离子体介质中电磁波模拟的数值色散误差，这对长时间、电大尺寸区域的模拟至关重要。
• 将此前在真空中取得成功的 m-自适应技术扩展至线性色散介质（如冷等离子体），由于引入了额外的动力学行为，该问题在冷等离子体中具有非平凡性。
• 在拟谱有限差分格式族中识别最优自由参数，以最小化数值色散，同时保持显式时间推进。
• 证明指数时间差分（ETD）能够实现更高阶色散抑制，而标准跳跃蛙时间推进则不能。
• 开发一种模板尺寸小于现有四阶方法的方法，从而提高计算效率。

提出的方法

• 在矩形网格上，采用最低阶边值型拟谱有限差分（MFD）方法对冷等离子体中的麦克斯韦方程组进行空间离散化。
• 应用广义质量凝聚形式，以消除每时间步所需的线性求解，同时保留自由参数用于优化。
• m-自适应过程通过傅里叶域中的色散关系分析，优化这些自由参数以最小化数值色散。
• 采用指数时间差分（ETD）进行时间积分，实现对刚性线性项的精确积分并获得更高阶精度。
• 由此产生的指数时间拟谱有限差分（ETMFD）格式实现了四阶数值色散误差，而标准 ETMFD 格式的基准误差为二阶。
• 该方法保持了 E 场与 J 场相对于 H 场的显式交错结构，避免了半隐式时间平均。

实验结果

研究问题

• RQ1 m-自适应能否成功地从真空扩展到包含极化电流辅助微分方程的冷等离子体模型？
• RQ2 为何在冷等离子体的 MFD 家族中，跳跃蛙时间推进无法将数值色散抑制提升至二阶以上，而 ETD 却可以？
• RQ3 在广义质量凝聚 MFD 格式中，使冷等离子体中数值色散最小化的最优自由参数集是什么？
• RQ4 ETMFD 方法是否实现了四阶数值色散误差？该结果能否通过独立模拟验证？
• RQ5 优化后的 ETMFD 方法的模板尺寸与文献中其他四阶方法相比如何？

主要发现

• m-自适应过程在结合指数时间差分（ETD）时，成功地将数值色散误差从二阶降低至四阶（关于网格尺寸）。
• ETMFD 方法在数值色散误差上实现了四阶收敛，该结论通过理论分析与数值模拟均得到证实。
• ETMFD 格式中的最优自由参数与局部质量矩阵与真空中情况完全相同，表明不同模型间具有结构一致性。
• ETMFD 方法保持了紧凑的模板尺寸，相较于依赖更高阶基础离散化的其他四阶方法，具有计算优势。
• 标准跳跃蛙时间推进在冷等离子体模型中无法通过 m-自适应实现更高阶色散抑制，凸显了时间积分方法选择的关键作用。
• 对特解的数值模拟显示，电场与电流密度的 L² 误差均呈现四阶收敛，验证了理论结果。