1次元定常熱伝導
理論と物理
1次元定常熱伝導の基礎
先生、1次元の定常熱伝導って、実際どんな場面で使うんですか？
平板・円筒・球殻のように一つの方向にだけ温度勾配がある問題に適用する。壁体の断熱評価、配管保温材の設計、電線の許容電流計算などが典型例だ。
支配方程式
1次元定常・内部発熱ありの支配方程式はこうなる。
$k$ が定数で発熱がなければ $\frac{d^2T}{dx^2} = 0$ となり、温度は直線分布になる。
一様発熱があるとどうなりますか？
両面温度固定 $T(0)=T_1$, $T(L)=T_2$ で $\dot{q}_v$ 一様なら
二次曲線が重畳される。最高温度は必ずしも中央ではなく、$T_1 \neq T_2$ の場合は偏る。
熱抵抗の概念
1次元熱伝導では電気回路のアナロジーが非常に有効だ。平板の熱抵抗は
対流の熱抵抗は $R_{conv} = \frac{1}{hA}$ だ。これを直列・並列に接続して全体の温度降下を見積もる。
まさにオームの法則ですね。$\Delta T = qR$ がV=IRに対応する。
その通り。この熱抵抗ネットワークの考え方は、FloTHERMのCompact Thermal Modelや、JEDECのDELPHIモデルの基盤にもなっている。
Fourierの熱方程式、獄中で着想
Joseph Fourier（1768–1830）は1798年のエジプト遠征から帰国後、投獄期間中にも研究を続けた。1822年刊行の『熱の解析的理論』で1次元熱伝導方程式を完成させたが、王立アカデミーは初稿（1807年）を「厳密性が欠ける」として12年間却下し続けた。
各項の物理的意味
- 蓄熱項 $\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$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
有限差分法による離散化
1次元なら手計算で解けそうですけど、あえて数値解法を使う理由は？
温度依存の熱伝導率 $k(T)$ や非一様発熱がある場合は解析解が得られない。FDMで等間隔格子に離散化すると
ここで $k_{i+1/2}$ はセル界面の調和平均で評価する。$k_{i+1/2} = \frac{2k_i k_{i+1}}{k_i + k_{i+1}}$
調和平均を使うのはなぜですか？
異種材料の界面で熱流束の連続性を保つためだ。算術平均だと界面で熱流束が不連続になりうる。これはCFDソルバーでも同じ取り扱いがされている。
FEMによる1次元解析
1次元FEMでは2節点線形要素を使う。要素熱伝導マトリクスは
構造のバネ要素と全く同じ形だ。対流境界は追加項 $hA\begin{bmatrix}0&0\\0&1\end{bmatrix}$ を右端に加える。
1次元だとExcelでも実装できそうですね。
そうだ。教育目的には10要素程度のExcel実装が非常に有効で、FEMの本質を理解できる。実務ではもちろん汎用ソルバーを使うが、検証用の自作ツールとして手元に持っておくと役立つ。
有限差分法の原点はRichardson
Lewis Fry Richardson（1910年）が流体方程式を差分近似した手法が1次元熱伝導の数値解法の原型だ。第一次世界大戦中、気象予報の手計算に差分格子を使い、64人の「人間コンピュータ」を動員する壮大な実験を行ったことで有名。計算精度の概念を数値熱力学に持ち込んだ先駆者。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
設計計算での活用
1次元モデルって実務でも本当に使いますか？3Dがあるのに。
概念設計段階では1次元モデルが圧倒的に効率的だ。壁体の断熱厚さ決定、配管保温材の選定、電線サイズの決定は1次元計算で十分な精度が出る。
実務例: 配管保温材の設計
外径50mmの蒸気配管（150度C）にグラスウール保温材（$k=0.04$ W/(m K)）を巻く場合、外気25度C、外表面熱伝達係数 $h=10$ W/(m2K) として
保温厚さ50mmなら $r_i=25$mm, $r_o=75$mm で、単位長さ当たり熱損失は約20 W/m になる。
これ、手計算で10秒で出ますね。3Dシミュレーションを回すまでもない。
そうだ。ただし配管エルボや分岐部、フランジ部は2次元・3次元効果が出るので、全体の熱損失量計算には1D、局所の温度評価には3Dと使い分ける。
結果の検証
3D解析の結果検証に1D理論解を使うのは非常に有効だ。
| 検証項目 | 1D理論値 | 3D解析値 | 判定基準 |
|---|---|---|---|
| 最高温度 | 理論式から算出 | ソルバー出力 | 差が5%以内 |
| 熱流量 | $q=kA\Delta T/L$ | 面積分 | 差が2%以内 |
| 温度勾配 | $dT/dx = -q/(kA)$ | パスプロット | 分布の一致 |
1Dが分かっていると3Dの結果をすぐに検証できるんですね。
「桁が合っているか」の確認だけでも、設計ミスの80%は防げる。これが1次元モデルの最大の価値だ。
炉壁設計の黄金律
製鉄高炉の炉壁は1次元定常近似で設計されることが多い。新日鐵住金（現日本製鉄）の高炉では耐火レンガ→断熱キャスタブル→鉄皮の3層構成を採用し、内側1600℃・外側60℃を維持するためにλと厚さを最適化している。炉壁単位面積あたりの熱損失低減が年間エネルギーコストに直結する。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
ツール別の取り扱い
1次元熱伝導を各ソフトでやる方法を教えてください。
汎用FEMソルバーでの実装方法を比較しよう。
| ツール | 要素タイプ | 設定のポイント |
|---|---|---|
| Ansys Mechanical | LINK33（2節点熱伝導バー）、SOLID70の1D使用 | 断面積をREALで定義 |
| Abaqus | DC1D2（2節点線形）、DC1D3（3節点二次） | *SOLID SECTION で断面積指定 |
| COMSOL | Heat Transfer in Solids → 1D component | 幾何的に1D線を作成 |
電子機器の熱抵抗ネットワークだと別のツールを使いますよね？
回路シミュレータ的なアプローチなら、FloTHERMのSmartPartやAnsys Icepakのネットワークモデル、あるいは6SigmaETの簡易モデルが適している。これらは内部的にはR-Cネットワーク（定常ではRのみ）で1次元熱抵抗を並べている。
JEDEC標準との関連
IC部品の熱特性はJEDEC JESD51シリーズで規定されている。$\theta_{JA}$（ジャンクション-環境間熱抵抗）や $\Psi_{JT}$（ジャンクション-ケーストップ間）は本質的に1次元熱抵抗モデルだ。
データシートに書いてある $\theta_{JA}$ はそういう意味だったんですね。
ただし $\theta_{JA}$ は測定条件（JEDEC標準基板、自然対流）に依存する。実基板の実装条件とは異なるので、実設計では3Dシミュレーションが必須だ。1Dモデルはあくまでも初期見積もりと位置づける。
MATLABのpdepe関数の誕生
MathWorksはMATLAB 5.2（1997年）でpdepe関数を搭載し、1次元放物型・楕円型PDEを数行コードで解けるようにした。それ以前はユーザーがCrank-Nicolson法を手実装するのが一般的で、境界条件のコーディングミスによるバグが横行していた。pdepeの登場は1D熱伝導教育の敷居を大幅に下げた。
選定で最も重要な3つの問い
- 「何を解くか」：1次元定常熱伝導に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
非線形1次元問題
1次元でも先端的な話題ってあるんですか？
温度依存熱伝導率のKirchhoff変換が重要なテクニックだ。変数変換
を施すと、非線形方程式が線形のラプラス方程式 $\frac{d^2\theta}{dx^2}=0$ に帰着する。
非線形問題が線形になるのはすごいですね。
半導体デバイスの熱設計で、シリコンの $k(T) = k_0(T/T_0)^{-1.3}$ のような強い温度依存性がある場合に特に有効だ。
機能性傾斜材料（FGM）
熱伝導率が空間的に連続変化するFGMも1次元で基礎的な知見が得られる。
このような指数分布を仮定すると解析解が得られ、TBCコーティング最適化の指針になる。
逆問題：熱伝導率の同定
温度測定データから熱伝導率を逆算するのって1次元でもできますか？
可能だ。2点の温度測定値と熱流束測定値があれば、$k = qL/(A\Delta T)$ で直接算出できる。温度依存kの同定には、複数温度での測定と最適化アルゴリズム（Levenberg-Marquardt法）を組み合わせる。
COMSOLのOptimization Moduleでは、実測温度との差を目的関数として $k(T)$ のパラメータを自動同定できる。1次元問題なら計算が軽いので、反復最適化も高速に回せる。
Coffee Break よもやま話
薄膜の量子熱伝導
1次元の古典Fourier則は膜厚がフォノン平均自由行程（Siで約300 nm@室温）を下回ると破綻する。2001年にDavid Cahillらが発表した時間領域サーモリフレクタンス（TDTR）法は、10 nm以下の薄膜でフォノン輸送を実測する手法を確立し、ナノエレクトロニクスの熱設計に革命をもたらした。
温度測定データから熱伝導率を逆算するのって1次元でもできますか？
可能だ。2点の温度測定値と熱流束測定値があれば、$k = qL/(A\Delta T)$ で直接算出できる。温度依存kの同定には、複数温度での測定と最適化アルゴリズム（Levenberg-Marquardt法）を組み合わせる。
COMSOLのOptimization Moduleでは、実測温度との差を目的関数として $k(T)$ のパラメータを自動同定できる。1次元問題なら計算が軽いので、反復最適化も高速に回せる。
薄膜の量子熱伝導
1次元の古典Fourier則は膜厚がフォノン平均自由行程（Siで約300 nm@室温）を下回ると破綻する。2001年にDavid Cahillらが発表した時間領域サーモリフレクタンス（TDTR）法は、10 nm以下の薄膜でフォノン輸送を実測する手法を確立し、ナノエレクトロニクスの熱設計に革命をもたらした。
トラブルシューティング
よくある問題
1次元なら問題は少なそうですけど、それでもハマるポイントはありますか？
1次元特有の落とし穴がある。
1. 座標系の選択ミス
円筒座標系の1次元問題で直交座標系の式を使ってしまうミスが多い。円筒座標では
であり、$r$ の項を落とすと温度分布が全く変わる。平板では線形分布だが円筒では対数分布 $T(r) = C_1 \ln r + C_2$ だ。
球座標でも同じような注意が必要ですか？
そうだ。球座標では $\frac{1}{r^2}\frac{d}{dr}(kr^2\frac{dT}{dr})=0$ で、解は $T(r) = C_1/r + C_2$ になる。座標系を間違えると温度降下の見積もりが大きくずれる。
2. 複合壁の界面温度
複合壁で界面温度を手計算するとき、各層の熱抵抗の分配を間違えるケースが多い。
全体の $q$ を先に求めてから各界面温度を順次計算する手順を守れば間違いにくい。
3. 3Dモデルとの不一致
1Dの見積もりと3Dの結果が大きく違うときは何を疑えばいいですか？
以下を確認する。
| 不一致要因 | 確認方法 |
|---|---|
| 2D/3D的な広がり効果 | ヒートスプレッディングの寄与を推算 |
| 接触熱抵抗 | 界面のギャップコンダクタンスを確認 |
| 放射の影響 | 高温部で放射が無視できない |
| 対流係数の見積もり | 1Dで仮定したhが3D流れ場と乖離 |
ヒートスプレッディング効果は特に見落としやすい。小さな発熱源から広い平板に熱が広がる場合、1Dモデルでは捉えられない2D/3D効果で実際の温度上昇は1Dの見積もりより低くなる。Song-Lee-Au のスプレッディング抵抗の近似式で補正するとよい。
熱電対位置のずれが誤差を生む
1次元定常解析の検証実験で、熱電対先端が壁面から0.5 mmずれるだけで測定値が3〜8℃変わることが多い。NELの不確かさ評価ガイド（NPL Measurement Note）では、熱電対接合部の位置決め誤差を系統誤差の主要因として挙げており、実験との対比には取付精度の詳細記録が必須とされる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——1次元定常熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告