磁気流体力学(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)が使われる。実験データは限られているので、解析解での検証が特に重要だ。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース(下向きの空力的な力)を発生します。つまり理論上、天井に貼り付けて走れる! チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
実務者のための直感的理解
この解析分野のイメージ
CFDって、要は「デジタル風洞」です。自動車メーカーが巨大な風洞実験設備に何億円もかけるところを、PCの中で再現できる。でも1つ注意——風洞実験なら「風を当てれば結果が出る」けど、CFDでは「メッシュの品質」と「乱流モデルの選択」という見えない品質要因がある。ここを手抜きすると、きれいなコンター図が出ても中身はデタラメ…なんてことになりかねません。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め(計算領域)、水の入り口と出口を設計し(境界条件)、ポンプの強さを設定する(流量条件)。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか? 意外と直感的ではありませんか?
初心者が陥りやすい落とし穴
「y+って何ですか?」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない!」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た!」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか(低速)、全開にするか(高速)。でもCFDではもう一つ——「どのくらい暴れた水を出すか」(乱流強度)も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね? CFDでも入口条件のミスは下流全体に波及します。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「磁気流体力学(MHD)をもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →