[论文解读] Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes

Zhiliang Xu, Dinshaw S. Balsara|arXiv (Cornell University)|Oct 4, 2011

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

一句话总结

本文提出了一种在非结构化三角形网格上求解理想磁流体动力学（MHD）方程的高阶精度、散度自由有限体积格式。该方法结合了约束传输（CT）框架、一种新型加权通量插值方法以及一种散度自由的WENO重构技术，以精确保持磁感应强度的散度自由，实现了三阶精度，并在包含激波和湍流的复杂MHD测试案例中表现出鲁棒性。

ABSTRACT

In this paper, we introduce a high-order accurate constrained transport type finite volume method to solve ideal magnetohydrodynamic equations on two-dimensional triangular meshes. A new divergence-free WENO-based reconstruction method is developed to maintain exactly divergence-free evolution of the numerical magnetic field. A new weighted flux interpolation approach is also developed to compute the z-component of the electric field at vertices of grid cells. We also present numerical examples to demonstrate the accuracy and robustness of the proposed scheme.

研究动机与目标

• 开发一种在非结构化三角形网格上求解理想MHD的高阶精度、散度自由有限体积方法，其中先前的高阶方法仅限于结构化网格。
• 解决在非结构化网格上保持离散散度自由磁感应强度的挑战，由于缺乏算术平均，单元中心磁感应强度重构并非显而易见。
• 通过采用正则性修复策略，防止在强激波和高磁压区域出现负压，从而确保模拟的鲁棒性。
• 将约束传输（CT）方法推广至非结构化三角形网格，结合WENO重构与加权通量插值，实现高阶精度。
• 在基准MHD问题（包括转子问题、爆炸波问题和Orszag-Tang涡流问题）上，展示该方法的准确性、稳定性及非振荡行为。

提出的方法

• 采用约束传输（CT）框架，将磁感应强度置于单元边上，电场置于单元顶点，以强制实现磁感应强度的离散散度自由演化。
• 提出一种新型加权通量插值方法，用于计算单元顶点处电场的z分量，确保一致性和高阶精度。
• 开发一种新型散度自由WENO重构方法，从基于边的磁感应强度值重构单元中心的磁感应强度，保持共享边上法向分量的连续性。
• 对单元中心的守恒变量（密度、动量、能量）采用标准WENO重构，以实现高阶精度。
• 在因数值误差导致出现负压的区域，应用局部正则性修复策略，基于熵更新总能量密度，以保持物理一致性。
• 在非结构化三角形网格上实现有限体积格式，结合高阶空间重构与Runge-Kutta时间积分，确保守恒性与精度。

实验结果

研究问题

• RQ1 能否在非结构化三角形网格上，通过高阶（三阶）有限体积格式，精确保持理想MHD模拟中的磁感应强度散度自由？
• RQ2 在非结构化网格的约束传输框架中，如何准确且一致地计算单元顶点处电场的z分量？
• RQ3 能否构建一种散度自由的WENO重构方法，以确保在非结构化三角形网格中共享边上法向磁感应强度分量的一致性？
• RQ4 所提出的格式在强激波和湍流等复杂MHD流动中，如何保持压力的物理正性并避免非物理振荡？
• RQ5 该方法在多大程度上可推广至三维及一般非结构化网格？

主要发现

• 所提出的二阶与三阶精度格式在光滑流动测试中达到了预期的收敛阶数，其中三阶格式在间断附近表现出更优的分辨率。
• 散度自由WENO重构确保了共享边上磁感应强度法向分量的连续性，使离散散度自由条件在机器精度范围内得以保持。
• 加权通量插值方法实现了对顶点处电场的精确且一致的计算，这对CT框架中保持散度自由性质至关重要。
• 在高磁压、低气压的爆炸波问题中，通过基于熵的局部正则性修复策略，成功防止了负压的出现，维持了物理一致性。
• Orszag-Tang涡流模拟表明，该格式能够无振荡地捕捉复杂的湍流状流动结构，并在整个计算域内保持正密度与正压强。
• 转子问题模拟显示，密度未出现退化，且压力始终保持为正，无需使用正则性修复，表明该格式在强旋转流中具有优异鲁棒性。