蓄熱材（PCM）の熱シミュレーション

いい質問だ。PCM（Phase Change Material）は固体から液体に変わるとき——つまり融けるときに大量の熱を吸収する材料だ。この「潜熱」を利用して温度をコントロールするのが蓄熱の本質なんだ。

例えば、氷が0°Cで融けるとき、温度は上がらないまま334 kJ/kgもの熱を吸収するだろう？ PCMはこの原理を利用する。パラフィン系なら融点28°C付近で潜熱200 kJ/kg、塩水和物なら融点32°C付近で潜熱260 kJ/kgくらいだ。電池パックの温度を40°C以下に保つために、セル間にPCMを充填する設計が近年急増している。

固体と液体の境界——融解フロントがどう移動するか、いつどこまで融けるか、全部融けたあとに温度がどう上昇するかを予測する。これが「ステファン問題」と呼ばれる古典的な移動境界問題だ。CAEではエンタルピー法を使って、境界を明示的に追跡せずに解くのが主流なんだよ。

そうだ。1次元のステファン問題を考えてみよう。片面が高温に加熱されたPCMスラブがあるとする。固液界面の位置を $s(t)$ とすると、界面では「潜熱を吸収しながら界面が前進する」という条件が成り立つ：

その通り。だから実務ではほとんどの場合、界面追跡を回避するエンタルピー法を使うんだ。「界面がどこにあるかは知らなくてもいい、エネルギー収支さえ合っていれば温度場は正しく出る」という発想だよ。

エンタルピー法では、体積あたりのエンタルピー $H(T)$ を変数とするエネルギー保存式を解く：

正解！ 純物質なら融点で垂直にジャンプする。パラフィンや塩水和物のように融解温度範囲がある材料では、傾斜した階段状になる。このジャンプの高さが潜熱 $L$ に相当するんだ。シミュレーション上は、この急な変化が非線形性の原因になるから、時間刻みとメッシュには注意が必要だよ。

もう一つ有名な方法が等価比熱法（Apparent Heat Capacity Method）だ。潜熱を融解温度範囲の見かけの比熱 $c_{p,eff}$ としてモデリングする：

実務的にはエンタルピー法が圧倒的に有利だ。等価比熱法は時間刻みが粗いと、ちょうど融解温度域を飛び越えてしまって潜熱を取り逃す——いわゆる「潜熱スキップ」が起きる。Ansys FluentやCOMSOLの相変化モジュールは内部的にエンタルピー法を使っていて、等価比熱法はFEMの標準熱解析ソルバーを使い回すときの簡易手法だね。

ステファン数は顕熱と潜熱の比を表す無次元数だ。PCMシミュレーションの「難しさ」を定量化する指標になる：

ざっくり言うと：

大きく3つに分類される。シミュレーションに必要な物性値とともに整理しよう：

そう、これがPCMシミュレーション最大の実務的課題だ。熱伝導率が低いと融解速度が極端に遅くなる。だから実際の製品では、アルミフィンを入れたり、カーボンナノチューブを分散させたりして実効的な熱伝導率を上げる。シミュレーションでは「PCM＋フィン」の複合体として扱うか、均質化した実効物性値を使うか、ここが設計者の腕の見せどころだね。

FVM（有限体積法）なら、セルごとのエンタルピーを保存量として直接扱えるから相性がいい。Ansys FluentやSTAR-CCM+はこの方式だ。FEM（有限要素法）では、弱形式に変換してから離散化する。

PCMの場合は陰解法（後退Euler）が鉄則だ。理由は2つ：

固体でも液体でもない中間状態——固液が共存する領域だ。パラフィン混合物だと融点が28〜32°Cの範囲に広がるから、その温度帯全体がマッシーゾーンになる。

実は $C_{mush}$ の選定はPCMシミュレーションの「職人芸」の部分だ。純パラフィンなら $10^5$ 〜 $10^6$ が多く、実験での融解フロント位置と合わせて調整する。$C_{mush}$ が大きすぎるとマッシーゾーンが硬くなりすぎて融解が遅くなり、小さすぎると固相領域で非物理的な流れが出る。論文では「$C_{mush}$ 感度解析」を必ず行っているから、3水準以上試すことを勧めるよ。

めちゃくちゃ動く。これがPCMシミュレーションの難しいところだ。融けたパラフィンは温度差で密度差が生じて、浮力駆動の自然対流が発生する。Rayleigh数が $10^5$ を超えると対流効果が支配的になって、純伝導の予測より2〜5倍速く融解が進むことがある。

