複合壁の熱伝導
理論と物理
複合壁の基礎理論
複合壁って、異なる材料が積み重なった壁体のことですよね？
そうだ。建築壁体、冷蔵庫の断熱壁、炉壁ライニングなど、実用上の壁体はほぼ全て複合壁だ。各層の熱抵抗を直列接続して全体の熱通過を計算する。
支配方程式
n層の複合壁の全熱抵抗は各層の熱抵抗の和だ。
熱通過率（U値）は
U値って建築の省エネ基準でよく聞く値ですね。
まさにそうだ。日本の省エネ基準（平成28年基準）では、外壁のU値が地域ごとに規定されている。例えば東京（6地域）では外壁 U ≤ 0.53 W/(m2K) だ。
界面温度の算出
各界面の温度は、全体の熱流量 $q$ を求めてから順に計算する。
結露判定にも使えるんですか？
そうだ。各界面温度が露点温度を下回ると結露が発生する。断熱材の位置を変えるだけで結露発生面が変わるので、複合壁の層構成設計は非常に重要だ。
熱抵抗アナロジーの起源
電気回路のオームの法則（1827年）と熱流のアナロジーはFourier（1822年）の熱解析論より後に整備された。1940年代に米国建築学会（ASHRAE）が壁体のU値計算に電気回路的抵抗直列モデルを標準化し、今日の断熱計算の根幹を作った。
各項の物理的意味
- 蓄熱項 $\rho c_p \partial T/\partial t$：単位体積あたりの熱エネルギー蓄積率。【日常の例】鉄のフライパンは熱しにくく冷めにくいが、アルミ鍋は熱しやすく冷めやすい——これは密度 $\rho$ と比熱 $c_p$ の積（熱容量）の違い。熱容量が大きい物体は温度変化が緩やかになる。水は比熱が非常に大きい（4,186 J/(kg·K)）ため、海沿いの気温は内陸より安定する。非定常解析ではこの項が温度の時間変化速度を決める。
- 熱伝導項 $\nabla \cdot (k \nabla T)$：フーリエの法則に基づく熱伝導。温度勾配に比例した熱流束。【日常の例】金属スプーンを熱い鍋に入れると持ち手まで熱くなる——金属は熱伝導率 $k$ が高いため、高温側から低温側へ素早く熱が伝わる。木製スプーンが熱くならないのは $k$ が小さいから。断熱材（グラスウール等）は $k$ が極めて小さく、温度勾配があっても熱が伝わりにくい。「温度差のあるところに熱が流れる」という自然の傾向を数式化したもの。
- 対流項 $\rho c_p \mathbf{u} \cdot \nabla T$：流体の運動に伴う熱輸送。【日常の例】扇風機に当たると涼しく感じるのは、風（流体の流れ）が体表面近くの暖かい空気を運び去り、新鮮な冷たい空気を供給するから——これが強制対流。暖房で部屋の天井付近が暖かくなるのは、暖められた空気が浮力で上昇する自然対流。PCのCPUクーラーのファンも強制対流で放熱している。対流は熱伝導よりも桁違いに効率的な熱輸送手段。
- 熱源項 $Q$：内部発熱（ジュール熱、化学反応熱、放射線吸収等）。単位: W/m³。【日常の例】電子レンジは食品内部のマイクロ波吸収（体積発熱）で加熱する。電気毛布のヒーター線はジュール発熱（$Q = I^2 R / V$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
FEMでの複合壁モデリング
複合壁をFEMでモデル化するポイントは何ですか？
各層を個別の材料領域として定義し、界面でノードを共有させる（またはTied Contact）。重要なのは各層の厚さ方向に十分な要素分割を確保すること。
| 層の厚さ | 推奨要素数（厚さ方向） | 理由 |
|---|---|---|
| 断熱層（低k） | 最低4層 | 温度勾配が急 |
| 構造層（高k） | 2〜3層 | 温度勾配が緩やか |
| 接触層（TIM） | 1〜2層 | 非常に薄い（アスペクト比に注意） |
TIMって厚さ数十umなのに要素を入れるんですか？
TIMを実体メッシュで入れるとアスペクト比が極端に悪化する。代わりにAnsysのThermal Contact（ギャップコンダクタンス）やAbaqusの*GAP CONDUCTANCEで厚さゼロの界面要素として扱う方が効率的だ。
並列熱抵抗
壁体に窓や柱がある場合、並列熱抵抗も考慮する。面積加重平均で
ただしこれは2D/3D効果（熱橋）を無視した近似なので、精密評価にはFEMが必要だ。
鉄骨造の柱が熱橋になるって聞きました。
木造の柱（k=0.15）と鉄骨の柱（k=50）では300倍以上の差がある。鉄骨柱周辺では局所的に温度が下がり結露の原因になる。2D FEMで断面解析するのが標準的な評価手法だ。
境界条件の選び方が鍵
複合壁の定常解析では、表面の対流熱伝達係数hの設定が全体精度を左右する。1980年代のASHRAE 90.1規格では室内側h=8.3 W/m²·Kを規定していたが、2010年版では詳細CFDで計算した局所値の使用を認め、省エネ設計の精度が大幅に向上した。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
設計実務での活用
複合壁の計算は手計算で十分ですか？
1次元の直列抵抗計算は手計算で即座にできるし、Excel化しておくと断熱材の厚さ最適化に便利だ。ただし熱橋の影響評価には2D/3D FEMが必要。
設計例: 冷蔵倉庫の壁体
内部-30度C、外気35度Cの冷蔵倉庫壁体の例:
| 層 | 材料 | 厚さ [mm] | k [W/(mK)] | R [m2K/W] |
|---|---|---|---|---|
| 外壁 | 鋼板 | 0.6 | 50 | 0.000012 |
| 断熱 | ウレタンフォーム | 150 | 0.024 | 6.25 |
| 防湿 | PE膜 | 0.2 | 0.33 | 0.0006 |
| 内壁 | ステンレス | 0.8 | 16 | 0.00005 |
合計R = 6.25 m2K/W（断熱層が支配的）、U = 0.16 W/(m2K)
断熱層のRが99%以上を占めてますね。
そうだ。つまり断熱層以外の層は熱的にはほとんど無視できる。設計のポイントは断熱層の厚さと施工品質（隙間、圧縮）だ。
経年劣化の考慮
断熱材は経年で性能が低下する。ウレタンフォームのkは初期0.024から20年後に0.030程度まで上昇する。設計時には経年劣化係数（通常1.1〜1.3）を見込んでおく。
実際の材料選定では初期性能だけ見てはいけないんですね。
その通り。長期性能保証が重要で、ISO 11561に準拠した加速劣化試験データを確認する。
パッシブハウスの壁厚計算
1991年にドイツのDarmstadtで建設された世界初のパッシブハウスは、複合壁の等価U値を0.1 W/m²·K以下に設計した。これは標準的な日本の外壁（約0.5〜1.0 W/m²·K）の5〜10倍の断熱性能で、年間暖房エネルギーを15 kWh/m²以下に抑えた。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
ツール別対応
複合壁の解析に適したツールを教えてください。
用途別の推奨ツールだ。
| 用途 | ツール | 特徴 |
|---|---|---|
| 建築壁体の省エネ計算 | THERM (LBNL), COMSOL | ISO 10211準拠の熱橋解析 |
| 工業炉壁 | Ansys Mechanical, Abaqus | 高温での非線形材料対応 |
| 電子機器筐体 | Ansys Icepak, FloTHERM | コンポーネントライブラリ |
| 簡易計算 | Excel, Python | 1D熱抵抗ネットワーク |
THERMっていうフリーソフトがあるんですね。
Lawrence Berkeley National Laboratoryが開発した2D定常熱伝導ソルバーだ。建築分野で広く使われており、窓枠の熱橋評価の業界標準になっている。ISO 10077にも準拠している。
Ansys Mechanicalでの設定
複合壁をAnsys Workbenchで組む手順:
1. DesignModeler/SpaceClaimで各層を個別ボディとして作成
2. Shareトポロジーでボディ間のノードを共有
3. 各ボディに材料を割り当て
4. Steady-State Thermal解析タイプを選択
5. 外面にConvection、内面にTemperatureを設定
6. Solve → Temperature分布とTotal Heat Fluxを確認
Shareトポロジーを忘れると界面で断熱になってしまうんですね。
そうだ。Shareトポロジー忘れは最もよくあるミスの一つだ。温度コンターに界面で明確な不連続が見えたら真っ先にここを確認する。
THERM（LBNL）は無償の定番ツール
Lawrence Berkeley National Laboratoryが1994年から無償公開しているTHERMは、複合壁・窓枠の2D熱伝導解析の事実上の国際標準ツールとなっている。2023年時点でも毎月数万件のダウンロードがあり、欧米の建築エネルギー規制（EU EPBD等）への適合確認に広く使われている。
選定で最も重要な3つの問い
- 「何を解くか」：複合壁の熱伝導に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
真空断熱パネル（VIP）
最新の断熱技術にはどんなものがありますか？
真空断熱パネル（VIP）は内部を真空にして気体の熱伝導を排除する。有効熱伝導率は0.004 W/(m K)程度で、従来断熱材の1/6だ。冷蔵庫、建築外壁、LNG配管に使われている。
薄くても高断熱が実現できるんですね。
ただしVIPは釘打ち不可、カット不可という施工上の制約がある。経年で真空度が低下しkが上昇するリスクもある。解析ではkの経年変化を考慮した感度解析が重要だ。
エアロゲル断熱
エアロゲル（k=0.013〜0.018 W/(m K)）は薄膜でも高断熱を実現でき、配管保温やEVバッテリーの断熱に急速に普及している。Ansys Fluentの多孔質媒体モデルでエアロゲル内の熱輸送を解析できる。
複合壁の最適設計
各層の厚さと材料の組み合わせを最適化する問題は、制約付き最適化として定式化できる。
COMSOLのOptimization ModuleやAnsys DesignXplorerのDirect Optimizationで自動化できる。
コストと性能のトレードオフを自動で最適化できるんですね。
DOE（Design of Experiments）+ 応答曲面法が実用的だ。断熱厚さ、材料種類、層順序の組み合わせを系統的に評価し、コスト最小の設計を見つける。
熱橋（サーマルブリッジ）の3D効果
2次元以上の熱橋補正係数ψ（プサイ値）は1990年代にISO 10211で標準化された。スチールスタッドが入る壁では1Dモデルに比べ実効U値が40%以上増大することが多く、3D FEM解析なしに正確な断熱評価はできない。Schöckのアイソコルブのような熱橋対策製品の市場もこの認識から生まれた。
トラブルシューティング
よくある問題
複合壁の解析でありがちなミスを教えてください。
1. 層間接触の設定漏れ
問題: 各層がボディとして独立しており、Share Topologyやボンド接触が未設定だと、層間が断熱になる。温度分布が不連続になったり、一方の層の温度が変化しない。
対策: Ansys Workbenchなら「Connections > Contact Region > Thermal Conductance」を設定。完全接触なら極めて大きな値（例: 10^6 W/(m2K)）、実際のTIMなら測定値を入力。
2. 極薄層のメッシュ問題
TIMの厚さ50umを実体メッシュで入れたらエラーになりました。
原因: 隣接層との厚さ比が1:1000以上になるとアスペクト比が破綻する。
対策: 薄層はInterface（Shell Conduction）かContact Conductanceで表現する。Ansysでは「Thermal Contact」の「Conductance Value」= k/t で設定する。例: k=3 W/(mK), t=50um なら Conductance = 60000 W/(m2K)。
3. 放射の見落とし
高温炉壁（500度C以上）では壁面間の放射が無視できない。複合壁内部のエアギャップがある場合、空気の伝導だけでなく放射熱伝達も考慮する。
放射の分をh_radとして対流項に加算するんですね。
1Dモデルではそうする。3D FEMでは面間のS2S放射モデルを入れる。Ansys MechanicalのRadiation BoundaryまたはAbaqusの*RADIATIONキーワードで設定する。
接着剤層の見落としで温度分布が狂う
実験棟の複合壁解析でシミュレーション値と実測が15℃ずれた事例がある。原因は厚さ0.3 mmの接着剤層（λ≈0.2 W/m·K）の入力漏れで、薄いが熱抵抗への寄与は無視できなかった。CAEモデルでは膜厚が小さい層ほど意識的にチェックリストに加える運用が重要。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——複合壁の熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告