2次元定常熱伝導
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世紀の物理学者はこの統一性を「自然の深い美しさ」と称賛した。
2次元定常熱伝導の数値計算手法
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⁹倍以上向上した。
2次元定常熱伝導の実務適用
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次元定常熱伝導のソフトウェア比較
各ツールでの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か国以上で利用されている。
2次元定常熱伝導の先端研究
等価異方性モデル
2次元の先端的な話題としてはどんなものがありますか？
多層PCBやラミネート複合材料を、等価異方性連続体として2Dモデル化する手法が実務で重要だ。面内と面直の有効熱伝導率は
面内は算術平均（並列）、面直は調和平均（直列）で計算する。
XFEM/GFEMの適用
材料界面やクラックを含む2次元熱伝導問題に拡張有限要素法（XFEM）が適用される。メッシュを界面に合わせなくても界面の温度跳びを表現できる。
界面に合わせなくていいのは楽ですね。
特にクラックや剥離が進展する問題で有効だ。Abaqusでは*XFEMキーワードで利用可能だが、熱解析への適用はまだ研究段階の要素もある。
AIによるメッシュ最適化
最近では、2次元熱伝導問題のメッシュ適応細分化を強化学習で自動化する研究がある。温度勾配の大きい領域を自動検出し、最小の要素数で所望の精度を達成するメッシュを生成する。
AIがメッシュを切ってくれる時代が来るんですね。
Ansys 2025R1ではMesh Morpher AI機能が搭載され、境界層メッシュの自動最適化が実装されている。2D問題は計算コストが低いので、AIメッシュ最適化のテストベッドとしても最適だ。
等温線と流線の直交性
2次元ポテンシャル流と同様に、2次元定常熱伝導の等温線と熱流束線は常に直交する。この性質を利用した「フラックスプロット法」は1940年代にSchlichtingが教育用に普及させ、FEMが普及する前の1970年代まで手書きで複雑形状の熱流束を推定する実用技術として工学部で教えられた。
2次元定常熱伝導のトラブル対応
よくある問題
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年版）では再入角近傍への特殊要素（クォータポイント要素）の適用を推奨している。
なった
詳しく
報告