[论文解读] High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD

James A. Rossmanith|arXiv (Cornell University)|Oct 16, 2013

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

一句话总结

本文提出了一种用于理想磁流体动力学（MHD）的高阶不连续伽辽金有限元方法，通过一种新型约束传输方法实现全局无散磁场。通过引入基于DG格式更新的单元中心磁矢势，并以高阶精度重构边/面法向磁场分量，该方法确保了局部无散约束与法向磁通量的全局连续性，实现了无需椭圆求解或自由参数的稳定、高阶精度模拟。

ABSTRACT

The modification of the celebrated Yee scheme from Maxwell equations to magnetohydrodynamics is often referred to as the constrained transport approach. Constrained transport can be viewed as a sort of predictor-corrector method for updating the magnetic field, where a magnetic field value is first predicted by a method that does not preserve the divergence-free condition on the magnetic field, followed by a correction step that aims to control these divergence errors. This strategy has been successfully used in conjunction with a variety of shock-capturing methods including WENO, central, and wave propagation schemes. In this work we show how to extend the basic CT framework to the discontinuous Galerkin finite element method on both 2D and 3D Cartesian grids. We first review the entropy-stability theory for semi-discrete DG discretizations of ideal MHD, which rigorously establishes the need for a magnetic field that satisfies the following conditions: (1) the divergence of the magnetic field is zero on each element, and (2) the normal components of the magnetic field are continuous across element edges/faces. In order to achieve such a globally divergence-free magnetic field, we introduce a novel CT scheme that is based on two ingredients: (1) we introduce an element-centered magnetic vector potential that is updated via a discontinuous Galerkin scheme on the induction equation; and (2) we define a mapping that takes element-centered magnetic field values and element-centered magnetic vector potential values and creates on each edge /face a representation of the normal component of the magnetic field; this representation is then mapped back to the elements to create a globally divergence-free element-centered representation of the magnetic field. For problems with shock waves, we make use of so-called moment-based limiters to control oscillations in the conserved quantities.

研究动机与目标

• 开发一种用于理想MHD的高阶不连续伽辽金有限元方法，严格强制磁感应强度满足无散条件。
• 将约束传输（CT）框架扩展至二维和三维结构化笛卡尔网格上的DG方法，同时保持高阶精度。
• 确保每个单元内部磁感应强度局部无散，且单元界面处法向磁感应强度分量连续。
• 避免因磁感应强度散度误差导致的数值不稳定，此类误差可能引发非物理的负压强或负密度。
• 集成基于矩量的限制器以控制激波附近振荡，提升非线性MHD流中模拟的鲁棒性。

提出的方法

• 引入基于DG格式求解感应方程更新的单元中心磁矢势。
• 利用磁矢势值构建单元边（二维）或面（三维）上法向磁感应强度分量的高阶表示。
• 将重构的边/面法向磁感应强度值映射回单元，生成全局无散、单元中心的磁感应强度表示。
• 采用预测-校正策略：首先通过标准DG演化预测磁感应强度，然后通过基于磁矢势的重构进行校正，以全局强制满足∇·B = 0。
• 使用基于矩量的限制器控制激波附近守恒量的虚假振荡，维持稳定性。
• 通过满足必要条件（每个单元内散度为零，界面处法向磁通量连续）确保熵稳定性。

实验结果

研究问题

• RQ1 约束传输框架能否成功扩展至笛卡尔网格上理想MHD的高阶不连续伽辽金方法？
• RQ2 在不使用非结构化网格或椭圆求解的情况下，如何在高阶DG框架中构建全局无散磁感应强度？
• RQ3 磁矢势在实现高阶、全局无散磁感应强度重构中起到何种作用？
• RQ4 所提方法在涉及激波和复杂流场结构的标准MHD测试算例上的表现如何？
• RQ5 该方法能否在保持间断附近振荡控制的同时维持高阶精度？

主要发现

• 通过强制每个单元内局部无散条件及单元界面处法向磁感应强度分量连续，实现了全局无散磁感应强度。
• 从磁矢势重构边/面法向磁感应强度分量，即使使用高阶多项式近似，也能保证法向通量的二阶空间精度。
• 在光滑区域保持高阶精度，并通过基于矩量的限制器有效控制激波附近的振荡。
• 在标准MHD测试算例上的数值结果表明方法具有鲁棒性和稳定性，未因散度误差观察到负压强或负密度。
• 该方法避免了投影方法所需的椭圆求解或超分方法中的自由参数，提供了一种稳定、无自由参数的替代方案。
• 该框架适用于二维和三维结构化笛卡尔网格，且提供了在边和面上重构高阶磁感应强度分量的显式公式。