エンタルピー法による相変化解析
理論と物理
概要
先生、エンタルピー法って聞いたんですけど、普通の熱伝導解析と何が違うんですか？
固体が溶けたり固まったりするときは、温度が変わらないのに熱を吸収・放出するだろう？ その「潜熱」を温度場の方程式に自然に組み込む手法がエンタルピー法だよ。鋳造の凝固解析やPCM（相変化材料）の蓄熱設計で必須の技術だ。
温度が変わらないのに熱が動くって…なんか不思議ですね。普通の熱伝導方程式ではダメなんですか？
ダメなんだ。純物質の融点では温度が一定のまま潜熱が出入りするから、通常の $\rho c_p \partial T/\partial t$ だけでは熱収支が合わない。エンタルピー $H$ を主変数にするか、有効比熱 $c_{eff}$ に潜熱を吸収させるかの2通りのアプローチがある。
具体的にはどんな問題に使われることが多いですか？
代表的な応用は3つある。
- 鋳造凝固解析：溶湯がどこから固まるか、引け巣（収縮欠陥）がどこに発生するかを予測。自動車のエンジンブロック・シリンダーヘッドの品質管理に直結。
- PCM蓄熱システム：室温付近で相変化するパラフィン・塩水和物などを使った熱エネルギー貯蔵（LTES）の設計。データセンター冷却・建材一体型蓄熱に活用。
- 溶接・AM（積層造形）の溶融池：レーザーやアークで生じる局所溶融池の動的挙動予測。スパッタ・残留応力・変形の制御に重要。
支配方程式
エンタルピー法の基本方程式ってどんな形になるんですか？
エンタルピー $H$ を使うと、支配方程式はこうなる。
ここで $H$ は顕熱と潜熱の両方を含む総エンタルピー。液相分率 $f_l$ を導入すると：
$L_f$ は融解潜熱（J/kg）、$f_l$ は液相分率（0〜1）。
$f_l$ って0から1の間を動くんですか？
そうだよ。純物質なら融点で0→1に不連続に変わるが、合金や樹脂では $T_s$（固相線温度）から $T_l$（液相線温度）の間で線形に変化するMushy Zone（糊状域）が存在する：
有効比熱法では潜熱をこの $f_l$ の温度微分で比熱に乗せる：
離散化手法
この方程式をFEMやFVMで解くときの注意点はありますか？
最大の難点は $c_{eff}$ が温度に対して非常に急峻な変化をする点だ。純金属では融点でデルタ関数的なスパイクになる。時間刻み $\Delta t$ と要素サイズを小さくしないと、融点をまたぐ要素で温度が「飛び越える」アンダーシュートが発生する。
- 時間積分：Crank-Nicolson法（θ=0.5）または完全陰的（θ=1.0）が安定
- 反復収束：各タイムステップで $T \leftrightarrow H$ の相互更新が必要（ソースタームのライン化）
- メッシュ：固液界面付近を細かく、Mushy Zone幅の1/5以下の要素サイズ
「温度が融点を飛び越えるアンダーシュート」というのを防ぐ具体的な基準はありますか？
Stefan数（Ste）が参考になる：
$Ste \ll 1$ のとき（潜熱が支配的）は時間刻みを小さく取る必要がある。目安として1タイムステップでのセル内の温度変化が Mushy Zone幅（$T_l - T_s$）の1/10以下になるように $\Delta t$ を設定する。アルミ合金の鋳造では $T_l - T_s \approx 50$℃、$\Delta t < 0.01$s（凝固フロント速度に依存）が典型値だ。
数値解法と実装
数値手法の詳細
エンタルピー法と有効比熱法って、どっちを使えばいいんですか？
ざっくり言うと、実務では有効比熱法（または潜熱源項法）が主流だ。理由はシンプルで、既存の熱伝導ソルバーに手を入れやすいから。商用コードのほとんどはこちらを採用している。
| 手法 | 主変数 | メリット | デメリット |
|---|---|---|---|
| エンタルピー法 | $H$ | 純物質の不連続融点に対応 | $T(H)$ の逆変換が必要 |
| 有効比熱法 | $T$ | 既存ソルバーに統合しやすい | 融点スパイクで数値誤差が大きい |
| 潜熱源項法 | $T$ | バランスが良い | 収束に反復が必要 |
収束性と安定性
計算が発散したりオシレーションしたりするのを防ぐコツはありますか？
実務でよく使われる安定化テクニックを紹介するね。
- 温度スムージング（スミアリング）：Mushy Zone幅を実測より±3〜5℃広げる。比熱スパイクの高さを下げる効果がある。
- アンダーリラクゼーション：$f_l$ の更新に緩和係数（0.7〜0.9）を適用する。
- タイムステップ制御：相変化が起きているセルでは $\Delta t$ を自動縮小するアダプティブ法。
- 予熱解析：まず定常解析で温度場を初期化してから過渡解析に移行する。
Mushy Zone内の流動モデル
溶湯が凝固していく過程で、液体の流れ（対流）はどう扱うんですか？ 固まりかけのMushy Zoneも流れる？
Fluent等の流体ソルバーでは、Mushy Zoneを多孔質体としてモデル化するCarman-Kozeny法が標準的だよ：
$A_{mush}$ は糊状域定数（典型値 $10^4$〜$10^7$）、$\varepsilon$ は小数（ゼロ除算防止）。液相では $f_l = 1$ でソース項がゼロ（自由流れ）、固相では $f_l = 0$ でソース項が無限大に近くなり流速がゼロに強制される。この手法でRANS流れと相変化を同じソルバーで一括解析できる。
$A_{mush}$ の値はどう決めればいいですか？ 感度が高そうですね。
$A_{mush}$ は大きい（$10^7$）ほど固体領域をより剛体に近く扱える。ただし大きすぎると収束が悪化する。実際には凝固速度や引け巣位置の実験結果と照合して感度解析でキャリブレーションするのが実務的なアプローチだ。純アルミ系では $A_{mush} = 10^5$〜$10^6$、耐熱合金系では $10^6$〜$10^7$ が参考値として使われることが多い。
実践ガイド
鋳造凝固解析のフロー
鋳造の凝固解析を初めてやるんですけど、どこから手をつければいいですか？
典型的な解析フローを紹介しよう。
- 材料データ収集：固相線・液相線温度 $T_s, T_l$、融解潜熱 $L_f$、温度依存の熱伝導率 $k(T)$ と比熱 $c_p(T)$。
- 界面熱抵抗（HTC）の設定：鋳型と溶湯間のエアギャップ形成を考慮した時間依存HTCを用意する。これを見落とすと凝固速度が大幅に狂う。
- メッシュ生成：Mushy Zoneが移動する領域を細かく（1〜2mm以下）。鋳型側も含めた連成メッシュが必要。
- タイムステップ設定：初期は $\Delta t = 0.01$〜0.1sで開始し、収束挙動を見ながら調整。
- 結果検証：凝固フロント速度が実験値と一致するか確認。熱収縮・引け巣の位置を観察。
よくある失敗と対策
先輩が「温度が0℃付近で振動する」って困ってたんですが、これはエンタルピー法のトラブルですか？
典型的な症状だね。純水（融点0℃）の凝固解析で有効比熱法を使うと、融点でのスパイクが原因で温度がオシレーションする。対策は3つ：(1) Mushy Zone幅を0.5〜1℃に設定してスパイクを鈍らせる、(2) 潜熱源項法に切り替える、(3) 時間刻みを1/10以下に小さくする。実務では(1)+(2)の組み合わせが多い。
界面熱抵抗（HTC）を「時間依存」と言っていましたが、具体的にどう設定すればいいですか？
鋳造では金型（ダイカスト・重力鋳造）と溶湯が接触してから凝固・収縮するにつれてエアギャップが形成され、HTCが低下する。典型的なモデルとしては：
- 接触期（凝固開始まで）：HTC = 3000〜5000 W/(m²K)（金属接触に近い）
- エアギャップ形成期：HTCを指数関数的に低下：$h(t) = h_0 e^{-\alpha t}$
- 固化後：HTC = 200〜500 W/(m²K)（空気ギャップ相当）
実験でしか正確な値は得られないが、引け巣の位置・形状が予測と合わない場合はHTCのキャリブレーションが最初の調整ポイントになる。
実践チェックリスト
- 材料データ：Ts, Tl, Lf、温度依存k(T)、cp(T)を揃える
- 手法選択：有効比熱法または潜熱源項法を選択し、Mushy Zone幅を適切に設定
- メッシュ：Mushy Zone移動領域の要素サイズ ≤ Mushy Zone幅の1/5
- 界面熱抵抗（HTC）：鋳型と溶湯の界面に時間依存HTCを設定
- タイムステップ：Steファン数基準でΔtを設定（1ステップでの温度変化 ≤ ΔTmushy/10）
- 流動連成（必要時）：Carman-KozenyモデルのAmush値を検討
- 収束確認：各タイムステップの残差を監視（温度/エンタルピーの収束基準）
- 結果検証：凝固フロント位置・凝固時間を実測値と照合
ソフトウェア比較
主要なCAEツールでのエンタルピー法の実装はどう違うんですか？
代表的なツールの比較を見てみよう。
| ツール | 手法 | Mushy Zone | 備考 |
|---|---|---|---|
| Ansys Fluent | 潜熱源項法 | Carman-Kozeny多孔質モデル | SOLIDIFICATION/MELTINGモジュール |
| Ansys Mechanical | 有効比熱法 | 線形補間 | 材料DBで自動スミアリング（v2020以降） |
| COMSOL | 両方選択可 | 任意関数定義 | Phase Change Interfaceが直感的 |
| OpenFOAM | 潜熱源項法 | buoyancyFoam系ソルバー | solidificationMeltingSourceクラス |
| ProCAST | FEM + エンタルピー法 | 実測相図連携 | 鋳造専用、熱収縮・引け巣予測に強い |
鋳造専用ツールって強みがあるんですね。汎用CAEと使い分けるんですか？
そう。ProCAST・MAGMASOFTのような鋳造専用ツールは相図データベースや引け巣予測アルゴリズムが充実している。一方、Fluent/COMSOLは流体との連成（溶湯の対流）を扱いやすい。複雑な流動凝固連成には汎用CFDを使い、単純な凝固プロセスなら専用ツールの方がセットアップが速い。
先端技術
エンタルピー法の最新研究トレンドってどんなものがありますか？
いくつか面白いトピックがある。
- 位相場法（Phase-Field法）との融合：エンタルピー法の効率と位相場法の界面追跡精度を組み合わせる。デンドライト成長のマルチスケール解析に応用。
- 機械学習ベースのMushy Zone推定：多元合金の相図をPINNで学習し、より正確な $f_l(T)$ を予測。
- 格子ボルツマン法（LBM）との統合：格子レベルの相変化を自然に扱えるLBMへのエンタルピー法の組み込みが活発。
- AM（積層造形）の溶融池解析：レーザーPBF・DED工程での急速溶融・凝固のマルチスケール解析。
積層造形（AM）の溶融池解析では何が難しいんですか？
スケールのギャップが最大の問題だ。レーザービーム径は10〜100μm、溶融池サイズは100〜500μm程度だが、この溶融池が動き回って部品全体（100mm〜1mスケール）の残留応力・変形に影響する。全域を溶融池サイズで解くと計算量が天文学的になる。そのためマルチスケールアプローチが必須で：
- 溶融池スケール：エンタルピー法＋Marangoni対流でメルトプール形状・スパッタを予測
- 積層スケール：均質化モデル（固有ひずみ法）で残留応力・変形を計算
- AIによるブリッジ：溶融池特性からマクロひずみへの写像をNNで学習
この分野は2020年代の製造業DXの中で最も研究投資が盛んな領域のひとつだ。
エンタルピー法の誕生秘話
エンタルピー法はKillworth（1973年）とVoller & Cross（1981年）によって発展した。当時、Stefanの移動境界問題を陽的に解くと界面位置の追跡が非常に困難だった。「界面を追わずに、エンタルピーという量を解けばよい」というアイデアは単純だが革命的で、今日の相変化CAEの礎を築いた。
トラブルシューティング
エンタルピー法特有のトラブルはどんなものがありますか？
詳しいトラブルシューティングガイドを用意しているよ。主要な問題パターンを確認してほしい。
| 症状 | 主な原因 | 対策 |
|---|---|---|
| 融点付近での温度オシレーション | 有効比熱法の比熱スパイク | Mushy Zone幅を広げる、潜熱源項法に変更、Δtを縮小 |
| 凝固フロントが遅すぎ/速すぎ | HTC設定ミス（鋳型-溶湯界面） | HTCを実測値ベースで再設定、感度解析 |
| 引け巣位置が実際と一致しない | Amush（多孔質定数）の不適切な値 | Amushの感度解析（10^4〜10^7の範囲で確認） |
| 収束が遅い・発散する | flの更新発散 | アンダーリラクゼーション係数0.5〜0.8に設定 |
| 計算時間が想定より大幅に超過 | Mushy Zone内の要素数過多 | Mushy Zone外は粗いメッシュに、適応細化（AMR）を使用 |
温度オシレーション、収束失敗、界面スメアリングなどのトラブル別解決策