[论文解读] An Unsplit, Cell-Centered Godunov Method for Ideal MHD

Robert Crockett, Phillip Colella+3 more|University of North Texas Digital Library (University of North Texas)|Sep 11, 2003

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

一句话总结

本文提出了一种用于理想磁流体动力学（MHD）的二阶、无分裂、单元中心Godunov方法，通过离散投影与滤波强制满足磁场的散度自由约束。该方法在光滑流动中实现二阶精度，能正确捕捉激波，且在非无旋磁场存在时仍保持稳定，无性能退化，展现出在多维测试问题中的鲁棒性。

ABSTRACT

We present a second-order Godunov algorithm for multidimensional, ideal MHD. Our algorithm is based on the unsplit formulation of Colella (J. Comput. Phys. vol. 87, 1990), with all of the primary dependent variables centered at the same location. To properly represent the divergence-free condition of the magnetic fields, we apply a discrete projection to the intermediate values of the field at cell faces, and apply a filter to the primary dependent variables at the end of each time step. We test the method against a suite of linear and nonlinear tests to ascertain accuracy and stability of the scheme under a variety of conditions. The test suite includes rotated planar linear waves, MHD shock tube problems, low-beta flux tubes, and a magnetized rotor problem. For all of these cases, we observe that the algorithm is second-order accurate for smooth solutions, converges to the correct weak solution for problems involving shocks, and exhibits no evidence of instability or loss of accuracy due to the possible presence of non-solenoidal fields.

研究动机与目标

• 开发一种二阶、多维的Godunov格式用于理想MHD，确保在非无旋磁场存在时仍保持精度与稳定性。
• 评估不同算法组件——MAC投影、近似投影与滤波——在强制实现磁场散度自由约束方面的有效性。
• 确定是否存在一种组件组合，可在广泛的线性与非线性MHD测试问题中确保精度与稳定性。
• 将所提出的无分裂格式与现有方法（特别是八波MHD格式）进行比较，评估其在正确捕捉激波跳跃与波传播方面的性能。

提出的方法

• 该方法基于Colella的无分裂Godunov格式，采用二阶预测-校正形式，所有主变量（包括磁场）均位于单元中心。
• 通过Poisson求解器对中间面中心磁场执行MAC投影步骤，以满足到求解器容差范围内的散度自由条件。
• 使用近似投影方法，通过另一Poisson求解，对单元中心磁场施加二阶散度自由约束。
• 在每个时间步结束时应用滤波，以抑制单元中心磁场中的单极源，减少截断误差的影响。
• 通过一系列一维与二维问题对格式进行测试，包括旋转线性波、MHD激波管、低β通量管与磁化转子问题。
• 修正方程分析表明，截断误差可能导致特征向量缺陷与不稳定性；投影与滤波步骤通过平滑非无旋分量来缓解这些影响。

实验结果

研究问题

• RQ1 当通过投影与滤波强制实施散度自由约束时，无分裂、单元中心Godunov方法是否在光滑MHD流动中保持二阶精度？
• RQ2 所提出的算法是否能正确捕捉弱激波与强激波，且不会因非无旋磁场的存在而出现不稳定或精度损失？
• RQ3 在八波MHD格式中，MAC投影步骤是否对准确捕捉激波跳跃至关重要，特别是在多维设置中？
• RQ4 不同强制策略——MAC投影、近似投影与滤波——如何影响各类MHD测试问题中的精度与稳定性？
• RQ5 修正方程中截断误差强迫项在多大程度上导致不适定性或不稳定性？能否通过算法组件加以缓解？

主要发现

• 该方法在所有测试的线性与非线性问题中（包括旋转平面波与低β通量管）均实现了光滑解的二阶精度。
• 对于涉及激波的问题（如MHD激波管与磁化转子问题），该算法能正确收敛至正确的弱解。
• 在八波MHD格式中，MAC投影对准确捕捉激波跳跃至关重要；若省略该步骤，将导致激波结构错误。
• 仅使用滤波即可稳定小振幅平面波并保持二阶精度，表明其在缓解截断误差影响方面具有显著效果。
• MAC投影与滤波的组合能有效抑制非无旋分量并防止不稳定性，即使在基础格式中散度自由约束仅满足至截断误差水平时亦然。
• 该方法在多维模拟中未表现出因非无旋磁场导致的不稳定性或精度损失，证实了其鲁棒性。