2次元定常熱伝導
理論と物理
2次元熱伝導の基礎
2次元の定常熱伝導って、1次元からどう拡張されるんですか？
温度が2方向に変化する問題だ。電子基板の面内温度分布、金型の断面温度場、建築壁体の熱橋評価などが代表例だ。
支配方程式
等方性材料の2次元定常熱伝導方程式はこうなる。
$k$ が一定で発熱がなければラプラス方程式 $\nabla^2 T = 0$ に帰着する。
解析解はあるんですか？
矩形領域で辺ごとに温度を指定する場合は、変数分離法で級数解が得られる。例えば3辺が0度C、1辺が $T_0$ の正方形領域では
この解はFEMコードの検証ベンチマークとして非常に重要だ。
形状係数法
2次元の幾何学的効果を1次元的に扱う便法として形状係数 $S$ がある。
例えば地中埋設管の熱損失は $S = 2\pi L / \cosh^{-1}(D/r)$ で評価できる。教科書（Incropera等）に主要形状の $S$ がまとめられている。
形状係数を使えば、2Dの問題も手計算できるんですね。
そうだ。ただし複雑形状では形状係数が未知なので、FEMで解くことになる。形状係数の値をFEMの結果から逆算して、類似設計に再利用することも実務ではよくある。
Laplace方程式と解析解の美しさ
2次元定常熱伝導は∇²T=0（Laplace方程式）で記述される。1822年にFourierが示した矩形板の解析解はsin級数で表され、今日でも検証ベンチマークとして使われる。Laplace方程式は電位・流速ポテンシャルとも同形で、19世紀の物理学者はこの統一性を「自然の深い美しさ」と称賛した。
各項の物理的意味
- 蓄熱項 $\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による2次元離散化
2次元問題のFEMではどんな要素を使いますか？
2次元熱伝導には三角形要素と四角形要素がある。
| 要素 | 節点数 | 温度分布 | 精度 |
|---|---|---|---|
| 3節点三角形 | 3 | 線形 | 低（温度勾配が一定） |
| 6節点三角形 | 6 | 二次 | 高 |
| 4節点四角形 | 4 | 双線形 | 中 |
| 8節点四角形 | 8 | 二次 | 高 |
3節点三角形は温度勾配が要素内で一定なので、細分化しないと精度が出ない。実務では6節点三角形か8節点四角形が推奨だ。
四角形と三角形のどっちがいいですか？
同じ節点数なら四角形の方が精度が高い。しかし複雑形状の自動メッシュでは三角形（テトラ）の方が生成しやすい。Ansysのメッシャーは六面体優先（hex-dominant）を選ぶとバランスが良い。
有限差分法（2次元）
等間隔格子での中心差分は
$\Delta x = \Delta y$ なら有名な5点ステンシルになる。Excelや自作Pythonコードで検証用に実装できる。
商用ソルバーの結果を自作コードで検証するっていうアプローチ、堅実ですね。
特に2次元の検証はExcel VBAやPythonの数十行で書けるので、やらない理由がない。
FDMの中央差分は2次精度
2次元FDMの中央差分スキームは空間2次精度で、格子幅を半分にすると誤差が1/4になる。1960年代のIBM 704上でNASAが実施したロケットノズル2D熱解析では、128×128格子でも数時間を要したが、現代のPCなら0.01秒以下で同じ計算が終わる。50年間で演算速度は10⁹倍以上向上した。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
2Dモデルの活用場面
3D解析が主流なのに、わざわざ2Dモデルを使う利点は何ですか？
計算コストが桁違いに小さい。パラメトリックスタディで数十ケース回す場面や、断面方向の温度勾配が支配的な問題では2Dが最適だ。
実務例: PCB基板の熱橋解析
多層PCBの断面をモデル化する場合、各層（銅箔、プリプレグ、コア）を異なる $k$ を持つ帯として並べる。銅箔率が80%のL1層と20%のL3層では等価熱伝導率が大きく異なる。
| 層 | 厚さ [um] | 面内k [W/(mK)] | 面直k [W/(mK)] |
|---|---|---|---|
| L1 (銅80%) | 35 | 318 | 1.1 |
| プリプレグ | 100 | 0.3 | 0.3 |
| L2 (銅50%) | 35 | 199 | 0.55 |
| コア | 800 | 0.3 | 0.3 |
面内と面直で2桁以上違うんですね。
PCB基板は極端な異方性材料だ。2D断面解析で面直方向の温度勾配とビアの効果を評価し、面内方向はスプレッディング抵抗で補正する。これがFloTHERMやIcepakの内部モデリングでも使われる手法だ。
対称性の活用
対称性を利用すれば計算規模を大幅に削減できる。
- 1軸対称: 1/2モデル、対称面に断熱条件
- 2軸対称: 1/4モデル
- 周期対称: 繰り返し単位を1つだけモデル化
対称面が断熱条件になる理由は？
対称面を挟んで温度場が鏡像対称なら、面を横切る温度勾配がゼロ、すなわち熱流束がゼロだからだ。物理的に考えれば当然の帰結だ。
CPUダイの2D温度マップ
Intelは2003年のPentium 4「Prescott」コア（90 nmプロセス、TDP 84 W）で2D定常熱解析を駆使してホットスポット管理を行った。ダイ上の最大温度差が20℃以上になり、電流密度集中箇所を2Dマップで特定して配線層を再設計した。この手法はその後のマルチコア設計の標準となった。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
各ツールでの2D解析
2D熱伝導を各ソフトでやるにはどうしますか？
主要ツールの2D熱解析要素を一覧にする。
| ツール | 要素 | 備考 |
|---|---|---|
| Ansys Mechanical | PLANE55(4節点), PLANE77(8節点) | 平面/軸対称を選択 |
| Abaqus | DC2D4(四角形), DC2D3(三角形) | DCAX4で軸対称 |
| COMSOL | 2D component → Heat Transfer | 任意形状に自動メッシュ |
| Ansys Fluent | 2D/Axi-symmetric mesh | FVMベース |
軸対称モデルは2Dで3Dを表現できるから便利ですね。
その通り。円筒形ヒートシンクや回転体の温度場は軸対称2Dで十分な精度が得られる。3Dの1/100以下の計算コストで同等の結果を得られる場合も多い。
COMSOLでの設定例
COMSOLで2D定常熱伝導を組む手順:
1. Model Wizard → 2D → Heat Transfer in Solids
2. Geometry: Rectangle等で領域作成
3. Materials: Built-in materialまたはUser defined
4. Heat Transfer: Thermal insulation(デフォルト), Temperature, Heat Flux等を設定
5. Mesh → Physics-controlled (Normal)
6. Study → Stationary
COMSOLの強みはGUI上で弱形式を直接編集できる点で、カスタム方程式の追加が容易だ。
研究用途にはCOMSOLが向いていそうですね。
実験検証を伴う研究では非常に使いやすい。一方、大量の設計ケースを自動で回す場合はAnsys WorkbenchのDesign Explorerの方がワークフローが整っている。
COMSOL Multiphysicsの2D熱モジュール
COMSOLは2005年にFEMLABからブランド名を変更し、2D熱伝導解析のGUIを大幅に改善した。特にパラメトリックスイープ機能により、材料λを1行の設定で100通り変えた感度解析が可能になり、複合材の最適設計に用いる研究者が急増した。現在は世界180か国以上で利用されている。
選定で最も重要な3つの問い
- 「何を解くか」：2次元定常熱伝導に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
等価異方性モデル
2次元の先端的な話題としてはどんなものがありますか？
多層PCBやラミネート複合材料を、等価異方性連続体として2Dモデル化する手法が実務で重要だ。面内と面直の有効熱伝導率は
面内は算術平均（並列）、面直は調和平均（直列）で計算する。
XFEM/GFEMの適用
材料界面やクラックを含む2次元熱伝導問題に拡張有限要素法（XFEM）が適用される。メッシュを界面に合わせなくても界面の温度跳びを表現できる。
界面に合わせなくていいのは楽ですね。
特にクラックや剥離が進展する問題で有効だ。Abaqusでは*XFEMキーワードで利用可能だが、熱解析への適用はまだ研究段階の要素もある。
AIによるメッシュ最適化
最近では、2次元熱伝導問題のメッシュ適応細分化を強化学習で自動化する研究がある。温度勾配の大きい領域を自動検出し、最小の要素数で所望の精度を達成するメッシュを生成する。
AIがメッシュを切ってくれる時代が来るんですね。
Ansys 2025R1ではMesh Morpher AI機能が搭載され、境界層メッシュの自動最適化が実装されている。2D問題は計算コストが低いので、AIメッシュ最適化のテストベッドとしても最適だ。
Coffee Break よもやま話
等温線と流線の直交性
2次元ポテンシャル流と同様に、2次元定常熱伝導の等温線と熱流束線は常に直交する。この性質を利用した「フラックスプロット法」は1940年代にSchlichtingが教育用に普及させ、FEMが普及する前の1970年代まで手書きで複雑形状の熱流束を推定する実用技術として工学部で教えられた。
材料界面やクラックを含む2次元熱伝導問題に拡張有限要素法（XFEM）が適用される。メッシュを界面に合わせなくても界面の温度跳びを表現できる。
界面に合わせなくていいのは楽ですね。
特にクラックや剥離が進展する問題で有効だ。Abaqusでは*XFEMキーワードで利用可能だが、熱解析への適用はまだ研究段階の要素もある。
最近では、2次元熱伝導問題のメッシュ適応細分化を強化学習で自動化する研究がある。温度勾配の大きい領域を自動検出し、最小の要素数で所望の精度を達成するメッシュを生成する。
AIがメッシュを切ってくれる時代が来るんですね。
Ansys 2025R1ではMesh Morpher AI機能が搭載され、境界層メッシュの自動最適化が実装されている。2D問題は計算コストが低いので、AIメッシュ最適化のテストベッドとしても最適だ。
等温線と流線の直交性
2次元ポテンシャル流と同様に、2次元定常熱伝導の等温線と熱流束線は常に直交する。この性質を利用した「フラックスプロット法」は1940年代にSchlichtingが教育用に普及させ、FEMが普及する前の1970年代まで手書きで複雑形状の熱流束を推定する実用技術として工学部で教えられた。
トラブルシューティング
よくある問題
2次元熱伝導でハマりがちなポイントを教えてください。
1. 平面モデルと軸対称モデルの取り違え
問題: Ansys PLANE55で軸対称を指定すべきところを平面のままにすると、円筒形状の熱抵抗が全く異なる値になる。平面問題は単位厚さあたり、軸対称は2πr を考慮する。
対策: KEY OPT設定を必ず確認。Abaqusでは要素名が異なる（DC2D4 vs DCAX4）ので間違いにくい。
2. 異方性材料の座標系
異方性材料で座標系を間違えるとどうなりますか？
材料座標系とグローバル座標系が不一致だと、面内kと面直kが入れ替わる。PCB基板で $k_{in}=30$, $k_{through}=0.5$ W/(m K) が逆になったら温度上昇が数十倍変わる。Ansys WorkbenchではCoordinate Systemの割り当てを要素ごとに確認すること。
3. コーナー特異点
直角コーナーの内角部では理論的に温度勾配（熱流束）が無限大になる。メッシュを細かくするほど値が増大し、収束しない。
対策: 局所的な最大熱流束値は使わず、ある程度の領域で平均化した値を設計に用いる。または微小なフィレットを付けて特異性を除去する。
構造解析の応力特異点と同じ対処法ですね。
まさに同じだ。温度特異点は実物ではフィレットや面取りで緩和されるので、解析モデルにも微小R（0.1mm程度）を入れるのが現実的だ。
角部の特異性とメッシュ依存性
2次元形状の内角が鋭い再入角では理論的に熱流束が無限大になる特異解が生じる。FEMでこの部分を細かくメッシュ分割するほど温度勾配が大きくなり、収束しない現象が起きる。Abaqusマニュアル（2018年版）では再入角近傍への特殊要素（クォータポイント要素）の適用を推奨している。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——2次元定常熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告