建築用PCM（相変化蓄熱材）の熱シミュレーション
理論と物理
PCMとは何か
先生、建築にPCMを使うって聞いたんですけど、壁が溶けるんですか？ なんか怖いんですけど...
壁自体が溶けるわけじゃないよ。PCM（Phase Change Material＝相変化蓄熱材）っていうのは、固体から液体に変わるときに大量の熱を吸収する材料のこと。建築では、パラフィンワックスの一種をマイクロカプセル化して直径5〜50μmの粒にして、石膏ボードやモルタルに練り込むんだ。
マイクロカプセル？ つまりパラフィンが殻に閉じ込められてて、その中でだけ溶けたり固まったりするってことですか？
その通り。例えばBASFのMicronal DSシリーズは融点23°C付近のパラフィンをメラミン樹脂の殻で包んでいる。日中、室温が23°Cを超えると壁の中のPCMが溶け始めて太陽熱を潜熱として蓄える。夜になって室温が下がると今度は固まって放熱する。この「蓄熱→放熱」サイクルで、室温変動を2〜4°C抑制できるんだ。
へぇ、壁が「熱の緩衝材」になるわけですね。でも普通のコンクリートでも蓄熱しますよね？ PCMのメリットって何ですか？
いい質問だね。コンクリートの蓄熱は顕熱——つまり温度が上がることで蓄える。これだと蓄熱密度が低い。水の顕熱は1°C上昇あたり4.18 kJ/kgだけど、パラフィンC18の溶解潜熱は約244 kJ/kg。同じ質量で比較すると、PCMは温度変化たった1〜2°Cの範囲で、コンクリートが50°C上昇して蓄える量に匹敵する熱を蓄えられる。しかも室温の「快適帯」付近で作動するから、冷暖房負荷の削減に直結するんだ。
PCMの蓄放熱における全エンタルピー変化は、固相の顕熱＋潜熱＋液相の顕熱で表される。
ここで $m$ は質量、$c_s, c_l$ は固相・液相の比熱、$T_m$ は融点、$L_f$ は溶解潜熱、$T_1, T_2$ は初期・最終温度である。建築用パラフィン系PCMの代表的な物性値を以下に示す。
| 物性 | 値 | 単位 |
|---|---|---|
| 融点 $T_m$ | 21〜26 | °C |
| 溶解潜熱 $L_f$ | 170〜250 | kJ/kg |
| 密度 $\rho$（固体） | 850〜950 | kg/m³ |
| 熱伝導率 $k$ | 0.15〜0.25 | W/(m·K) |
| 比熱 $c_p$（固/液） | 1.8〜2.4 | kJ/(kg·K) |
Stefan問題と移動境界
PCMが溶けたり固まったりするとき、固相と液相の境界が動くんですよね？ それってどう定式化するんですか？
それがStefan問題だ。1891年にオーストリアの物理学者Jozef Stefanが定式化した、相変化を伴う移動境界問題。固相領域と液相領域それぞれで熱伝導方程式を解き、境界面 $x = s(t)$ でのエネルギー保存条件（Stefan条件）を満たす。1次元だとこうなる。
固相領域 ($0 < x < s(t)$):
液相領域 ($s(t) < x < L$):
Stefan条件（界面 $x = s(t)$ でのエネルギーバランス）:
界面の位置 $s(t)$ が未知数で、しかも時間とともに動くんですね。これって解析的に解けるんですか？
半無限固体の片側加熱みたいな特殊条件なら、Neumann解 $s(t) = 2\lambda\sqrt{\alpha t}$ という解析解がある。でも建築壁体のような有限厚さ・両面温度変動・複数材料層の問題では解析解は無理。数値解法が必須で、そこで登場するのがエンタルピー法だ。
エンタルピー法の支配方程式
エンタルピー法って、Stefan問題を解く便利な方法なんですか？ どういう発想なんでしょう？
Stefan問題の厄介なところは「界面位置の追跡」だ。エンタルピー法はこの問題を回避する天才的なアイデアで、温度 $T$ の代わりに体積エンタルピー $H$ を従属変数にする。相変化の前後で $H(T)$ 関数が潜熱分だけジャンプするから、固相・液相・混合相（mushy zone）を単一の方程式で統一的に扱える。界面を明示的に追跡する必要がないんだ。
エンタルピー法の支配方程式は以下のように書ける。
ここで体積エンタルピー $H(T)$ は温度の関数として定義される。
建築用PCMでは、相変化が単一温度ではなく $T_s$ から $T_l$ の温度範囲（通常2〜4°C幅）にわたって起こる。この範囲がmushy zoneであり、固相と液相が共存する領域だ。
なるほど！ 温度で方程式を書くと界面で不連続になって困るけど、エンタルピーで書けば滑らかに扱えるってことですね。
その通り。mushy zone幅 $\Delta T = T_l - T_s$ がゼロに近づくと理想的な「鋭い界面」のStefan問題に帰着する。実際の建築用PCMは $\Delta T = 2\sim4$°Cの有限幅を持つから、むしろエンタルピー法が自然な定式化なんだ。
PCM壁体の伝熱メカニズム
実際の建築壁体って、PCM層だけじゃなくて断熱材とかコンクリートとか色々な層がありますよね。全体としてどうモデル化するんですか？
典型的なPCM壁体の構成を例で説明しよう。外側から「外装材（10mm）→ 断熱材（50mm）→ PCM含有石膏ボード（12.5mm）→ 空気層→ 内装仕上げ」みたいな多層構造だ。各層の1次元非定常熱伝導を連立して解く。
多層壁体の各層 $i$ について:
層間の接続条件（温度の連続性と熱流束の保存）:
外面・内面の境界条件には、対流（$h_\mathrm{ext}$, $h_\mathrm{int}$）と日射吸収（$\alpha_s I_\mathrm{sol}$）、長波放射（$\varepsilon \sigma (T^4 - T_\mathrm{sky}^4)$）を含める。
PCM層だけ $c_p(T)$ が温度依存になって非線形になるわけですね。他の層は通常の熱伝導方程式のままで済むと。
そうだね。だから計算上のポイントは「PCM層をどれだけ丁寧に離散化するか」に集約される。断熱材やコンクリートは粗い分割でも十分だけど、PCM層は相変化の温度幅 $\Delta T$ の中で急激にエンタルピーが変化するから、細かい分割が必要になる。
Stefan問題の歴史
Stefan問題の名は、1891年にボルツマンの師匠でもあったオーストリアの物理学者Jozef Stefan（1835-1893）に由来する。彼が取り組んだのは「北極海の氷がどれだけの厚さに成長するか」という問題だった。氷の表面温度と海水温度が与えられたとき、氷の厚さ $s(t)$ が $\sqrt{t}$ に比例して成長するというNeumann解は、今日でもPCMシミュレーションの検証用ベンチマーク問題として使われている。建築分野への応用が本格化したのは1980年代で、米国エネルギー省（DOE）のDoris Boehmらの研究がきっかけだ。
各項の物理的意味
- エンタルピー蓄積項 $\partial H / \partial t$：単位体積あたりのエンタルピー変化率。顕熱（温度変化）と潜熱（相変化）の両方を含む。PCMの融点付近では、わずかな温度変化に対して $H$ が急激に増大する——これが「熱を吸い込む」効果。日常の例で言えば、氷水が0°Cで長時間留まるのは、氷が溶ける潜熱を吸収し続けているから。
- 熱伝導項 $\nabla \cdot (k \nabla T)$：フーリエの法則に基づく熱拡散。PCMの熱伝導率は0.15〜0.25 W/(m·K)と低く、コンクリート（1.4 W/(m·K)）の1/6程度。このため、PCM層は「熱を蓄えるのは得意だが、伝えるのは遅い」という特性を持つ。壁体設計では、PCMの蓄熱能力を活かしつつ、内面側への放熱経路を確保するバランスが重要になる。
- 熱源項 $Q$：内部発熱（一般にPCM壁体では $Q=0$）。ただし、電気ヒーター内蔵の能動型PCMパネルではジュール発熱 $Q = I^2R/V$ を考慮する場合がある。
仮定条件と適用限界
- 1次元伝熱の仮定：壁体の面積に比べて厚さが十分小さい場合、面方向の熱伝導は無視でき、厚さ方向のみの1D問題として扱える。窓枠周辺やコーナー部はヒートブリッジとなるため2D/3D解析が必要。
- 自然対流の無視：マイクロカプセル型PCMでは液相が微小カプセル内に閉じ込められているため、マクロな自然対流は生じない。ただし、バルクPCM（タンク蓄熱等）では液相の対流が支配的になり、Navier-Stokes方程式との連成が必要。
- 体積変化の無視：パラフィンは相変化時に約10%の体積変化を伴うが、マイクロカプセルの弾性変形で吸収されるため、構造的影響は通常無視できる。
- ヒステリシスの無視：実際のPCMは融解と凝固で異なる温度-エンタルピー曲線（ヒステリシス）を示すが、多くのBESツールでは無視される。精密な評価にはヒステリシスモデルが必要。
次元解析と無次元数
| 無次元数 | 定義 | 物理的意味 | PCM壁体での典型値 |
|---|---|---|---|
| Stefan数 $\mathrm{Ste}$ | $c_p \Delta T / L_f$ | 顕熱と潜熱の比 | 0.02〜0.1（潜熱支配） |
| Biot数 $\mathrm{Bi}$ | $h L / k$ | 表面対流と内部伝導の比 | 0.5〜5 |
| Fourier数 $\mathrm{Fo}$ | $\alpha t / L^2$ | 無次元時間（拡散の進行度） | 年間サイクルでは $\gg 1$ |
数値解法と実装
エンタルピー法の離散化
エンタルピー法を実際にコンピュータで解くには、どう離散化するんですか？
建築のBES（Building Energy Simulation）ツールでは、壁体を厚さ方向に有限差分法（FDM）で離散化するのが主流だ。EnergyPlusのCondFDアルゴリズムを例にすると、PCM層を $N$ 個のノードに分割して、各ノード $j$ で以下の陰解法差分式を解く。
$H$ と $T$ の両方が未知数だから、$H(T)$ の関係式で繋ぐわけですね。でもこれ非線形ですよね？
鋭いね。各時間ステップで $H^{n+1}$ から $T^{n+1}$ を逆関数 $T = H^{-1}(H)$ で求める必要がある。具体的には、まず前の時間ステップの温度を初期推定にして、$H(T)$ 曲線上で温度を更新する反復計算を行う。mushy zone 内では $T = T_s + (H - H_s)\Delta T / (\rho L_f)$ で明示的に解けるから、反復は1〜3回で収束することが多い。
見かけ比熱法との比較
「見かけ比熱法」というのも聞いたことがあるんですが、エンタルピー法とどう違うんですか？
見かけ比熱法（Apparent Heat Capacity法）は、潜熱を温度依存の等価比熱 $c_\mathrm{eff}(T)$ として表現する方法だ。
見かけ比熱法を使えば、通常の熱伝導方程式 $\rho c_\mathrm{eff}(T) \partial T / \partial t = \nabla \cdot (k \nabla T)$ として解けるから、既存のソルバーに実装しやすい。ただし大きな落とし穴がある——時間刻みが粗いと、相変化の温度範囲を「飛び越えて」しまい、潜熱が計算から丸ごと抜け落ちる可能性があるんだ。
え、それは怖い！ つまり1ステップで $T_s$ から $T_l$ を飛び越えると、潜熱分のエネルギーが消えるってことですか？
その通り。これが見かけ比熱法の最大の弱点。一方、エンタルピー法は $H$ を直接追跡するから、時間刻みが多少粗くても潜熱の保存が保証される。実務では以下のように使い分ける。
| 手法 | メリット | デメリット | 推奨場面 |
|---|---|---|---|
| エンタルピー法 | エネルギー保存が厳密、ロバスト | $H \to T$ の逆変換が必要 | BESツール（EnergyPlus等） |
| 見かけ比熱法 | 既存ソルバーに実装容易 | 時間刻み依存、潜熱の見落としリスク | CFD/FEM（COMSOL、Fluent等） |
時間刻みとメッシュ分割
PCM壁体の解析で適切な時間刻みとメッシュの細かさってどのくらいですか？
これは結構シビアなポイントだ。PCMの熱拡散率は $\alpha = k/(\rho c_p) \approx 8 \times 10^{-8}$ m²/s と非常に小さい。12.5mm厚のPCM石膏ボードに対する特性時間は $\tau = L^2/\alpha \approx 2400$ 秒（40分）。これを基準にすると:
- 時間刻み: 1〜3分（年間シミュレーションでは3分が実用的、精密評価では1分）
- 空間分割: PCM層は最低4分割（$\Delta x \le 3$ mm）、推奨は6〜8分割
- Fourier数の制約: $\mathrm{Fo} = \alpha \Delta t / \Delta x^2 \le 0.5$（陽解法の場合。陰解法では制約なし）
年間シミュレーションだと 365日 $\times$ 1440分 = 525,600ステップ... 3分刻みでも175,200ステップですね。計算量はどうなんですか？
1D壁体の有限差分だから、1ステップの計算コストは微小。EnergyPlusで1ゾーンモデルの年間シミュレーションなら、PCM壁体を含めても数十秒〜数分で終わる。多ゾーン建物全体でも10分程度だ。COMSOLで2D/3D詳細解析する場合はケタ違いに重くなるけど、壁体内部の自然対流や界面近傍の詳細を見たい場合以外は1D-BESで十分。
非線形収束の扱い
PCM層の非線形性って、収束に問題を起こしませんか？
確かに、$c_\mathrm{eff}(T)$ が融点付近で急激に変化するから、Newton-Raphson法で解くとヤコビアンの条件数が悪化する。実務的な対策は3つある。
- mushy zone幅を広げる: $\Delta T$ を実測値（通常2〜4°C）に忠実に設定する。人為的に0.1°Cなど狭くすると発散しやすい
- エンタルピー法を採用する: 見かけ比熱法よりもロバスト
- 緩和係数の導入: $T^{n+1}_\mathrm{new} = \omega T^{n+1}_\mathrm{calc} + (1-\omega) T^n$、$\omega = 0.5\sim0.8$
PCM壁体のたとえ
PCM壁体は「壁に埋め込まれた氷嚢（ひょうのう）」だと思えばよい。夏の暑い日、氷嚢を額に当てると氷が溶けきるまで涼しさが持続する——これが潜熱の効果。PCM壁体も同じで、日中の太陽熱で壁内のPCMが「溶ける」ことで室温上昇を抑え、夜に外気温が下がると「固まって」蓄えた熱を外に放出する。ただし氷嚢と違って室温付近（23°C前後）で作動するのがポイント。「快適温度で働く壁の中の小さな氷嚢」が建築用PCMの本質だ。
実践ガイド
EnergyPlusでのPCM設定
EnergyPlusでPCM壁体をモデル化するには、具体的にどうすればいいですか？
EnergyPlus（v24.1以降）では、MaterialProperty:PhaseChange オブジェクトで温度-エンタルピー曲線を定義する。手順はこうだ。
- 壁体構成の定義:
Constructionオブジェクトで多層壁を定義。PCM層は通常のMaterialとして基本物性を入力 - PCM物性の追加:
MaterialProperty:PhaseChangeで温度-エンタルピーの対（ペア）を8〜16点指定 - CondFDの有効化:
HeatBalanceAlgorithmをConductionFiniteDifferenceに設定（デフォルトのCTFではPCMを扱えない） - 時間刻みの設定:
Timestepを20（= 3分間隔）以上に設定
温度-エンタルピー曲線の入力って、具体的にはどんな数値を入れるんですか？
BASF Micronal DS 5040 X（融点23°C、潜熱110 kJ/kg）を例にすると、こんな感じだ。
| 温度 [°C] | エンタルピー [J/kg] | 説明 |
|---|---|---|
| -20 | 0 | 基準点 |
| 15 | 31,500 | 固相の顕熱蓄積 |
| 20 | 40,500 | 相変化開始直前 |
| 21 | 55,000 | mushy zone 前半 |
| 23 | 120,000 | 融点（潜熱ピーク） |
| 25 | 155,000 | mushy zone 後半 |
| 26 | 162,000 | 相変化完了 |
| 50 | 210,000 | 液相の顕熱蓄積 |
重要なのは、この曲線が必ず単調増加でなければならないこと。エンタルピーが減少する点があるとEnergyPlusが警告なしに異常な結果を出す場合がある。あと、CondFDの空間ノード数は CondFD:SpaceDiscretizationConstant で制御するけど、デフォルト値3だとPCM層の分割が粗すぎる。PCM含有層は0.5以下に設定して細かく分割しよう。
TRNSYSでのPCMモデル
TRNSYSでPCMを扱う場合はどうですか？ EnergyPlusとの違いはありますか？
TRNSYSではType 399（Kuznikら開発）やType 1270（TREFFTZモジュール）でPCM壁体をモデル化する。特徴的なのはヒステリシス対応だ。TRNSYSのPCMモジュールは融解と凝固で異なるエンタルピー曲線を設定できる。
ヒステリシスって、溶ける温度と固まる温度が違うってことですか？
その通り。実際のPCMでは凝固が融解より1〜3°C低い温度で始まる（過冷却）。これを無視すると年間の蓄放熱量を最大15%過大評価してしまうという報告がある。特に中間期（春・秋）で完全な融解-凝固サイクルが成立しない日に影響が大きい。TRNSYSならこの効果を再現できるから、精密な評価が必要な場合はTRNSYSが有利だ。
COMSOL/Ansysでの詳細解析
1D-BESでは解けない問題ってどういうケースですか？
2D/3D解析が必要になるのは主に3つのケースだ。
- ヒートブリッジの評価: 窓枠・バルコニー接合部・鉄骨柱周辺でのPCMの局所的挙動
- バルクPCMの自然対流: 床下蓄熱槽やPCMパネル内部の液相対流が支配的な場合
- マイクロカプセルレベルの伝熱解析: カプセル-マトリクス間の局所温度分布、充填率の最適化
COMSOLなら「Heat Transfer with Phase Change」インターフェースで見かけ比熱法ベースのPCMモデルが直接使える。Gaussianタイプの $c_\mathrm{eff}(T)$ 分布をmushy zone幅とともに自動生成してくれるから設定は楽だ。Ansys Fluentでは Solidification/Melting モデルを使い、VOF法で界面を追跡する。ただし建築用途の年間シミュレーションには重すぎるから、あくまで詳細現象の理解と簡易モデルの検証に使う位置づけだね。
年間エネルギー評価のワークフロー
実際にPCM壁体の省エネ効果を評価するとき、どういう手順で進めるんですか？
典型的なワークフローはこうだ。
- ベースモデルの構築: PCMなしの建物モデルを作成し、年間冷暖房負荷を算出
- PCMパラメトリックスタディ: 融点（18〜28°C）、含有率（10〜30wt%）、設置位置（内壁/外壁/天井）をパラメータとして系統的に変化
- 気象データの選定: ASHRAE TMY3やEA気象データ。少なくとも3都市（暑熱・温暖・寒冷）で比較
- 感度分析: PCM物性（特に潜熱と融点）の±10%変動に対する冷暖房負荷の応答を確認
- 経済性評価: ΔkWh × 電力単価 vs. PCM追加コスト → 回収年数を算出
融点の最適値って気候によって変わりますか？
大きく変わる。ざっくり言うと、冷房負荷削減を狙うなら融点は冷房設定温度 - 1〜2°C（日本の住宅なら24〜26°C）、暖房負荷削減なら暖房設定温度 + 1〜2°C（20〜22°C）が目安。両方を同時に狙うなら融点23°C前後が妥協点になる。東京のように冷暖房の両方が必要な気候では、2種類の融点を持つカスケードPCMが研究されているけど、コストが倍になるから実用化はまだ先だね。
| 気候区分 | 都市例 | 推奨融点 | 主な効果 | 冷房負荷削減 |
|---|---|---|---|---|
| 亜熱帯 | 那覇 | 25〜27°C | 冷房ピークカット | 15〜25% |
| 温暖湿潤 | 東京 | 23〜25°C | 冷房+中間期快適性 | 10〜20% |
| 温帯 | 札幌 | 20〜22°C | 暖房負荷削減 | 5〜15% |
| 乾燥高温 | リヤド | 26〜28°C | 冷房大幅削減 | 20〜30% |
サウジアラビアのPCM壁体実証
キングアブドラ科学技術大学（KAUST）では2015年、外壁にMicronal DS 5040 X（融点23°C）を混合したコンクリートパネルを設置し、2年間の実測を実施した。外気温45°C超の砂漠環境で室内最高温度を3〜5K低減し、空調負荷を12%削減した。EnergyPlus PCMモジュールとの比較では予測精度±0.8°Cを達成している。日本では竹中工務店が2018年に都内オフィスビルでPCM天井パネルの実証試験を行い、夏季ピーク電力の8%削減を報告した。
初心者が陥りやすい落とし穴
「CTFアルゴリズムのままPCMを入れてしまう」——EnergyPlusのデフォルト伝熱アルゴリズムはCTF（Conduction Transfer Function）で、これは線形定常問題向け。PCMの非線形性は扱えず、MaterialProperty:PhaseChange を定義しても潜熱が完全に無視される。CondFDアルゴリズムへの切り替えを忘れると、PCMが全く効かないモデルになる。エラーも警告も出ないから気づきにくい。
検証のコツ
PCMモデルの妥当性確認には、Stefan問題のNeumann解析解と比較するのが定石だ。半無限固体の片面加熱で界面位置 $s(t) = 2\lambda\sqrt{\alpha t}$ と比較すれば、離散化の精度を定量的に確認できる。EnergyPlusの場合は、さらにエネルギーバランスチェックとして、年間の「PCM蓄熱量の積算 = PCM放熱量の積算」を確認する。これがずれている場合、時間刻みかメッシュ分割の問題がある。
ソフトウェア比較
BES（建築エネルギーシミュレーション）ツール
建築のPCMシミュレーションに使えるツールって、どんなものがありますか？
大きく分けて2系統ある。BES（Building Energy Simulation）ツールとCFD/FEMツールだ。まずBES系を見よう。
| ツール | PCMモデル | 手法 | ヒステリシス | ライセンス |
|---|---|---|---|---|
| EnergyPlus | MaterialProperty:PhaseChange | エンタルピー法（CondFD） | 非対応 | 無償（DOE開発） |
| TRNSYS | Type 399 / Type 1270 | エンタルピー法（FDM） | 対応 | 商用 |
| ESP-r | Special material module | 見かけ比熱法 | 非対応 | 無償（OSS） |
| WUFI Plus | PCM material property | エンタルピー法 | 対応 | 商用 |
CFD/FEMツール
CFD/FEM系はどうですか？
| ツール | PCMモデル | 手法 | 自然対流 | 用途 |
|---|---|---|---|---|
| COMSOL Multiphysics | Heat Transfer with Phase Change | 見かけ比熱法（Gaussian） | 対応 | 詳細壁体解析、ヒートブリッジ |
| Ansys Fluent | Solidification/Melting | エンタルピー-多孔体法 | 対応 | バルクPCM、蓄熱槽 |
| OpenFOAM | chtMultiRegionFoam + PCM拡張 | 見かけ比熱法 | 対応 | 研究用途 |
| Ansys Mechanical | Phase Change material property | エンタルピー法 | 非対応 | 熱応力連成 |
ツール選定の指針
結局、どのツールを使えばいいんですか？ 判断基準を教えてください。
目的で使い分けるのが正解だ。
- 年間省エネ評価 → EnergyPlus（無償、実績豊富、1D壁体で十分）
- ヒステリシスを含む精密評価 → TRNSYS（商用だがPCMモデルが高機能）
- ヒートブリッジの2D/3D解析 → COMSOL（マルチフィジクス対応）
- 蓄熱槽内の流動解析 → Ansys Fluent or OpenFOAM
実務ではEnergyPlusで年間評価 → 問題箇所をCOMSOLで詳細解析という2段階アプローチが多い。
建築用PCMの主要メーカー
BASF（ドイツ）のMicronal DSシリーズはマイクロカプセル型で世界最大シェア。Rubitherm（ドイツ）はパラフィン系バルクPCMで研究機関向けに強い。BioPCM（Phase Change Energy Solutions社、米国）は大豆油由来のバイオベースPCMを展開。国内ではDNP大日本印刷がPCMシートを手掛ける。価格はパラフィン系で1kg当たり500〜2,000円、無機塩水和物系で200〜800円が相場で、断熱材の10〜30倍のコスト差が普及の壁となっている。ただし、電力ピークカットによるデマンド料金の削減を考慮すると、オフィスビルでは5〜10年で回収できるケースも出てきている。
ツール選定のたとえ
BESツール（EnergyPlus、TRNSYS）は「年間の家計簿」——毎日の光熱費がいくらかかるかを1年通して集計する道具。CFD/FEMツール（COMSOL、Fluent）は「顕微鏡」——壁体の中で何が起きているかをミクロに観察する道具。家計簿で「8月の電気代が高い」と分かったら、顕微鏡で「なぜ西壁のPCMが効いていないのか」を調べる——この使い分けが効率的だ。
先端技術
カスケードPCMと複合融点設計
PCMの先端技術って、どんな方向に進化しているんですか？
注目されているのがカスケードPCMだ。融点の異なる複数のPCMを壁体の厚さ方向に積層する。例えば外壁側に融点28°CのPCM、内壁側に融点21°CのPCMを配置すると、夏季は外側で日射熱を吸収し、冬季は内側で室内熱を保持する。年間を通した省エネ効果が単一融点PCMより15〜25%向上するという報告がある。
面白い！ でもシミュレーションは複雑になりそうですね。
層ごとに異なる $H(T)$ 曲線を割り当てるだけだから、原理的にはそう難しくない。EnergyPlusでも各PCM層に別々の MaterialProperty:PhaseChange を定義すれば対応できる。ただし最適化は厄介で、融点の組み合わせ、層厚比、設置位置を網羅的に探索する必要がある。最近は遺伝的アルゴリズムや機械学習を使った最適設計が研究されている。
機械学習による最適配置
機械学習をPCMの最適化にどう使うんですか？
典型的なアプローチは、EnergyPlusの年間シミュレーションを数百〜数千ケース実行して訓練データを作り、サロゲートモデル（Gaussian Process Regressionやランダムフォレスト）を構築する。入力変数はPCMの融点、潜熱、含有率、壁体内位置、気候条件。出力は年間冷暖房負荷削減量。サロゲートモデルが完成すれば、ベイズ最適化で最適パラメータを瞬時に探索できる。
最新の研究ではPINN（Physics-Informed Neural Network）を使ってエンタルピー方程式そのものをニューラルネットワークに学習させるアプローチも出てきている。これなら数値解法を使わずにリアルタイムで温度分布を予測できる可能性がある。ただしまだ研究段階で、実務に降りてくるのはもう少し先だ。
ナノ粒子強化PCM
PCMの弱点って熱伝導率が低いことですよね？ それを改善する方法はありますか？
良い着眼点だ。パラフィンの $k = 0.2$ W/(m·K) はアルミ（237 W/(m·K)）の1000分の1以下。そこでナノ粒子強化PCM（NePCM）が研究されている。カーボンナノチューブ、グラフェンナノプレート、Al₂O₃ナノ粒子を1〜5wt%添加すると、熱伝導率が50〜300%向上する。ただしトレードオフがあって、ナノ粒子の添加で潜熱が5〜15%減少し、粘度が増大して充填性が悪化する。
シミュレーション上は、NePCMの有効熱伝導率をMaxwell-Garnett理論やBruggeman式で推定して $k_\mathrm{eff}$ として入力する。
ここで $\phi$ はナノ粒子の体積分率、$k_\mathrm{np}$ はナノ粒子の熱伝導率。
トラブルシューティング
収束失敗と温度振動
先生、PCMシミュレーションでよくあるトラブルって何ですか？ 先輩が「PCMを入れると計算が暴れる」って言ってたんですけど...
最も多いのが融点近傍での温度振動だ。見かけ比熱法で $\Delta T$（mushy zone幅）を狭く設定しすぎると、$c_\mathrm{eff}$ のピークが鋭くなりすぎて数値的不安定を引き起こす。具体的には温度が融点を挟んで毎ステップ振動する。
| 症状 | 原因 | 対策 |
|---|---|---|
| 融点付近で温度が毎ステップ振動 | mushy zone幅が狭すぎる / 時間刻みが粗い | $\Delta T$ を2°C以上に設定、時間刻みを短縮 |
| 潜熱が効いていない（PCMなしと同じ結果） | CTFアルゴリズムのまま（EnergyPlus） | HeatBalanceAlgorithm を CondFD に変更 |
| エネルギーバランスエラーが大きい | 空間分割が粗い | PCM層の分割数を増やす（最低4分割） |
| 年間蓄熱量 $\ne$ 放熱量 | 完全な融解-凝固サイクルが未達成 | 融点が室温帯から外れていないか確認 |
| COMSOLで「ill-conditioned matrix」 | $c_\mathrm{eff}$ の急峻なピーク | Gaussianスムージング幅を広げる |
ヒステリシスの無視による誤差
さっきヒステリシスの話が出ましたけど、無視するとどのくらい誤差が出るんですか？
BrijらのJournal of Building Performance Simulation（2019年）の報告では、ヒステリシスを無視すると年間冷暖房負荷の削減効果を8〜15%過大評価する。特に問題になるのは中間期（春・秋）だ。日較差が小さい日は、融解は起きても凝固が完了しない「不完全サイクル」が頻発する。ヒステリシスがあると凝固開始温度が融解温度より低いから、翌朝になってもPCMが部分的に液体のまま——蓄熱余力がない状態で次の日を迎えてしまう。
対策としては、ヒステリシス対応のツール（TRNSYSやWUFI Plus）を使うか、EnergyPlusなら保守的に見積もる（凝固側のエンタルピー曲線で評価する）という感じですか？
その通り。EnergyPlusを使う場合の実務的な対処法は、DSCの実測から凝固側のエンタルピー曲線を使う（保守的推定）か、融解・凝固の平均曲線を使う（中庸的推定）だ。いずれにせよ、DSCのデータシートに融解・凝固の両曲線が載っているかをメーカーに確認することが第一歩だね。
過冷却の取り扱い
「過冷却」って、PCMが融点以下になっても固まらない現象ですよね？ これもシミュレーションに影響しますか？
パラフィン系PCMは過冷却がほとんどない（1°C以下）から実用上問題にならない。問題になるのは無機塩水和物系（例: CaCl₂·6H₂O、Na₂SO₄·10H₂O）で、3〜10°Cの過冷却を示すことがある。過冷却が大きいと、夜間に室温が融点以下に下がっても凝固が始まらず、蓄えた潜熱を放出できない。
シミュレーション上の対処は、凝固開始温度 $T_\mathrm{nucleation}$ を融点より低く設定すること。エンタルピー曲線を凝固側で $T_\mathrm{nucleation}$ からの放熱として再定義する。ただし核生成は確率的現象なので、決定論的シミュレーションでは完全には再現できない。無機塩系PCMを使う場合は、核生成剤を添加して過冷却を1°C以下に抑制した材料を選定するのが実務的だ。
なるほど。材料選定の段階で過冷却問題を潰しておくのが一番ということですね。シミュレーションで無理やり対処するより、良い材料を選ぶ方が確実だと。
その通り。建築用PCMの実務では「正確なDSCデータを入手する」「メーカーと物性値を擦り合わせる」というアナログな作業が成功の鍵だ。どんなに高度なソルバーを使っても、入力する物性データの質が低ければ garbage in, garbage out。PCMシミュレーションは、材料データの質が結果の信頼性を直接左右する分野なんだ。
なった
詳しく
報告