Development of a Quasi-Static MHD Solver for Turbulent Flows with Non-Uniform Applied Magnetic Fields in OpenFOAM

技術要約：非一様な印加磁場を伴う乱流のための準静的MHDソルバーの開発（OpenFOAMにおける）

問題提起
標準的なOpenFOAMの磁気流体力学（MHD）ソルバーは、通常、ローレンツ力をマクスウェル応力テンソルを介して評価する磁気誘導定式化を利用している。この手法は、一様な印加磁場に対しては堅牢であるが、印加磁場（$B_0$）が非一様である場合には深刻な離散的キャンセル誤差に苦しむ。連続体極限において、マクスウェル応力の発散と磁気圧の勾配は正確に相殺される。しかし、空間的に変化する場を伴う離散的な有限体積設定では、この相殺は不完全である。その残差は、$B_0^2$に比例するスプリアス（偽の）体積力として作用する。現実的な工学的構成（例：0.15 mのダクト内における10 Tの磁場）では、この残差は$\sim 10^3$ m s$^{-2}$もの加速度に達することがあり、これは物理的なローレンツ加速度よりも2桁大きい。この誤差は、メッシュの細分化やタイムステップの削減によっても軽減されないため、数値的な解像度の問題ではなく、定式化自体の根本的な欠陥であることを示している。

手法
これに対処するため、著者は準静的（低磁気レイノルズ数、$R_m \ll 1$）電位定式化に基づく有限体積ソルバーであるmhdturbFoamQSを開発した。コアとなる手法は以下の通りである：

1. 定式化の転換： マクスウェル応力テンソルからローレンツ力を計算する代わりに、電荷保存則（$\nabla \cdot j = 0$）から導かれる電位（$\phi$）に関するポアソン方程式を解く。電流密度（$j$）はオームの法則（$j = \sigma(-\nabla\phi + u \times B_0)$）から代数的に再構成され、ローレンツ力は直接 $j \times B_0$ として評価される。これにより、大きな$O(B_0^2)$項の減算を回避している。
2. 一般化された場の取り扱い： 印加磁場 $B_0$ は定数ではなく、一般的な volVectorField として扱われ、大きさや方向の空間的な変化を許容する。
3. 整合的な境界条件： 離散的なポアソン方程式の両側に、絶縁壁条件（$j \cdot n = 0$）の特定のインプリメンテーションが適用される。これにより、運動起電力のフラックスが右辺から誤って除去されることを防ぎ、意図しない電極条件の賦課を防止している。
4. 暗黙的な磁気減衰： ローレンツ力の分解から生じる硬い線形減衰項（$-\sigma|B_0|^2 u$）は、運動量行列内で fvm::Sp を用いて暗黙的に処理される。これにより、磁気減衰に伴う制限的なタイムステップ制約（本来は$1/B_0^2$に比例する）を取り除いている。
5. 乱流結合： 本ソルバーは、輸送方程式を変更することなく、標準的なOpenFOAMの乱流モデル（Spalart–Allmaras, $k$–$\epsilon$, $k$–$\omega$ SST）と統合されている。これにより、電磁気的な定式化を乱流モデルの変化から隔離している。

主要な貢献

• 特定と定量化： 本論文は、非一様な場における誘導ベースのソルバーにおける離散的キャンセル誤差を特定・定量化し、その誤差がメッシュの細分化に伴って減衰するのではなく、飽和することを実証した。
• ソルバーの開発： 準静的電位ルート、整合的な絶縁壁処理、および暗黙的な磁気シンク処理を組み込んだ mhdturbFoamQS（344行のC++コード）の実装。
• 検証フレームワーク： 以下の包括的な検証スイートを含む：
  • 平面チャネルの厳密解（Hartmann解）との比較（摩擦係数誤差 < 0.17%、速度プロファイル誤差 < 0.5%）。
  • 絶縁壁を持つ円管におけるShercliff–Gold解との比較。
  • 電流が存在しない非一様な印加磁場を用い、流体が初期状態で静止している状態での「ヌルテスト」。準静的ソルバーは正しいゼロ速度を返すのに対し、誘導ソルバーはスプリアスな運動を生成する。
  • ポテンシャルおよび再構成された電流の両方について、2次精度を確認する製造解（Manufactured solution）テスト。
• 再現性： ソルバー、検証ケース、および後処理スクリプトは、Zenodo上のconcept DOIの下でGPL v3として公開されており、報告されたすべての数値が独立して再現可能であることを保証している。

結果

• 検証： 本ソルバーは、層流の解析解を高精度で再現する。ヌルテストにおいて、誘導ソルバーは流体が静止しているはずの場において最大0.36 m s$^{-1}$（厳密な壁の値を用いた場合は3.7 m s$^{-1}$）のスプリアス速度を生成したが、準静的ソルバーはマシン精度でのゼロを維持した。
• 乱流の実証： 本ソルバーを、非一様な交互磁石場（6回対称の方位角構造）を伴う乱流ダクト流（$Re \approx 10^4$）に適用した。シミュレーションは、磁場強度最大10 Tにおいて、3つの異なる乱流モデルを用いて実行された。
  • ソルバーは安定しており、誘導ベースのアプローチで見られる発散を示さなかった。
  • 結果は、体積平均圧力および乱流運動エネルギーが磁場強度とともに単調に増加することを示した。
  • 乱流モデルの選択に対する顕著な感度が観察された（例：10 Tにおいて、圧力上昇は $k$–$\epsilon$ の+17.3%から Spalart–Allmaras の+30.1%まで変動した）。これは、電磁気的な結合が乱流モデリングの不確実性と切り離されていることを浮き彫りにしている。

意義と主張
本論文は、準静的電位定式化が、空間的に変化する磁場を伴う乱流MHDフローをシミュレートするための実用的かつ堅牢な手法を提供し、マクスウェル応力定式化に固有の離散的キャンセル病理を効果的に排除することを主張している。本研究は、電磁気的な結合と数値的定式化を評価することに重点を置いており、乱流物理を新たに提案するのではなく、標準的なモデルを修正せずに使用することで、数値的影響を分離している。本研究は、低$R_m$近似が成立する限り、既存の誘導ベースの実装を悩ませるスプリアス体積力なしに、工学規模の非一様な磁場を扱えることを確立している。著者は、制限事項として、乱流モデルにおける明示的な磁気項の欠如、絶縁壁への限定、および時間ステップ分割における1次の分裂誤差を挙げており、これらを今後の課題と定義している。