[論文レビュー] 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

ひとこと要約

本稿では、磁場の発散を全領域でゼロに保つための新規な制約伝搬（constrained transport）手法を用いて、理想磁気流体力学（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に適用する。
• 2次元および3次元直交グリッド上でのDG法に制約伝搬（CT）フレームワークを拡張し、高次精度を維持する。
• 各要素内で局所的に発散ゼロの磁場を確保するとともに、要素界面を越えて法線磁束が連続であることを保証する。
• 磁場発散誤差に起因する数値的不安定性（物理的に不自然な負圧力や密度の出現）を回避する。
• 衝撃波付近での振動を制御するためのモーメントベースのリミッタを統合し、非線形MHD流れにおける安定性を向上させる。

提案手法

• 磁場誘導方程式を解く非連続ガラーキンスキームにより、要素中心の磁場ベクトルポテンシャルを更新する。
• ベクトルポテンシャル値を用いて、要素のエッジ（2次元）または面（3次元）における磁場の法線成分を高次精度で再構築する。
• 再構築されたエッジ／面法線磁場値を要素に戻して、全領域で発散ゼロとなる要素中心の磁場表現を生成する。
• 予測補正戦略を採用：まず標準的なDG時間積分により磁場を予測し、その後ベクトルポテンシャルに基づく再構築により全領域で∇·B = 0 を強制する。
• 衝撃波付近での保存量の不自然な振動を制御するため、モーメントベースのリミッタを用い、安定性を確保する。
• 各要素で発散がゼロであり、要素界面を越えて法線磁束が連続であることを満たすことで、エントロピー安定性を確保する。

実験結果

リサーチクエスチョン

• RQ1 制約伝搬フレームワークは、直交グリッド上での高次非連続ガラーキン法に成功裏に拡張可能か？
• RQ2 ステガレートグリッドや楕円型方程式の解法を必要とせず、高次DGフレームワーク内で全領域発散ゼロの磁場を構築する方法は何か？
• RQ3 磁場ベクトルポテンシャルが、高次精度かつ全領域発散ゼロの磁場再構築を可能にする役割は何か？
• RQ4 本手法は、衝撃波や複雑な流れ構造を含む標準的なMHDテストケースにおいて、どのように性能を発揮するか？
• RQ5 不連続性の存在下でも、振動を制御しながら高次精度を維持できるか？

主な発見

• 各要素内で局所的に発散ゼロの条件を満たし、要素界面を越えて法線磁束が連続であることを保証することで、全領域で発散ゼロの磁場を実現する。
• ベクトルポテンシャルから再構築されたエッジ／面法線磁場成分は、高次多項式近似を用いても、空間的に2次精度を達成する。
• 滑らかな領域では高次精度を維持し、衝撃波付近ではモーメントベースのリミッタにより振動を効果的に制御する。
• 標準的なMHDテストケースにおける数値結果から、発散誤差に起因する負圧力や密度の発生は観測されず、安定的かつ強固な性能が示された。
• 投影法における楕円型方程式の解法や、双曲型発散クリーニング法における自由パラメータを必要とせず、安定的かつパラメータフリーな代替手法を提供する。
• 2次元および3次元直交メッシュに適用可能であり、エッジおよび面における高次磁場成分の再構築に明示的な式が提供されている。