フィン伝熱解析
フィン伝熱の理論基礎
フィンの伝熱メカニズム
先生、フィンの中では伝導と対流が同時に起きてるんですよね？
そう。フィン根元からベース温度の熱が伝導で先端方向に運ばれ、フィン表面から対流で周囲流体に放出される。この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を示す。
フィン伝熱の数値計算手法
境界条件と解析解
フィンの境界条件にはどんな種類がありますか？
根元は $\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を達成する条件を計算で求められる。
フィン伝熱の実務適用
フィン設計の最適化
フィンの長さや厚みの最適値はどう決めるんですか？
$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は圧力損失と熱抵抗のトレードオフで決定される実務的な値だ。
フィン伝熱のソフトウェア比較
商用ツールでのフィン解析
フィン解析に適したツールはどれですか？
解析の目的で使い分ける。
| 目的 | 推奨ツール | 理由 |
|---|---|---|
| 伝導のみ（固定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は無償でフィン効率・有効度の一括パラメトリック計算が行える。
フィン伝熱の先端研究
非一様断面フィン
テーパフィンや三角形フィンの方が効率がいいんですよね？
理論的にはそうだ。断面積が先端に向かって減少するフィンは、材料を根元に集中配置するため同じ体積で高い効率が得られる。
| 形状 | 相対材料量 | $\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）機能でこの問題を自動的に解消できる。
関連トピック
なった
詳しく
報告