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

ひとこと要約

本稿では、離散的投影とフィルタリングを通じて磁場の発散自由制約を満たす、2次精度で、分解けない、セル中心のゴドノフ法を用いた理想磁気流体力学（MHD）の手法を提示する。この手法は滑らかな流れに対して2次精度を達成し、衝撃波を正しく捉え、非ソリノイダルな磁場によって劣化することなく安定性を保つ。多次元のテスト問題においても頑健性を示している。

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.

研究の動機と目的

• 非ソリノイダルな磁場が存在する状況下でも、精度と安定性を維持する2次精度の多次元ゴドノフスキームを、理想MHDに開発すること。
• MAC投影、近似投影、およびフィルタリングといった異なるアルゴリズム的要素が、磁場の発散自由制約をどのように強制するかを評価すること。
• これらの要素の1つの組み合わせが、線形および非線形MHDテスト問題の広範なスイートにおいて、精度と安定性を保証できるかを特定すること。
• 提案された分解けないスキームの性能を、特に8波MHD定式化と比較し、正しい衝撃波ジャンプと波の伝播を捉える能力を評価すること。

提案手法

• すべての主要変数（磁場を含む）をセル中心に配置した、コールラの分解けないゴドノフスキームに基づく2次精度の予測・修正形式を用いる。
• 発散自由条件を解法の許容誤差まで満たすために、ポissonソルバを用いて中間の面中心磁場にMAC投影ステップを適用する。
• セル中心磁場に対する2次精度の発散自由制約を、別のポissonソルバを用いて近似投影により強制する。
• 各タイムステップの終了時にフィルタを適用し、セル中心磁場内のモノポール源を抑制することで、切り捨て誤差の影響を緩和する。
• 1次元および2次元の問題のスイートを用いてスキームをテストする。これには、回転した線形波、MHD衝撃波管、低βフラックスチューブ、磁場付きローター問題が含まれる。
• 修正方程式解析により、切り捨て誤差が固有ベクトルの欠落や不安定性を引き起こす可能性があることが示された。投影とフィルタリングのステップにより、非ソリノイダル成分のスムージングが行われ、これらの影響が緩和される。

実験結果

リサーチクエスチョン

• RQ1 発散自由制約が投影とフィルタリングによって強制される場合、分解けないセル中心ゴドノフスキームは、滑らかなMHD流れにおいて2次精度を維持するか？
• RQ2 非ソリノイダルな磁場によって引き起こされる不安定性や精度の損失が生じない状況で、提案されたアルゴリズムは弱いおよび強いMHD衝撃波を正しく捉えることができるか？
• RQ3 8波MHD定式化において、特に多次元設定下で、MAC投影ステップは正確な衝撃波ジャンプを実現するために不可欠であるか？
• RQ4 MAC投影、近似投影、およびフィルタリングといった異なる強制戦略は、多様なMHDテスト問題において精度と安定性にどのように影響を与えるか？
• RQ5 修正方程式に現れる切り捨て誤差駆動項が、不適切な問題設定や不安定性を引き起こす程度はどの程度であり、アルゴリズム的要素によってこれを緩和できるか？

主な発見

• 回転した平面波や低βフラックスチューブを含む、すべてのテストされた線形および非線形問題において、滑らかな解に対して2次精度が達成された。
• MHD衝撃波管や磁場付きローター問題といった衝撃波を含む問題において、正しい弱解に収束することが確認された。
• 8波MHD定式化において、MAC投影は正確な衝撃波ジャンプを実現するために不可欠であり、その欠如は誤った衝撃構造を引き起こす。
• フィルタのみで、小振幅平面波が安定化され、2次精度が維持されることから、切り捨て誤差の影響を緩和する有効性が示された。
• MAC投影とフィルタの組み合わせにより、非ソリノイダル成分が効果的に抑制され、ベーススキームで発散自由制約が切り捨て誤差の範囲でのみ満たされる場合でさえも、不安定性が防止された。
• 非ソリノイダルな磁場による不安定性や精度の損失の兆候は一切認められず、多次元シミュレーションにおける頑健性が確認された。