拡大伝熱面（フィン）
理論と物理
拡大伝熱面の基礎
先生、ヒートシンクのフィンってなぜ効くんですか？
対流熱伝達の基本式 $q = hA(T_s - T_\infty)$ で、$h$ を大きくするのが難しいとき面積 $A$ を増やすのが拡大伝熱面（Extended Surface）の戦略だ。ベース面だけでは不十分な放熱を、フィンで表面積を10〜100倍に拡大して達成する。
単純に面積を増やせばいいんですね。
ただしフィンの先端に行くほど温度が下がるため、全面がベース温度で放熱するわけではない。フィン効率 $\eta_f$ で実効性を評価する。
フィンの支配方程式
一様断面の直線フィンの温度分布は、エネルギー保存から
ここで $\theta = T(x) - T_\infty$、$m = \sqrt{hP/(kA_c)}$。$P$ はフィン周長、$A_c$ はフィン断面積だ。
一般解は $\theta(x) = C_1 e^{mx} + C_2 e^{-mx}$ で、境界条件によって定数が決まる。
$m$ はフィンパラメータですね。$m$ が大きいと温度が急激に下がる。
そう。$m$ は「フィンの細さ」を表す。薄くて長い（$A_c$ 小、$P$ 大）フィンほど $m$ が大きく、先端温度が低い。
境界条件と解
| 先端条件 | 温度分布 | 放熱量 |
|---|---|---|
| 断熱先端 | $\theta = \theta_b \frac{\cosh m(L-x)}{\cosh mL}$ | $q = \sqrt{hPkA_c}\,\theta_b \tanh mL$ |
| 一定温度先端 | 双曲関数の線形結合 | 場合に応じた式 |
| 対流先端 | 修正長さ $L_c = L + A_c/P$ で近似 | 断熱先端の式を $L_c$ で適用 |
| 無限長フィン | $\theta = \theta_b e^{-mx}$ | $q = \sqrt{hPkA_c}\,\theta_b$ |
実務では断熱先端の近似が多いですか？
先端からの放熱は全体の数%程度なので、修正長さ $L_c$ を使った断熱先端近似で十分な精度が得られる。
フィン誕生の歴史
拡大伝熱面（フィン）の概念は1922年にアルフレッド・ハーパーが体系化した。今やIntelのCPUクーラーにも採用され、面積を最大10倍に拡張することで放熱量を劇的に増やす技術の根幹をなしている。
各項の物理的意味
- 蓄熱項 $\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$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
フィン効率
フィン効率って具体的に何を表すんですか？
フィン効率 $\eta_f$ はフィン全面がベース温度だった場合の最大放熱量に対する実際の放熱量の比だ。
$$\eta_f = \frac{q_{\text{actual}}}{q_{\text{max}}} = \frac{q_{\text{actual}}}{hA_f \theta_b}$$
フィン効率って具体的に何を表すんですか？
フィン効率 $\eta_f$ はフィン全面がベース温度だった場合の最大放熱量に対する実際の放熱量の比だ。
断熱先端の直線フィンなら $\eta_f = \tanh(mL)/(mL)$。$mL < 1$ でフィン効率90%以上、$mL > 3$ で急激に低下する。
$mL$ がフィン設計のキーパラメータなんですね。
$mL \approx 1$ がコスト対効果の最適点とされている。これ以上長くしても材料の割に放熱が増えない。
FEMでのフィン解析
フィンのFEM解析は薄い構造なのでシェル要素やビーム要素でモデル化できる場合もあるが、熱解析では通常ソリッド要素を使う。
| アプローチ | メリット | デメリット |
|---|---|---|
| 3Dソリッド | 精度最高 | メッシュ数が多い |
| 2D断面 | 繰返し構造を効率的に解析 | 3D効果を無視 |
| 1D解析解 | 高速、パラメスタ容易 | 2D/3D熱流を無視 |
ヒートシンクのフィンが100枚あったら全部モデル化するんですか？
対称性を使えば1枚分のモデルで済む。フィン間のピッチで対称面を設定し、対称面に断熱条件を課す。FloTHERMやIcepakではフィンアレイをパラメトリックモデルとして自動生成する機能がある。
全面フィン効率
ヒートシンク全体の性能はフィン効率とベース面を組み合わせた全面効率 $\eta_o$ で評価する。
$A_f$ はフィン面の合計面積、$A_t$ は全表面積（フィン+ベース露出部）だ。
ベース露出部はフィン効率100%として扱うんですね。
その通り。全面効率を使えば全体の放熱量を $q = \eta_o h A_t \theta_b$ と簡潔に表現できる。
フィン効率の計算手順
フィン効率ηは双曲線関数tanh(mL)/(mL)で表される。mはフィンパラメータで√(hP/kA)。アルミ製ピンフィン（k=237 W/m·K）では効率95%超が得られ、銅（k=401 W/m·K）でも同様に高効率を実現できる。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
ヒートシンクの設計指針
ヒートシンクのフィンを設計するとき、何を基準にすればいいですか？
設計パラメータと影響を整理する。
| パラメータ | 増やすと | トレードオフ |
|---|---|---|
| フィン高さ | 面積増加で放熱向上 | $mL$ 増大でフィン効率低下 |
| フィン数 | 面積増加で放熱向上 | 流路が狭くなり圧力損失増大 |
| フィン厚 | $A_c$ 増加でフィン効率向上 | 重量増加、流路面積減少 |
| フィンピッチ | — | 最適値あり（自然対流で6〜12mm） |
フィンピッチに最適値があるんですね。
自然対流では狭すぎると空気が流れず、広すぎると面積が足りない。Bar-CohenとRohsenowの最適ピッチの相関式は
$L$ はフィン高さ、$\text{Ra}_L$ はレイリー数だ。
材料選定
| 材料 | $k$ [W/(m K)] | 密度 [kg/m$^3$] | 用途 |
|---|---|---|---|
| 銅 C1100 | 398 | 8,960 | 高性能ヒートシンク |
| アルミ A6063 | 200 | 2,700 | 汎用ヒートシンク |
| アルミ A1050 | 230 | 2,710 | ダイキャストフィン |
| グラファイト | 150〜400（面内） | 2,200 | 薄型放熱シート |
銅とアルミで $k$ は2倍ですが、密度は3倍以上違いますね。
重量あたりの放熱性能ではアルミが優れる。航空・車載ではアルミが主流で、データセンターのサーバ冷却のように重量制約が緩い場合に銅を使う。
製造方法と形状制約
製造方法によってフィン形状が制約される。
| 製法 | フィン厚 | アスペクト比 | コスト |
|---|---|---|---|
| 押出し | 1〜3mm | 〜8:1 | 低 |
| ダイキャスト | 1.5〜4mm | 〜6:1 | 中 |
| スカイブ加工 | 0.2〜0.5mm | 〜50:1 | 高 |
| 接合（ろう付） | 0.1〜0.3mm | 自由 | 高 |
スカイブ加工だとアスペクト比50:1のフィンが作れるんですね。
薄いフィンを高密度に配置できるので表面積が大きく取れる。サーバ用CPUクーラーで多用される。
CPUクーラーの設計実務
Ryzen 9 7950X（TDP 170W）向けヒートシンクでは、アルミニウム押出成形フィンを48枚配置し熱抵抗0.1 K/W以下を達成する。フィン間隔2mmで気流速度2 m/sの条件が実務的な最適点として多用される。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
商用ツールでのフィン解析
フィンの熱解析にはどのツールが適していますか？
用途に応じて使い分ける。
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| フィン単体の伝導解析 | Ansys Mechanical, COMSOL | 解析解との比較検証 |
| ヒートシンク+空気のCHT | Ansys Icepak, FloTHERM | 電子冷却に特化 |
| 高精度CFD | Ansys Fluent, STAR-CCM+ | 乱流モデルが充実 |
| パラメトリック最適化 | FloTHERM XT, Icepak | 形状パラメータの自動最適化 |
IcepakとFloTHERMではフィンの扱いが特殊なんですか？
Icepakではヒートシンクをパラメトリックオブジェクトとして定義できる。フィン数、フィン高さ、ピッチ、厚みを入力するだけで3Dモデルとメッシュが自動生成される。Design of Experiments（DOE）と組み合わせて最適形状を探索できる。
FloTHERMでのコンパクトモデル
FloTHERMではフィンアレイをコンパクトモデル（抵抗ネットワーク）で近似する機能がある。フィン1枚ずつメッシュを切らなくても全体の放熱特性を高速に評価できる。
板レベル設計でヒートシンクを100個配置する場合に便利ですね。
その通り。詳細モデルは最終的な形状確認に使い、初期設計段階ではコンパクトモデルで高速に回す。このアプローチがフロントローディング設計の基本だ。
COMSOLでのフィン最適化
COMSOLのOptimization Moduleを使えば、フィン形状（高さ、厚み、テーパ角度）を目的関数（熱抵抗最小化や重量最小化）のもとで自動最適化できる。Topology Optimizationでフィン形状を自由に探索することも可能だ。
直線フィンが最適とは限らないんですね。
トポロジー最適化すると樹枝状（デンドライト）の構造が出現することがある。アディティブマニュファクチャリングなら製造可能だが、従来工法では制約が大きい。
フィン解析対応ソフト比較
ANSYS Icepakはフィン形状のパラメトリック最適化に強く、2023年版ではAI補助メッシュ生成を搭載。Siemens Simcenter Flothermは基板レベルの実装解析に定評があり、どちらもJEDEC標準との連携機能を持つ。
選定で最も重要な3つの問い
- 「何を解くか」：拡大伝熱面に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
放射フィン
宇宙空間ではフィンからの放熱は放射だけですよね。
そう。宇宙用ラジエータは対流が使えないため、放射フィンが必須だ。支配方程式が非線形になる。
$T^4$ の非線形項があるため解析解は一般には得られず、数値解析が必要だ。
$T^4$ は厄介ですね。
線形化（$T^4 \approx 4T_m^3 T$ の近似）すれば対流フィンと同じ形になるが、温度差が大きいと精度が落ちる。宇宙機設計ではSINDA/FLUINTやThermal Desktopが標準ツールだ。
マイクロフィンとマイクロチャネル
マイクロスケール（フィン高さ0.1〜1mm）のフィンは電子デバイスの冷却で活用されている。水冷マイクロチャネルヒートシンクはフィン間の流路幅が0.05〜0.5mmで、$h$ が $10^4$〜$10^5$ W/(m$^2$ K) に達する。
自然対流の1000倍以上ですね。
Tuckerman & Peasの先駆的研究（1981年、MIT）以来、CPUクーラーやパワーモジュールの冷却に応用されている。設計にはAnsys FluentやCOMSOLのMicrofluidics Moduleが使われる。
フィンの構造最適化
最近のトレンドはアディティブマニュファクチャリング（金属3Dプリント）で実現する自由形状フィンだ。
| 構造 | 特徴 | $h$ 向上率 |
|---|---|---|
| ストレートフィン | 基本形。押出しで製造 | 基準 |
| ピンフィン | 円柱状。渦発生で混合促進 | +20〜40% |
| TPMS構造 | 三重周期極小曲面。高面積密度 | +50〜100% |
| ラティス構造 | 格子状。軽量かつ高伝熱 | +30〜80% |
TPMS構造って何ですか？
GyroidやSchwarz Pなどの数学的曲面で、体積あたりの表面積が非常に大きい。nTopologyやCOMSOLで設計し、金属3Dプリントで製造する。EOS M290やSLM Solutions等のL-PBF装置で実用化されている。
最適フィン形状の理論
シュミット（1926年）が証明したように、体積一定の条件で最大放熱量を得る最適フィン形状は放物線断面（パラボラフィン）である。現代のANSYS Fluentでは形状最適化ソルバーを用いてこの理論を数値的に検証できる。
トラブルシューティング
よくあるトラブルと対策
フィンの解析で注意すべき点は何ですか？
頻出の問題を整理しよう。
1. フィン効率の過大評価
原因: $h$ を一様と仮定しているが、実際はフィン根元と先端で異なる。フィン間の流れが十分に発達せず、$h$ が局所的に変化する。
対策: CHT解析（CFD+伝導の連成）で局所 $h$ を自動計算する。解析解は概算としてのみ使う。
2. ヒートシンクの放熱量が計算より低い
解析では100Wの放熱能力があるはずなのに、実測だと70Wしか出ないケースですね。
よくある原因はこうだ。
| 原因 | 影響度 |
|---|---|
| TIMの接触熱抵抗が過小評価 | 高 |
| 実際の風速が設計値より低い | 高 |
| フィン間のバイパス流 | 中 |
| 放射の寄与を無視 | 低〜中 |
| ベースプレートのスプレッディング抵抗 | 中 |
バイパス流って何ですか？
フィンの横から空気が漏れて、フィン間を通らずに迂回する流れだ。ヒートシンクをダクトで囲まないと風がバイパスする。Icepakではバイパス流を含めたシステムレベルの解析ができる。
3. メッシュ関連の問題
薄いフィン（0.5mm以下）のメッシングは難しい。
- フィン厚方向に最低2要素。1要素だと温度勾配を正しく捉えられない
- フィン根元の応力集中部はメッシュを細かく
- 接触面（ベース-フィン接合部）のメッシュ整合を確認
フィンが100枚あるとメッシュ数が膨大になりませんか？
対称条件で1ピッチ分だけモデル化するのが基本だ。全体モデルが必要な場合はFloTHERMのコンパクトモデルが有効だ。
フィン解析の誤差源
フィン基部と本体の接触熱抵抗を無視すると熱抵抗を最大30%過小評価する。Thermal Interface Material（TIM）の熱伝導率は0.5〜8 W/m·Kと幅広く、Bergquist GP3000（k≈3 W/m·K）などの実測値を必ず入力すること。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——拡大伝熱面の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告