[论文解读] Positivity-Preserving Finite Difference WENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations

Andrew Christlieb, Yuan Liu+2 more|arXiv (Cornell University)|Jun 19, 2014

Computational Fluid Dynamics and Aerodynamics参考文献 4被引用 4

一句话总结

本论文提出了一种用于理想磁流体动力学（MHD）的高阶正性保持有限差分WENO格式，采用约束传输（CT）方法以保持离散的无散磁场。通过将源自最大值原理保持技术的参数化通量限制器应用于高阶WENO-HCL通量，该方法在不降低CFL稳定性极限的前提下，确保了密度和压强的正性，从而在1D、2D和3D模拟中同时实现了高阶精度、无散磁场和物理解的正性。

ABSTRACT

In this paper, we utilize the maximum-principle-preserving flux limiting technique, originally designed for high order weighted essentially non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to develop a class of high order positivity-preserving finite difference WENO methods for the ideal magnetohydrodynamic (MHD) equations. Our schemes, under the constrained transport (CT) framework, can achieve high order accuracy, a discrete divergence-free condition and positivity of the numerical solution simultaneously. Numerical examples in 1D, 2D and 3D are provided to demonstrate the performance of the proposed method.

研究动机与目标

• 开发用于理想MHD的高阶有限差分WENO格式，以在低β等离子体区域下保持解的正性（密度和压强）。
• 在整个时间演化过程中保持离散的无散磁场，避免因非零散度引起的数值不稳定性。
• 在非交错框架下将正性保持与约束传输（CT）相结合，而不损害基础格式的CFL稳定性极限。
• 将原本用于标量双曲守恒律的最值原理保持通量限制技术，拓展至理想MHD方程的多组分、非线性系统。
• 在包含强激波、爆炸波和涡旋结构的1D、2D和3D基准问题中，展示所提格式的鲁棒性与精度。

提出的方法

• 该方法采用高阶WENO重构计算数值通量，并以一阶Lax-Friedrichs通量作为正性保持的构建模块。
• 在Runge-Kutta时间积分的最后阶段应用参数化通量限制器，仅修改最终通量以强制实现正性，同时保持基础格式的精度与稳定性。
• 约束传输（CT）框架通过高阶有限差分方法演化磁势，确保在离散层面磁场保持无散。
• 通过磁势演化实现散度清洁，空间导数通过WENO重构以保持高阶精度。
• 正性保持限制器源自Zhang和Shu针对可压缩Euler方程提出的框架，通过基于局部解状态限制通量，将其适配至MHD。
• 该方法在1D、2D和3D中基于结构化笛卡尔网格实现，采用出流边界条件，并对通量使用标准WENO重构。

实验结果

研究问题

• RQ1 用于理想MHD的高阶有限差分WENO格式能否在不降低CFL稳定性极限的前提下，保持密度和压强的正性？
• RQ2 在多维MHD模拟中，如何有效结合约束传输（CT）框架与高阶正性保持通量限制技术？
• RQ3 所提出的通量限制器是否在确保物理解边界的同时，保持了底层WENO格式的正式精度阶？
• RQ4 该格式在解析复杂MHD现象（如爆炸波、激波管和Alfvén波脉冲）时的性能如何，且散度误差最小？
• RQ5 在3D中，该方法能否防止低β等离子体区域中因负压强或负密度导致的数值崩溃？

主要发现

• 所提格式在所有1D、2D和3D测试案例中成功保持了密度和压强的正性，包括球形压强脉冲的3D爆炸波问题，其中未使用限制器时观察到负压强。
• 所有模拟中散度误差均保持在截断误差水平（O(h^p)），证实了约束传输框架在保持∇·B = 0方面的有效性。
• 在2D和3D爆炸问题中，相对总能量误差在整个时间过程中保持在1%以下，表明具有良好的守恒特性。
• 在光滑区域实现了高阶精度，收敛率与WENO重构的正式阶数一致（例如3阶或5阶）。
• CFL约束与底层一阶Lax-Friedrichs格式相同，意味着正性限制器未引入额外的稳定性限制。
• 在150³网格的3D爆炸问题中，密度和压强解表现出球对称性，与文献中其他高阶方法的结果一致。