フィン伝熱解析
理論と物理
フィンの伝熱メカニズム
先生、フィンの中では伝導と対流が同時に起きてるんですよね？
そう。フィン根元からベース温度の熱が伝導で先端方向に運ばれ、フィン表面から対流で周囲流体に放出される。この2つのバランスがフィン性能を決める。支配方程式はエネルギー保存から導かれる。
$A_c$ がフィン断面積、$P$ がフィン周長（濡れ縁長さ）だ。一様断面で $k$ が一定なら
$m$ の物理的な意味は何ですか？
$1/m$ がフィンの特性長さで、温度が $1/e$ に減衰する距離の目安だ。$m$ が大きいほど温度が急減衰する。つまり細くて薄い（$A_c$ 小）、表面が広い（$P$ 大）フィンほど $m$ が大きい。
各種断面の $m$ 値
| フィン断面 | $A_c$ | $P$ | $m$ |
|---|---|---|---|
| 矩形（幅$w$, 厚$t$） | $wt$ | $2(w+t)$ | $\sqrt{2h(w+t)/(kwt)}$ |
| 薄板（$w \gg t$） | $wt$ | $\approx 2w$ | $\sqrt{2h/(kt)}$ |
| 円柱（直径$d$） | $\pi d^2/4$ | $\pi d$ | $\sqrt{4h/(kd)}$ |
薄板フィンだと $m = \sqrt{2h/(kt)}$ で厚み $t$ にだけ依存するんですね。
フィン厚が支配パラメータだ。$t$ を半分にすると $m$ が $\sqrt{2}$ 倍、フィン効率は低下する。効率と材料使用量のトレードオフが設計の核心だ。
フィン効率の定義と意味
フィン効率ηは「実際の放熱量÷フィン全体が基部温度であった場合の放熱量」で定義される。1938年にHarper & Brownが定式化したこの無次元指標は、kが高いほど1に近づき銅フィンで一般に0.95〜0.99を示す。
各項の物理的意味
- 蓄熱項 $\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$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
境界条件と解析解
フィンの境界条件にはどんな種類がありますか？
根元は $\theta(0) = \theta_b = T_b - T_\infty$ で固定。先端の条件は4種類ある。
Case A: 断熱先端
Case B: 一定温度先端
Case C: 対流先端
修正長さ $L_c = L + t/2$（矩形）で断熱先端の解を適用。
Case D: 無限長フィン
実務ではCase Aの修正版（Case C）を使うことが多いですか？
そう。先端面積がフィン側面に比べて小さければ修正長さの近似で十分だ。先端を別途モデル化する必要はまずない。
FEMでの解法
FEMでフィンを解く場合、表面の対流条件は
を熱負荷ベクトルに追加する。対流境界の要素（Surface Effect要素：AnsysのSURF152など）が自動的にこの項を生成する。
解析解があるのにFEMで解く意味はありますか？
テーパフィン、放射を含むフィン、温度依存の $h$ や $k$ のような非線形問題は解析解では扱えない。FEMならこれらを自然に取り込める。
フィン有効度の計算方法
フィン有効度εはフィン付き面の熱流束÷フィンなし面の熱流束。εが2未満のフィンは費用対効果が乏しく実用に適さない。熱伝達係数h=50 W/m²Kのとき、厚さ1mmのアルミフィンでε≈15を達成する条件を計算で求められる。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
フィン設計の最適化
フィンの長さや厚みの最適値はどう決めるんですか？
$mL$ が設計の指標だ。
| $mL$ | フィン効率 $\eta_f$ | 評価 |
|---|---|---|
| 0.5 | 0.92 | 効率良好だがやや短い |
| 1.0 | 0.76 | コスト対効果の最適点 |
| 1.5 | 0.57 | 先端が冷えている |
| 2.0 | 0.48 | 先端1/3がほぼ無駄 |
| 3.0 | 0.33 | 長すぎ |
$mL = 1$ が最適なんですね。
「最適」の定義によるが、材料あたりの放熱量（$q_f / V_{\text{fin}}$）が最大になるのは $mL \approx 1.4$ だ。フィン数を増やす余裕があれば短いフィンを多数配置する方が効率的だ。
代表的な応用例
| 応用 | フィン材質 | 典型的な $mL$ | 備考 |
|---|---|---|---|
| CPU ヒートシンク | Al/Cu | 0.8〜1.2 | 強制対流、スカイブフィン |
| 空調フィンコイル | Al | 0.5〜1.0 | 薄板フィン+銅管 |
| 発電機冷却フィン | 鋳鉄 | 1.0〜2.0 | 放射+自然対流 |
| 宇宙用ラジエータ | Al/CFRP | 0.5〜1.5 | 放射のみ |
空調のフィンコイルは街中でよく見ますね。
エアコンの室外機のアルミフィンがまさにそれだ。フィンピッチ1.5〜2mm、板厚0.1〜0.15mmという極薄設計で、数百枚のフィンが銅管に圧入されている。
検証のポイント
フィン解析結果の検証は以下を確認する。
- 根元温度がベース温度と一致するか
- 先端温度が $T_\infty$ 以上か（$T_\infty$ 未満は非物理的）
- 放熱量が $M \tanh mL$ と概ね一致するか
- エネルギー収支（根元からの入熱 = 表面からの対流放熱）
解析解で検算できるのは安心ですね。
フィン問題は数少ない「解析解がある実用問題」なので、FEMの学習・検証に最適だ。
データセンターの冷却フィン設計
AWS東京リージョンのサーバーラックでは1Uあたり200W以上の発熱に対し、厚さ0.4mmのアルミフィンを40枚並べたヒートシンクで対応する。フィン間ピッチ2.5mmは圧力損失と熱抵抗のトレードオフで決定される実務的な値だ。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
商用ツールでのフィン解析
フィン解析に適したツールはどれですか？
解析の目的で使い分ける。
| 目的 | 推奨ツール | 理由 |
|---|---|---|
| 伝導のみ（固定h） | Ansys Mechanical, Abaqus | 解析解との検証に最適 |
| CHT（流体+固体） | Ansys Fluent, STAR-CCM+ | 局所hを自動計算 |
| 電子機器ヒートシンク | Icepak, FloTHERM | 部品ライブラリが充実 |
| パラメトリック最適化 | COMSOL, optiSLang | 形状パラメータの自動最適化 |
APDL実装例
矩形フィン1枚の解析コードはこうだ。
```
/PREP7
ET,1,SOLID70
MP,KXX,1,200 ! Al k=200
BLOCK,0,0.05,0,0.001,0,0.01 ! L=50mm, t=1mm, w=10mm
ESIZE,0.001
VMESH,ALL
/SOL
D,NODE(0,,,),,80 ! 根元80℃
SF,ALL,CONV,25,25 ! h=25, Tinf=25℃
SOLVE
```
SF,ALLで全面に対流条件を掛けるんですね。
根元面は温度拘束されているので、対流条件との競合はない（Dirichletが優先される）。実務ではDA(根元面)を対流条件から除外する方が厳密だが、結果への影響は微小だ。
FloTHERMでのヒートシンクモデル
FloTHERMではHeat Sinkオブジェクトとしてパラメトリック入力が可能だ。
- Base寸法、厚さ
- フィン数、フィン高さ、フィン厚
- フィンタイプ（平板、ピン、楕円）
メッシュはカルテシアン格子で自動生成される。コンパクトモデルに切り替えれば、板レベルの全体解析にも組み込める。
カルテシアン格子だとフィンの斜面は階段状になりませんか？
FloTHERMはSmart Cell技術で境界を補正する。FloTHERM XTではテトラメッシュも使えるので複雑形状にも対応できる。
フィン解析ソフトの選択指針
Mentor Flotherm XTはPCB搭載フィンのミルワット単位の解析に特化し、2024年版でIPC-2581フォーマットの直接インポートに対応した。一方、スタンドアロンのFinnpilotは無償でフィン効率・有効度の一括パラメトリック計算が行える。
選定で最も重要な3つの問い
- 「何を解くか」：フィン伝熱解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
非一様断面フィン
テーパフィンや三角形フィンの方が効率がいいんですよね？
理論的にはそうだ。断面積が先端に向かって減少するフィンは、材料を根元に集中配置するため同じ体積で高い効率が得られる。
| 形状 | 相対材料量 | $\eta_f$（$mL=1$時） | 解法 |
|---|---|---|---|
| 矩形 | 1.00 | 0.762 | 双曲関数 |
| 三角形 | 0.50 | 0.855 | ベッセル関数 |
| 放物線 | 0.33 | 0.903 | ベッセル関数 |
| 双曲線 | 可変 | 場合による | 数値解 |
材料半分で効率が上がるって凄いですね。
問題は製造コストだ。三角形フィンは押出しでは作りにくく、鍛造やアディティブマニュファクチャリングが必要になる。
放射を含むフィン
高温環境や宇宙用途では放射の寄与が大きい。放射を含むフィンの支配方程式は
$T^4$ の非線形項があるため、解析解は得られない。Newton-Raphson反復でFEMを解く。
放射の寄与はどのくらいの温度から重要になりますか？
目安として放射の線形化 $h_r = 4\varepsilon\sigma T_m^3$ で評価する。$T_m = 100$℃なら $h_r \approx 7$ W/(m$^2$ K)で自然対流と同程度。$T_m = 300$℃なら $h_r \approx 20$ で無視できない。
濡れフィン
空調のフィンコイルでは凝縮水がフィン表面を覆う「濡れフィン」状態が生じる。水膜の蒸発による潜熱移動が加わり、顕熱移動の2〜3倍の効果がある。
除湿効果も含めて評価しないといけないんですね。
ASHRAEのHandbook of FundamentalsにWet Fin Efficiencyの相関式が掲載されている。$m_{\text{wet}} = m \sqrt{1 + B'\omega}$ として修正 $m$ を使う。$\omega$ は湿度比、$B'$ は係数だ。
断面最適化フィンの理論
Duffin（1959年）は変分法を用いて「熱流束一定の条件下で最小体積となるフィン断面は三角形」であることを証明した。現代のトポロジー最適化ではこの解析解を初期推定値として用い、ANSYS Mechanicalで検証するアプローチが一般的だ。
トラブルシューティング
よくあるトラブルと対策
フィン解析のトラブルで多いのは何ですか？
典型的な問題を整理する。
1. FEMの放熱量が解析解と合わない
チェックリスト:
- 根元面に対流条件が重複して掛かっていないか
- 先端の境界条件（断熱か対流か）が解析解と一致しているか
- $h$ の単位（W/(m$^2$ K)）が正しいか
- メッシュがフィン厚方向に2要素以上あるか
2. フィン先端の温度が $T_\infty$ 以下になる
原因: 非物理的。対流条件の符号ミスか、$h$ の値が桁違いに大きい可能性。
対策: $h$ を確認。フィン先端温度は物理的に $T_\infty \leq T_{\text{tip}} \leq T_b$ の範囲内でなければならない。
3. 薄いフィンのメッシュ品質
板厚0.5mmのフィンだとメッシュが難しいですよね。
薄肉構造のメッシング戦略はこうだ。
| 方法 | 適用条件 | メリット |
|---|---|---|
| Sweep mesh | 一定断面フィン | 高品質六面体 |
| Inflation layer | フィン表面近傍 | 温度勾配を正確に捉える |
| Shell要素 | $t/L < 0.1$ | 要素数を大幅削減 |
Ansys Mechanicalでは薄肉部にSweep Meshを適用し、厚み方向を2〜3分割するのが標準だ。
4. ヒートシンク全体のCFDが収束しない
原因: フィン間の流路が狭く、メッシュ品質が悪い。特にフィン先端付近で高アスペクト比の要素が発生。
対策: フィン間にInflation Layerを適用。最小メッシュサイズをフィン間隔の1/10以下にする。IcepakではZero-Slackオプションでフィン間の最小メッシュ数を制御できる。
細かいフィンピッチだとCFDが重くなりますね。
対称条件で1ピッチだけモデル化するか、FloTHERMのコンパクトモデルを使って回避する。
フィン解析のメッシュ依存性
フィン先端付近は温度勾配が急峻なため、要素サイズをフィン厚さの1/10以下に細分化しないと熱流束誤差が10%を超えることがある。ANSYS 2024では適応メッシュ精細化（AMR）機能でこの問題を自動的に解消できる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——フィン伝熱解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告