フーリエの法則
理論と物理
フーリエの法則とは
先生、フーリエの法則って伝熱解析の一丁目一番地って聞いたんですけど、何がそんなに大事なんですか？
フーリエの法則は「温度勾配があるところに熱が流れる」という自然現象を数式にしたものだ。1822年にJoseph Fourierが著書『熱の解析的理論』で定式化した。全ての熱伝導解析はここから始まる。
200年前の法則が今でも現役ってすごいですね。
そう、ニュートン力学がF=maで始まるように、伝熱工学はフーリエの法則で始まる。ベクトル形式で書くとこうなる。
支配方程式
フーリエの法則の基本形は熱流束ベクトル $\mathbf{q}$ と温度勾配の関係式だ。
ここで $k$ は熱伝導率 [W/(m K)]、$\nabla T$ は温度勾配 [K/m] だ。マイナス符号は「熱は高温側から低温側に流れる」という熱力学第二法則を反映している。
1次元だともっとシンプルになりますよね？
そうだ。1次元の場合は偏微分が常微分になる。
これをエネルギー保存則と組み合わせると、定常熱伝導方程式が得られる。
$\dot{q}_v$ は単位体積あたりの内部発熱量 [W/m3] だ。kが一定なら $k \frac{d^2T}{dx^2} + \dot{q}_v = 0$ とさらに簡単になる。
内部発熱がない場合は温度分布が直線になるってことですか？
その通り。平板で両面の温度が $T_1$、$T_2$ なら $T(x) = T_1 + (T_2 - T_1)\frac{x}{L}$ という線形分布になる。これは解析解の検証でよく使う基本ケースだ。
3次元への拡張
一般的な3次元問題では、異方性材料を考慮してテンソル形式で書く。
等方性材料なら $k_{ij} = k \delta_{ij}$ で、ラプラス方程式 $\nabla^2 T = 0$（発熱なし）またはポアソン方程式 $\nabla^2 T + \frac{\dot{q}_v}{k} = 0$（発熱あり）に帰着する。
CFRPみたいな複合材料だと、繊維方向と直交方向でkが全然違いますよね？
そうだ。CFRPの場合、繊維方向の熱伝導率が5〜10 W/(m K) に対して直交方向は0.5〜1 W/(m K) 程度。テンソル形式でないと正しく扱えない。
境界条件
境界条件は3種類ある。
境界条件は3種類ある。
| 種類 | 名称 | 数式 | 物理的意味 |
|---|---|---|---|
| 第1種 | Dirichlet | $T = T_s$ | 表面温度を指定 |
| 第2種 | Neumann | $-k \frac{\partial T}{\partial n} = q_s$ | 熱流束を指定 |
| 第3種 | Robin | $-k \frac{\partial T}{\partial n} = h(T - T_\infty)$ | 対流熱伝達 |
断熱条件はNeumann条件で $q_s = 0$ にするってことですね。
そうだ。実務では第3種境界条件が最も頻繁に使われる。対流熱伝達係数 $h$ の見積もりが結果の精度を大きく左右するんだ。
フーリエの法則が生まれた年
ジョゼフ・フーリエは1822年の著書『熱の解析的理論』でq = −k∇Tを提唱した。この式は伝熱工学の出発点であり、200年以上経った現代のFEAソルバーでも基本支配方程式として変わらず使われている。
各項の物理的意味
- 蓄熱項 $\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$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
有限要素法による離散化
フーリエの法則をコンピュータで解くにはどうするんですか？
定常熱伝導問題の弱形式（ガラーキン法）から出発する。温度場 $T$ を形状関数 $N_i$ で近似する。
弱形式にして部分積分を適用すると、要素レベルの方程式が得られる。
ここで $K^e_{ij} = \int_{\Omega_e} k \nabla N_i \cdot \nabla N_j \, d\Omega$ が要素熱伝導マトリクス、$f^e_i = \int_{\Omega_e} \dot{q}_v N_i \, d\Omega + \int_{\Gamma_e} q_s N_i \, d\Gamma$ が要素熱負荷ベクトルだ。
構造解析の剛性マトリクスと形が似てますね。
いい着眼点だ。構造の $[K]\{u\}=\{F\}$ と数学的に同じ構造だ。ただし熱解析は未知数がスカラーの温度なので、節点あたりの自由度は1つだけ。構造解析より問題サイズがずっと小さくなる。
有限差分法・有限体積法
有限差分法（FDM）も定常熱伝導でよく使われる。中心差分で離散化すると、1次元では
有限体積法（FVM）はセル中心に温度を配置し、セル界面での熱流束を保存する。CFDソルバー（Ansys Fluent, STAR-CCM+）はFVMベースなので、共役熱伝達解析では自動的にFVMで固体側も離散化される。
方法によって精度は変わりますか？
FEMは複雑形状に強く高次要素で高精度が得られる。FVMは保存則を厳密に満たすので流体と固体の連成に向く。FDMは実装が容易だが構造格子に限定される。使い分けが重要だ。
行列解法
組み上がった全体方程式 $[K]\{T\}=\{f\}$ は、定常熱伝導の場合は線形なので一度の解法で済む。
| 解法 | 種類 | 特徴 | 推奨規模 |
|---|---|---|---|
| Cholesky分解 | 直接法 | 対称正定値に最適、正確 | 〜50万DOF |
| PCG法 | 反復法 | メモリ効率良好 | 50万〜1000万DOF |
| AMG前処理+CG | 反復法 | 大規模問題に強い | 1000万DOF以上 |
温度依存の熱伝導率がある場合はどうなりますか？
非線形問題になるので反復計算が必要だ。Newton-Raphson法で $k(T)$ を逐次更新する。収束判定は残差ノルム $10^{-6}$ 程度が標準だが、温度変化量で $10^{-3}$ K を基準にすることもある。
熱伝導率の測定手法
フラッシュ法（レーザーフラッシュ法）はASTM E1461規格に基づき、0.1〜2000 W/m·Kの広範囲を±3%精度で測定できる。ネッチュ社のLFA467装置は1回の試験で熱拡散率・比熱・熱伝導率を同時取得できる標準機として世界中で使われる。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
解析フロー
実際にフーリエの法則を使った熱伝導解析を回すとき、何から始めればいいですか？
定常熱伝導解析の標準フローはこうだ。
1. 形状簡略化: 対称面の活用、微小フィレットの除去、薄肉部のシェル化を検討
2. 材料定義: 熱伝導率 $k$ [W/(m K)] を温度依存テーブルで入力
3. メッシュ作成: 温度勾配の大きい領域を細かく、遠方を粗く
5. 求解・後処理: 温度分布と熱流束分布を確認
メッシュの粗さってどのくらいが目安ですか？
まず粗いメッシュで全体傾向を掴み、注目領域を2回細分化して結果が1%以内で変わらなければ収束したと判断する。定常熱伝導は構造解析ほどメッシュ感度が高くないので、比較的粗いメッシュでも妥当な結果が得られることが多い。
代表的な材料の熱伝導率
| 材料 | $k$ [W/(m K)] | 用途例 |
|---|---|---|
| 銅 (C1100) | 398 | ヒートスプレッダ |
| アルミ (A6063) | 200 | ヒートシンク |
| SUS304 | 16.3 | 構造部材 |
| エポキシ (FR-4) | 0.3 | PCB基板 |
| 空気 (25度C) | 0.026 | 自然対流間隙 |
| サーマルグリス | 1〜5 | TIM界面 |
銅とエポキシで1000倍以上差があるんですね。
だからこそ電子基板の熱設計ではビアやサーマルパッドで熱経路を確保することが重要なんだ。熱抵抗ネットワークで $R = L/(kA)$ を計算し、ボトルネックを見極める。
結果の検証方法
解析結果の妥当性は以下で確認する。
- エネルギー収支: 流入熱量と流出熱量が一致するか（誤差1%以内）
- 理論解との比較: 単純形状部分を理論解 $q = kA\Delta T/L$ で検算
- 温度の物理的妥当性: 境界条件の範囲内に温度が収まっているか
- 熱流束ベクトル: 断熱面に垂直成分がないか、対称面で条件を満たすか
エネルギー収支って具体的にどう確認するんですか？
Ansys Mechanicalなら「Heat Flow」のReaction Summaryで確認できる。全ての温度拘束面と対流面での熱量合計がゼロに近ければOKだ。COMSOLなら「Surface Integration」で面ごとの熱流量を計算する。
断熱材設計への応用
航空機の胴体断熱にはグラスウール（k≈0.04 W/m·K）を厚さ50mm充填し、−55°Cの外気から客室22°Cを維持する。フーリエの法則q=kA(ΔT/L)を用いると1m²あたり熱損失は約28Wと計算でき、燃料消費予測に直結する。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
商用ツールでの実装
フーリエの法則に基づく定常熱伝導をやるには、どのソフトがいいですか？
定常熱伝導はほぼ全ての汎用FEMソルバーで対応している。代表的なツールを比較しよう。
| ツール | 定常熱伝導の特徴 | 要素タイプ |
|---|---|---|
| Ansys Mechanical | SOLID70(8節点六面体)、SOLID87(10節点四面体)。APDLマクロでパラメトリック解析が容易 | FEM |
| Abaqus | DC3D8(六面体), DC3D4(四面体)。*HEAT TRANSFER, STEADY STATEステップで定義 | FEM |
| COMSOL | Heat Transfer in Solidsモジュール。GUI上で方程式をカスタマイズ可能 | FEM |
| Ansys Fluent | 固体領域もFVMで離散化。CHT解析の固体側に自動適用 | FVM |
Ansys MechanicalとFluent、同じAnsysなのに使い分けるんですね。
固体のみの熱伝導ならMechanicalが効率的だ。流体との連成があるならFluentやCFXで共役熱伝達（CHT）として解く方がワークフローがシンプルになる。
APDL実装例
Ansys APDLで平板の定常熱伝導を解く最小コードはこうだ。
```
/PREP7
ET,1,SOLID70
MP,KXX,1,200 ! アルミ k=200 W/(mK)
BLOCK,0,0.1,0,0.05,0,0.01
ESIZE,0.002
VMESH,ALL
/SOL
ANTYPE,STATIC
D,NODE(0,,,),,100 ! 左面 100℃
SF,AREA(0.1,,,),CONV,10,25 ! 右面 h=10, T∞=25℃
SOLVE
```
思ったよりシンプルですね。Abaqusだとどうなりますか？
Abaqusではinpファイルに STEP セクションで HEAT TRANSFER, STEADY STATE を指定する。境界条件は BOUNDARY で温度固定、FILM で対流条件を与える。
ツール選定の指針
選定基準をまとめるとこうなる。
| 判断基準 | 推奨ツール |
|---|---|
| 固体のみ・構造連成あり | Ansys Mechanical, Abaqus |
| 流体連成（CHT） | Ansys Fluent, STAR-CCM+ |
| マルチフィジックス | COMSOL |
| 電子機器専用 | Ansys Icepak, FloTHERM |
| 教育・研究 | COMSOL, OpenFOAM |
電子機器系だとIcepakやFloTHERMが出てくるんですね。
IcepakやFloTHERMはコンポーネントライブラリが充実していて、ファンカーブやヒートシンクのパラメトリック解析に特化している。内部的にはフーリエの法則+対流モデルで解いている点は同じだ。
フーリエ則実装のFEMソルバー
Abaqus 2024では熱伝導率テンソルの全成分を温度依存テーブルとして定義でき、CFRP積層板の異方性定常解析を高精度に実行できる。COMSOL Multiphysicsは材料ライブラリに6,000種超の熱物性値を収録し研究用途に広く使われる。
選定で最も重要な3つの問い
- 「何を解くか」：フーリエの法則に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
フーリエの法則の限界
フーリエの法則が使えない場面ってあるんですか？
ある。フーリエの法則は連続体・拡散的な熱輸送を仮定している。ナノスケールやフェムト秒レーザー加工では、フォノンの平均自由行程と比較して構造サイズが小さくなり、弾道的熱輸送が支配的になる。
そういう場合はどうするんですか？
ボルツマン輸送方程式（BTE）や分子動力学（MD）を使う。半導体デバイスのナノスケール熱解析ではBTEベースのツール（たとえばMITのAlmaBTE）が研究レベルで使われている。
温度依存性と非線形問題
多くの材料で $k$ は温度に依存する。たとえばアルミA6063は25度Cで200 W/(m K)だが、300度Cでは220 W/(m K) に上昇する。一方、SUS304は25度Cで16.3だが300度Cで21.5になる。
この温度依存性を考慮すると非線形問題になり、Newton-Raphson反復が必要だ。Ansys Mechanicalでは材料プロパティのテーブル入力（MPTEMP/MPDATA）で自動的に非線形ソルバーが起動する。
トポロジー最適化への応用
最近の研究トレンドってどんなのがありますか？
フーリエの法則に基づく熱伝導トポロジー最適化が注目されている。目的関数を「最大温度の最小化」や「熱コンプライアンスの最小化」として、材料分布を最適化する。
$\rho$ は各要素の材料密度（0〜1）だ。Ansys MechanicalやCOMSOLのTopology Optimizationモジュールで実装できる。アディティブマニュファクチャリングと組み合わせて最適な放熱経路を実現する事例が増えている。
フォノニック結晶と熱メタマテリアル
もう一つの先端テーマが、フォノニック結晶を使った指向性熱伝導だ。周期構造でフォノンのバンドギャップを作り、特定方向への熱流を制御する。これはフーリエの法則が成立するマクロスケールでも有効テンソルの異方性設計として応用できる。
熱を「操る」って感じですね。SF感がある。
まだ研究段階だが、将来的には電子デバイスの熱管理や宇宙機器の断熱設計に応用が期待されている。
フォノンと熱伝導の量子論
ナノスケールではフーリエの法則が破綻する。シリコンの熱伝導率はバルクの148 W/m·Kから膜厚10nmでは約20 W/m·Kまで低下する。2001年にChenが提唱したボルツマン輸送方程式（BTE）ベースの解析がIntelの3D積層半導体設計で活用されている。
トラブルシューティング
よくあるトラブルと対策
定常熱伝導解析で困ったとき、まず何を見ればいいですか？
よくある問題を整理しよう。
1. 温度が発散する・非現実的な値になる
原因: 境界条件の不足。定常熱伝導には最低1つの温度拘束（Dirichlet条件）または対流条件（Robin条件）が必要だ。全面にNeumann条件（熱流束）だけ指定すると、温度の基準が決まらず解が一意に定まらない。
対策: モデルのどこかに温度の基準点を設ける。物理的に温度が既知の面を見直す。
それは構造解析でいう剛体変位の問題と同じですね。
まさにそうだ。構造で拘束なしだと剛体移動するのと同様に、熱で温度拘束なしだと温度が浮動する。
2. エネルギー収支が合わない
原因: メッシュ不整合、接触面のギャップ、材料定義の誤り。
対策: 反力の合計値（Ansysなら FSUM, Abaqusなら Output > Energy）を確認。流入熱量と流出熱量の差が全体の1%未満であることを検証する。大きなずれがある場合、ノードのマージ漏れやTied Contactの設定ミスを疑う。
3. 温度分布が予想と大きく異なる
解析は回ったけど結果が実測と合わないときは？
チェックリストはこうだ。
| チェック項目 | よくあるミス |
|---|---|
| 単位系 | W/(m K) と W/(mm K) の混同（1000倍の差） |
| 熱伝達係数h | 自然対流で5〜10 W/(m2K) のところに100を入れている |
| 接触熱抵抗 | TIMやボルト締結部のギャップコンダクタンス未設定 |
| 材料テーブル | 温度依存テーブルの範囲外で外挿されている |
| 放射 | 高温部品で放射を無視している |
特に単位系の混同は致命的だ。Ansys MechanicalはSI単位系（m, kg, s）とmm単位系（mm, t, s）でkの値が変わる。mm系ではk=200がmW/(mm K)=0.2 W/(mm K)になる点に注意。
4. メッシュ収束しない
定常熱伝導で温度のメッシュ収束が悪い場合、特異点の近傍（角部、材料不連続点）を疑う。熱流束は温度勾配の微分なので、温度が収束していても熱流束は収束しにくい。
実務的にはどうすればいいですか？
温度そのものに着目すれば十分なケースが多い。熱流束の局所値が必要な場合は、パスに沿った分布を確認し、特異点から離れた位置で評価する。
異方性材料でのフーリエ則適用
炭素繊維強化プラスチック（CFRP）は繊維方向k≈7 W/m·K、垂直方向k≈0.7 W/m·Kと10倍の異方性を持つ。等方性を仮定して解析すると温度分布の最大誤差が30%に達することがあり、テンソル形式の熱伝導率行列の入力が必須だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——フーリエの法則の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告