磁気流体力学（MHD）
磁気流体力学（MHD）の理論基礎
磁気流体力学の基礎
先生、磁気流体力学（MHD）って、流体力学と電磁気学が合体したものですか？
その通り。MHDは導電性流体（液体金属、プラズマ、電解質溶液など）中での電磁場と流体運動の相互作用を扱う分野だ。電流が流れる流体に磁場が作用するとLorentz力が生じ、流れが変化する。逆に、導電性流体の運動は磁場を誘導する。
身近な応用例:
- 連続鋳造: 溶鋼の電磁ブレーキ、電磁攪拌
- アルミ電解精錬: Hall-Heraultセル内の溶融アルミの流動制御
- 核融合炉: プラズマの磁場閉じ込め
- MHDポンプ: 可動部なしで液体金属を搬送
- 宇宙推進: MPDスラスタ
MHDの支配方程式
MHDはNavier-Stokes方程式とMaxwell方程式を連成させた体系だ。
修正Navier-Stokes方程式（Lorentz力の追加）:
ここで $\mathbf{J} \times \mathbf{B}$ がLorentz力（体積力）。
オームの法則（運動する導電性流体中）:
$\sigma$ は電気伝導率、$\mathbf{E}$ は電場、$\mathbf{u} \times \mathbf{B}$ は流体運動による起電力。
磁場の誘導方程式:
右辺第1項が磁場の対流（凍結）、第2項が磁場の拡散。
流体の運動方程式に電磁力が加わり、同時に磁場の方程式にも流速が入る。完全な双方向連成ですね。
MHDの無次元数
MHD特有の無次元数を整理しよう。
| 無次元数 | 定義 | 物理的意味 |
|---|---|---|
| Hartmann数 Ha | $BL\sqrt{\sigma/\mu}$ | 電磁力/粘性力 |
| 磁気Reynolds数 Rm | $\mu_0 \sigma U L$ | 磁場対流/磁場拡散 |
| Stuart数（相互作用パラメータ）N | $\sigma B^2 L / (\rho U) = \text{Ha}^2/\text{Re}$ | 電磁力/慣性力 |
| 磁気Prandtl数 Pm | $\mu_0 \sigma \nu$ | 運動量拡散/磁場拡散 |
工業的な液体金属（溶鋼、溶融アルミなど）では $\text{Pm} \sim 10^{-6}$ と極めて小さい。つまり磁場の拡散が流体の運動量拡散よりはるかに速く、$\text{Rm} \ll 1$ となる。この場合、磁場は外部印加場でほぼ決まり、流体への影響だけを考えればよい（低Rm近似）。
Hartmann数が大きいほど磁場の影響が強い、ということですね。
そうだ。Ha = 0 で通常の流体力学、Ha → ∞ で流れが完全に磁場に拘束される。連続鋳造では Ha ∼ 100-1000 の範囲が典型的だ。
核融合炉のプラズマ閉じ込め——MHDが人類の夢を守っている
ITER（国際核融合実験炉）では、1億℃のプラズマをトカマク型の磁場で閉じ込めます。このプラズマの安定性解析の根幹がMHD理論です。プラズマが磁場から「脱走」しないか、「キンク不安定」や「バルーニング不安定」が発生しないか——これらをMHD方程式の固有値問題として解くことで、コイル設計や電流プロファイルを最適化します。核融合が「30年後に実現」と言われ続けているのは、このMHD安定性解析の難しさを物語っています。
磁気流体力学（MHD）の数値計算手法
MHDの数値解法
MHDの連成方程式をCFDでどうやって解くんですか？
低Rm近似（工業液体金属の大半）と高Rm問題（プラズマ、天体物理）で手法が異なる。
低Rm近似の解法
$\text{Rm} \ll 1$ の場合、誘導磁場は印加磁場に比べて無視でき、磁場は既知の外部場 $\mathbf{B}_0$ として扱える。電場の支配方程式は:
$\phi$ は電位。この Poisson方程式を解いて $\mathbf{E} = -\nabla\phi$ を求め、$\mathbf{J} = \sigma(-\nabla\phi + \mathbf{u} \times \mathbf{B}_0)$、$\mathbf{F}_L = \mathbf{J} \times \mathbf{B}_0$ を計算してN-S方程式のソース項に加える。
計算手順:
1. 流れ場から $\mathbf{u} \times \mathbf{B}_0$ を計算
2. 電位のPoisson方程式を解く
3. 電流密度 $\mathbf{J}$ を計算
4. Lorentz力 $\mathbf{F}_L$ をN-S方程式に加えて流れ場を更新
5. 収束まで1-4を繰り返す
電位の方程式が1つ増えるだけで、基本的にはN-Sの延長線上にあるんですね。
高Rm問題の解法
$\text{Rm} \gg 1$（プラズマ物理、天体物理）では磁場の誘導方程式をフルに解く必要がある。
主要な数値手法:
- Constrained Transport (CT) 法: $\nabla \cdot \mathbf{B} = 0$ をセルフェースでの磁束保存で厳密に維持
- Divergence Cleaning法: $\nabla \cdot \mathbf{B}$ の誤差をダンピングする補正方程式を追加
- ベクトルポテンシャル法: $\mathbf{B} = \nabla \times \mathbf{A}$ として $\nabla \cdot \mathbf{B} = 0$ を自動的に保証
$\nabla \cdot \mathbf{B} = 0$ の維持は数値MHDの根本的な課題だ。これが崩れると非物理的な磁気モノポール力が発生し、計算が破綻する。
Hartmann流れ — MHDの基本検証問題
MHDの最も基本的な解析解がHartmann流れだ。平行平板間の流れに壁面に垂直な磁場 $B_0$ を印加する。
Ha = 0 で放物線分布（Poiseuille流れ）、Ha → ∞ で壁面近傍のHartmann層（厚さ $\delta_H \sim H/\text{Ha}$）を除いてフラットな分布になる。
| Ha | 速度分布 | Hartmann層厚さ |
|---|---|---|
| 0 | 放物線 | なし |
| 10 | やや平坦化 | $H/10$ |
| 100 | 中心部ほぼフラット | $H/100$ |
| 1000 | Hartmann層のみに勾配 | $H/1000$ |
Ha = 1000 だとHartmann層がチャネル幅の1/1000しかないんですね。メッシュの要求が厳しそう。
その通り。Hartmann層を解像するには壁面近傍にHa数に比例した数のメッシュ層が必要になる。これがMHD計算の主要なコストの一つだ。
MHD数値解法の難関——「磁場の発散ゼロ」を保つ苦労
MHD解析の数値実装で最も悩ましいのが ∇·B=0（磁束の連続性）の維持です。有限体積法で磁場を離散化すると、数値誤差が蓄積して ∇·B≠0 になってしまい、非物理的な磁力が発生して解が歪む。これを防ぐ手法として「発散クリーニング（divergence cleaning）」や「制約輸送法（constrained transport）」が開発されました。特に太陽風の大規模MHDシミュレーションでは、発散誤差が1%を超えると数値的な「磁気モノポール」が生まれ、プラズマが誤った方向に飛んでいく。磁場の発散ゼロは物理の要請であり、数値解法の戦場です。
磁気流体力学（MHD）の実務適用
MHD解析の実務フロー
実際にMHDのCFD解析を行う手順を教えてください。
連続鋳造の電磁ブレーキ（EMBr）を例に説明しよう。
Step 1: 磁場の計算
- まず電磁場解析ソフト（Ansys Maxwell, COMSOL AC/DC, JMAGなど）で磁石・コイルの形状から磁場分布 $\mathbf{B}_0(\mathbf{x})$ を計算
- 結果をCFDソルバーにインポート（マッピング）
Step 2: 流体ドメインのメッシュ
- Hartmann層の解像: 壁面第一層厚さ ≈ $H/\text{Ha}$ の1/5以下
- 典型例: $H = 0.1\,\text{m}$, $\text{Ha} = 500$ → $\delta_H = 0.2\,\text{mm}$ → 第一層 = 0.04 mm
- 壁面からの成長率1.1-1.2でプリズム層を配置
Step 3: 物性値設定
| 物性値 | 溶鋼 (1550°C) | 溶融アルミ (700°C) |
|---|---|---|
| 密度 $\rho$ | 7000 kg/m$^3$ | 2385 kg/m$^3$ |
| 粘度 $\mu$ | 0.006 Pa·s | 0.0013 Pa·s |
| 電気伝導率 $\sigma$ | $7.14 \times 10^5$ S/m | $3.5 \times 10^6$ S/m |
| 熱伝導率 $k$ | 40 W/(m·K) | 100 W/(m·K) |
Step 4: 境界条件
- 入口: ノズルからの流入速度（通常 1-3 m/s）
- 出口: 圧力指定
- 壁面: no-slip + 電気的絶縁壁 or 導電壁
- 電気的境界条件が重要: $\mathbf{J} \cdot \mathbf{n} = 0$（絶縁壁）or $\phi = 0$（完全導電壁）
壁面の電気的境界条件がMHD特有ですね。導電壁と絶縁壁で結果が変わりますか？
大きく変わる。導電壁では壁面を通じて電流が帰還するため、Lorentz力のブレーキ効果が増大する。Hartmann層の厚さも異なる。絶縁壁では $\delta_H \sim H/\text{Ha}$、導電壁では $\delta_H \sim H/\sqrt{\text{Ha}}$ となる。
検証のポイント
MHD解析の検証チェックリスト:
- Hartmann流れの再現: 解析解との比較。速度分布と圧力損失が一致するか
- 圧力損失の検証: MHDによる追加圧力損失: $\Delta p_{\text{MHD}} \sim \sigma U B^2 L$（Ha²に比例）
- 電流の保存: $\nabla \cdot \mathbf{J} = 0$ がどの程度満たされているか確認
- 磁場のダイバージェンスフリー: $\nabla \cdot \mathbf{B}$ の最大値をモニター
- 実験との比較: UDV（Ultrasonic Doppler Velocimetry）による液体金属の速度計測データと比較
液体金属は不透明だからPIVが使えないんですね。超音波を使うのか。
その通り。液体金属の流速計測にはUDVやCIFT（Contactless Inductive Flow Tomography）が使われる。実験データは限られているので、解析解での検証が特に重要だ。
電磁攪拌（EMS）による連続鋳造品質改善——鉄鋼業のMHD応用
鉄鋼業の連続鋳造プロセスでは「電磁攪拌（EMS: Electromagnetic Stirring）」が品質改善の標準技術として普及している。鋳型や二次冷却帯に設置したコイルで回転磁場を生成し、溶鋼にローレンツ力による攪拌流れを誘起することで、①中心偏析の抑制、②介在物の浮上促進、③等軸晶率の向上を実現する。CFD-MHD連成解析では溶鋼の流動、温度場、凝固フロントを同時にシミュレーションし、EMS電流・周波数・配置の最適化に活用される。国内製鉄所の実証では、EMS条件をCFDで最適化することで中心偏析指数が0.05以上改善（0→Cが0.20→0.15へ）し、高強度鋼の歩留まりが5%向上した事例が鉄鋼学会誌で報告されている。
磁気流体力学（MHD）のソフトウェア比較
MHD対応ツールの比較
MHDのシミュレーションに使えるソフトを教えてください。
MHDは電磁場とCFDの連成が必要なので、ツール選定が重要だ。
Ansys Fluent（MHDモジュール）
FluentにはMHDモジュールが搭載されている。
- 対応: 低Rm近似（電位法）。外部磁場を読み込み、Lorentz力を自動計算
- 設定: Models > MHD で有効化。磁場はUser-Defined or Profile入力
- 連成: Ansys MaxwellやAnsys EMFからの磁場データをWorkbench経由でマッピング
- 乱流: MHDが乱流に与える影響のモデリング（Sommeria-Moreau, Knaepen-Moreau）
- 制限: 高Rm問題、自己誘導磁場は非対応
COMSOL Multiphysics
COMSOLはMHDに非常に強い。
- AC/DCモジュール + CFDモジュール の連成でMHDを解く
- 低Rm/高Rm両対応: 電位法とフル誘導方程式の両方をサポート
- Lorentz力の自動計算: 電磁場と流体場の双方向連成が標準機能
- AC磁場対応: 交流磁場による電磁攪拌のモデリング
- 弱点: 大規模3D問題（1000万セル超）ではメモリ・計算時間が課題
OpenFOAM
OpenFOAMにはMHDソルバーが含まれている。
- mhdFoam: 磁場の誘導方程式をフルに解くソルバー（高Rm対応）
- electrokinetics: 電場駆動の流れ（EOF、電気化学）
- カスタム実装: 低Rm近似の電位法ソルバーを自作する研究者が多い
- コミュニティ: Hartmann流れ、MHD乱流のチュートリアルが多数
専門ツール
| ツール | 用途 | 特徴 |
|---|---|---|
| ANSYS Maxwell | 電磁場計算 | 磁場分布の計算。CFDとの連成用 |
| JMAG | 電磁場計算 | 日本製。モーター・変圧器設計に強い |
| Plato/Ath ena | MHD専用（天体物理） | Godunov法ベースの高Rmソルバー |
| FLASH | 天体MHD | AMR対応、核融合・天体物理 |
| SFEMaNS | MHD専用 | スペクトル有限要素法。精度が高い |
連成ワークフロー
電磁場とCFDを別々のソフトで計算する場合のワークフローはどうなりますか？
典型的な連成ワークフロー（連続鋳造の場合）:
```
Maxwell/JMAG → 磁場分布B(x,y,z) → エクスポート
↓
Fluent/STAR-CCM+ で
CFDメッシュにマッピング
↓
MHD連成計算（定常 or 非定常）
↓
```
磁場のマッピング精度が重要だ。電磁場の計算メッシュとCFDメッシュは通常異なるため、補間誤差に注意が必要。Ansys Workbenchの場合、System Coupling でメッシュ間の自動マッピングが可能。
選定ガイドライン
| 用途 | 推奨 | 理由 |
|---|---|---|
| 連続鋳造EMBr | Fluent MHD + Maxwell | 実績が豊富 |
| アルミ電解槽 | COMSOL | AC磁場対応 |
| 核融合ブランケット | OpenFOAM | 高Ha数対応のカスタム実装 |
| 天体物理 | Athena++, FLASH | 高Rm、圧縮性MHD |
| 研究（汎用） | COMSOL | 学習コストが低い |
工業用途ならFluent+Maxwell、学術研究ならCOMSOLかOpenFOAMという使い分けですね。
MHDソルバー——オープンソースが先行する珍しい分野
MHD解析は商用ツールよりオープンソースが強い、数少ない流体分野の一つです。OpenFOAMのmhdFoamやGPU対応のAthena++は、核融合・宇宙物理コミュニティが長年開発を主導しており、商用コードより最新アルゴリズムの実装が早い。一方、ANSYS FluentのMHDモジュールは産業用途（電磁攪拌、アルミ鋳造）に特化した検証データが充実しています。「研究はオープンソース、産業はANSYS」という棲み分けが明確なのがMHD分野の特徴です。
磁気流体力学（MHD）の先端研究
MHD乱流
磁場が乱流に与える影響はどのようなものですか？
MHD乱流は通常の乱流とは質的に異なる振る舞いを示す。外部磁場は乱流の異方性を引き起こす。
磁場方向のスケールが大きく、磁場に垂直な方向のスケールが小さい「シガー型」の渦構造が形成される。これはJoule散逸が磁場に垂直な速度成分を減衰させるためだ。
MHD乱流のエネルギースペクトルは、通常のKolmogorov理論（$E(k) \propto k^{-5/3}$）から逸脱する。
- 弱MHD（$N \ll 1$）: Kolmogorovスペクトルに近い
- 中間（$N \sim 1$）: 異方的カスケード
- 強MHD（$N \gg 1$）: Iroshnikov-Kraichnan スペクトル $E(k) \propto k^{-3/2}$（Alfven波の相互作用）
磁場が強いとエネルギーカスケードの機構そのものが変わるんですね。
RANSモデルでこれを表現するために、MHD修正を加えた乱流モデルが提案されている。Knaepen & Moreau (2008) の修正 $k$-$\varepsilon$ モデルでは、Joule散逸項が追加される:
この項が乱流エネルギーの電磁的な散逸を表す。
液体金属電池
近年注目されているMHDの応用が液体金属電池だ。溶融金属の正極・負極と溶融塩の電解質で構成される大型蓄電池で、再生可能エネルギーのグリッドストレージ向けに開発されている。
MHDの課題:
- 充放電電流による自己誘導磁場が液体金属を攪拌し、層間の混合（短絡）を引き起こす
- Tayler不安定性: 軸方向電流が臨界値を超えると渦が発生
- Sloshing不安定性: 界面の波動と電磁力の連成
CFDではVOFまたはLevel-Setで3層（正極/電解質/負極）の界面を追跡しながら、電磁力との連成を解く。OpenFOAMベースの専用ソルバーが複数の研究グループで開発されている。
3つの液体金属層が磁場で乱されるのを抑制する設計にCFDが使われるんですね。
核融合炉のMHD
核融合炉のブランケットではPbLi（鉛リチウム合金）やFLiBe（フッ化リチウムベリリウム）などの液体金属が冷却材兼トリチウム増殖材として使われる。
磁場強度が5-10 Tと非常に強く、Hartmann数が $10^4 \sim 10^5$ に達する。このような極高Ha数でのHartmann層・Side層の解像は、通常のCFDメッシュでは不可能に近い。
対策:
- 壁関数的アプローチ: Hartmann層の解析解を壁面境界条件として組み込む
- 準2D（Sommeria-Moreau）モデル: 磁場方向に変化が小さい仮定で2Dに帰着
- スペクトル法: Hartmann層内の急激な変化を高次多項式で精度良く捕捉
最新の研究動向
MHDの最先端はどのあたりですか？
核融合から天体物理まで、MHDの守備範囲は本当に広いですね。
核融合炉のMHD——液体金属冷却材の磁場下流れと「MHDブレーキ効果」
核融合炉（ITER/DEMO）の第一壁冷却では液体金属（リチウムや共晶合金PbLi）が候補冷媒として研究されている。しかし強磁場（4〜12T）中を流れる液体金属にはローレンツ力が働き、流れを磁場方向に引き延ばし横断方向を抑制する「ハートマン効果（Hartmann Effect）」が生じる。これにより熱伝達率が大幅に低下する「MHDブレーキ効果」が発生し、設計熱流束1MW/m²の除熱が困難になる。CFD（MHD方程式とN-S方程式の連成）によるブランケットチャンネル設計が核融合工学の最前線課題で、ITER機構とKIT（ドイツ）が連携してBenchmark計算を公開している。ハートマン数Ha=B·L·√(σ/ρν)が1000を超える核融合条件では通常のRANSモデルが使えず専用MHDソルバが必要だ。
磁気流体力学（MHD）のトラブル対応
MHD計算の典型的トラブル
MHDの計算で特にハマりやすいポイントを教えてください。
MHD特有のトラブルをパターン別に整理しよう。
1. Hartmann層の未解像で圧力損失が過小
症状: MHDによる圧力損失が理論値より大幅に小さい
原因: Hartmann層（厚さ $\delta_H \sim H/\text{Ha}$）のメッシュ解像度が不足。壁面近傍の速度勾配と電流が正確に計算できていない。
対策:
- Hartmann層内に最低5セル以上配置
- Ha = 100, $H$ = 0.1 m の場合: $\delta_H$ = 1 mm、第一層 = 0.2 mm
- Ha = 1000 の場合: $\delta_H$ = 0.1 mm、第一層 = 0.02 mm
- 壁面近傍のプリズム層成長率を1.1以下に抑える
Hartmann数が大きいほどメッシュ要求が厳しくなるんですね。
2. 電流の保存が破れる
症状: $\nabla \cdot \mathbf{J}$ がゼロにならず、非物理的なチャージアップが発生
原因: 電位のPoisson方程式の解が不十分、または電流密度の計算スキームが保存型でない
対策:
- 電位方程式の収束基準を十分に厳しくする（残差 $< 10^{-8}$）
- 電流密度のフェース値を保存型で計算（勾配のフェース中心評価）
- 壁面の電気的境界条件を確認（絶縁壁: $\partial\phi/\partial n = (\mathbf{u} \times \mathbf{B}) \cdot \mathbf{n}$）
3. 磁場のdiv Bエラー（高Rm問題）
症状: 非物理的な磁気モノポール力が発生し、プラズマの挙動がおかしい
原因: 誘導方程式の離散化で $\nabla \cdot \mathbf{B} = 0$ が数値的に保証されていない
対策:
- Constrained Transport: 磁束をフェースで定義し、正確な離散Stokesの定理で更新
- Divergence Cleaning: $\psi$ を導入し、$\partial\mathbf{B}/\partial t + \nabla\psi = ...$ で磁気モノポールをダンプ
- Powell 8-wave法: ソース項に $-\nabla \cdot \mathbf{B}$ 比例の項を追加して安定化
div B = 0 の維持がMHD計算の命なんですね。
4. 乱流モデルが不適切
症状: MHDで流れが層流化するはずの領域で乱流粘度が過大に計算される
原因: 標準 $k$-$\varepsilon$ モデルはMHDによる乱流抑制効果を考慮していない
対策:
- MHD修正付き乱流モデルを使用（Fluent MHD module に含まれる）
- LESを使用（MHD効果がスケール分解で自然に考慮される）
- DNS（可能な場合）
- Stuart数 $N = \text{Ha}^2/\text{Re}$ が $> 1$ なら、乱流を無視して層流計算から始める
デバッグの手順
MHD計算のシステマティックなデバッグ手順:
1. 磁場なし（Ha = 0） で流れ場を収束させる
2. 弱い磁場（Ha = 10程度） を印加し、Hartmann流れとの比較で実装を検証
3. 磁場を段階的に目標値まで上げる
4. 各段階で圧力損失を理論値と比較: $\Delta p_{\text{MHD}}/\Delta p_0 \propto \text{Ha}^2$ （高Ha極限）
5. 電流密度の分布を可視化し、物理的に妥当か確認
また段階的なアプローチですね。Hartmann数をゼロから徐々に上げていく。
MHDは方程式の連成が強いので、いきなり高Ha数で計算を始めると発散するリスクが高い。Hartmann流れという確立された検証問題があるのは大きなアドバンテージだ。活用しない手はないよ。
MHD CFDが発散する——ローレンツ力項の陰的・陽的処理の選択ミス
MHD（磁気流体力学）CFDで頻出する発散問題：ローレンツ力項J×BをN-S方程式に陽的（Explicit）に加算すると、電流密度Jと速度場uの連成が不安定になりCFLが非常に小さくなければ解が発散する。ハートマン数Ha>100の強磁場条件では特に顕著で、タイムステップをHa²/Re倍に小さくしないと安定しない——これが現実的な計算時間内では不可能な小ささになることがある。解決策は「陰的磁気圧力補正（Implicit Lorentz Force）」スキームを使用するか、誘導方程式（Magnetic Induction Equation）を速度場と完全連成で解く「Full MHD」ソルバに切り替えること。低Rm（磁気Reynolds数<1）の液体金属条件では誘導磁場を無視した「Low-Rm近似」が適用でき、計算を大幅に単純化できる。
関連トピック
なった
詳しく
報告