磁気流体力学（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安定性解析の難しさを物語っています。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
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%を超えると数値的な「磁気モノポール」が生まれ、プラズマが誤った方向に飛んでいく。磁場の発散ゼロは物理の要請であり、数値解法の戦場です。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
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%向上した事例が鉄鋼学会誌で報告されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
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分野の特徴です。
選定で最も重要な3つの問い
- 「何を解くか」：磁気流体力学（MHD）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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特有のトラブルをパターン別に整理しよう。
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近似」が適用でき、計算を大幅に単純化できる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——磁気流体力学（MHD）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告