その通り。縦型のPCM容器だと特に顕著で、上部は対流で速く融けるのに下部は伝導だけでゆっくり融ける——非対称な融解パターンになる。「純伝導モデルで計算したら実験の2倍の時間がかかった」というのは、対流を無視した典型的な失敗パターンだよ。

PCM特有のチェックポイントを挙げよう：

EV用リチウムイオン電池パックのPCM冷却シミュレーションの典型的なフローを教えよう：

正解。融点が高すぎるとPCMが融けないから蓄熱が始まらない。低すぎると周囲温度で勝手に融けてしまって肝心なときに蓄熱容量が残っていない。夏場の車内温度（50〜60°C）を考慮すると、融点32〜35°Cが実務上のスイートスポットだ。

もう一つ重要なポイントがある。PCMが全部融けた後——完全融解後の温度上昇だ。パラフィンの比熱は液相で約2,000 J/(kg·K)と水の半分以下。潜熱を全部使い切ったら、あとは急激に温度が上がる。だからシミュレーションでは「何サイクルまでPCMが保つか」を必ず検証する。現場では「3サイクルの急速充電に耐えられるか」がよくある設計要件だよ。

PCM内蔵石膏ボードが代表例だ。壁に融点23°Cのマイクロカプセル化PCMを練り込んで、昼間の太陽熱を蓄え、夜間に放出して室温変動を抑える。EnergyPlusやTRNSYSで年間エネルギーシミュレーションを行うのが標準的なアプローチだ。

建築スケールではCFDは計算コストが見合わないからね。EnergyPlusのConFDアルゴリズムは1次元のエンタルピー法を壁層の各ノードに適用する簡易モデルだ。PCMの$H(T)$カーブをCSVで入力するだけで、年間8,760時間のシミュレーションが数分で終わる。冷房負荷の15〜25%削減が典型的な結果だよ。

PCM特有のメッシュ要件をまとめよう：

PCMの固液相変化に対応した主要ソルバーを比較しよう。ポイントは「エンタルピー法のネイティブ対応」と「液相の自然対流が解けるか」の2点だ：

Abaqusは構造解析がメインだからね。PCMの熱伝導だけを解くなら *LATENT HEAT キーワードで簡単に設定できるけど、液相の対流は無視される。薄型のPCMシート（厚さ1〜3 mm）なら対流の影響が小さいから問題ないが、大きなPCM容器だと融解時間を30〜50%過大に予測するリスクがある。FluentやSTAR-CCM+なら最初から流れと熱を連成で解くから安心だ。

OpenFOAM v2306以降では、buoyantPimpleFoamに相変化のソースタームを追加するアプローチが使える。具体的には fvOptions で solidificationMeltingSource を指定する：

急速に広がっている。特に注目されているのは以下の3方向だ：

Al₂O₃やCuOのナノ粒子をパラフィンに1〜5 wt%分散させると、実効熱伝導率が20〜50%向上する。シミュレーションでは Maxwell モデルで実効物性値を計算する：

その通り。論文では $\phi = 2\sim 3\%$ がスイートスポットという結論が多い。それ以上入れると粘度が急上昇して対流が死んでしまい、せっかくの熱伝導率向上が打ち消される。これを正しく予測するには、対流を含めた連成解析（Fluent/STAR-CCM+レベル）が必須だ。

あるあるだ。PCMシミュレーションの収束問題は9割が以下の3つに集約される：

いいところに気づいた。物理的にはそうだけど、数値的には $T_{solidus} = T_{liquidus}$ にするとエンタルピーのジャンプが無限に急峻になって、どんなに時間刻みを小さくしても安定しない。$\Delta T = 0.5$°C程度の「数値的スムージング」は実務上必須で、最終的な融解時間への影響は1%以下だよ。

エネルギーバランスをチェックするのが確実だ。全領域のエンタルピー変化量を時間積分して、入熱量と比較する。もし $\Delta H_{total}$ が $\int Q \, dt$ よりも $m \cdot L$ だけ少なかったら、潜熱が丸ごとスキップされている証拠だ。対策は：

それは「過冷却（supercooling）」で、バグではなく実際に起きる物理現象だ。特に塩水和物（CaCl₂·6H₂O など）は5〜15°Cの過冷却を示すことがある。つまり融点32°Cなのに20°Cまで冷えてからやっと凝固が始まる。

塩水和物を扱うなら入れるべきだ。パラフィンはほぼ過冷却しないから無視してOK。COMSOLには「Supercooling」オプションが用意されていて、核生成温度 $T_{nucleation}$ を指定できる。Fluentで再現するにはUDFで凝固開始条件を制御する必要があるから、ちょっと手間がかかる。