形状係数 — 理論と支配方程式

カテゴリ: 伝熱解析 | 2026-01-15
shape-factor-theory
理論と物理の世界へ

形状係数(Shape Factor)とは

🧑‍🎓

先生、形状係数って何に使うんですか?


🎓

2次元や3次元の定常熱伝導問題を、1次元の熱抵抗に帰着させるための概念だ。形状係数 $S$ を使えば


$$q = kS\Delta T$$

と表現でき、$\Delta T$ から直接放熱量が分かる。熱抵抗は $R = 1/(kS)$ だ。


🧑‍🎓

複雑な形状を1つの数値に集約できるんですね。


🎓

$S$ は形状と境界条件だけで決まる幾何学的な量で、単位は [m](2D問題では [m/m] = 無次元/単位奥行)だ。


代表的な形状係数

形状$S$適用条件
無限平板$A/L$基本形
同心円筒$2\pi L / \ln(r_2/r_1)$長さ $L$
同心球$4\pi r_1 r_2/(r_2 - r_1)$
埋設球(半無限体表面から深さ$z$)$4\pi r / (1 - r/(2z))$$z > r$
埋設円筒(半無限体)$2\pi L / \cosh^{-1}(z/r)$$z > r$, $L \gg r$
2つの平行円筒$2\pi L / \cosh^{-1}((d^2-r_1^2-r_2^2)/(2r_1 r_2))$中心間距離 $d$
🧑‍🎓

埋設球や埋設円筒の式が実用的そうですね。


🎓

地中埋設配管の放熱量計算、建物基礎からの地中伝熱、地熱ヒートポンプの設計などで使われる。手計算で概算が得られるのが形状係数の最大の利点だ。


形状係数の導出

🎓

ラプラス方程式 $\nabla^2 T = 0$ の解から導出する。同心円筒の場合


$$T(r) = T_1 + \frac{T_2 - T_1}{\ln(r_2/r_1)} \ln\frac{r}{r_1}$$
$$q = -kA\frac{dT}{dr}\bigg|_{r=r_1} = \frac{2\pi k L(T_1 - T_2)}{\ln(r_2/r_1)}$$

よって $S = 2\pi L / \ln(r_2/r_1)$ が得られる。


🧑‍🎓

解析解から形状係数を逆算するんですね。


🎓

解析解が得られない形状ではFEMで数値的に $S = q/(k\Delta T)$ を算出する。

Coffee Break よもやま話

チャレンジャー号事故とOリングの温度

1986年のスペースシャトル・チャレンジャー号の爆発事故は、低温でOリングのゴムが硬化し、シール機能を失ったことが原因。打ち上げ当日の気温は0°C付近——設計想定を大きく下回っていました。現代の熱-構造連成解析なら「0°Cでゴムの弾性率がどう変わるか」「シール面の接触圧が維持されるか」を事前に検証できます。温度依存材料特性の重要性を、最も痛ましい形で教えてくれた事故です。

各項の物理的意味
  • 蓄熱項 $\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$)で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
  • フーリエの法則:熱流束が温度勾配に比例する線形関係(極低温・超短パルス加熱では非フーリエ熱伝導が必要)
  • 等方性熱伝導:熱伝導率が方向に依存しない(複合材料・単結晶等では異方性を考慮)
  • 温度独立物性値(線形解析):物性値が温度に依存しない仮定(大温度差では温度依存性が必要)
  • 熱放射の扱い:表面間放射はビューファクタ法、参加媒体ではDO法やP1近似を適用
  • 適用外ケース:相変化(融解・凝固)では潜熱の考慮が必要。極端な温度勾配では熱応力連成が必須
次元解析と単位系
変数SI単位注意点・換算メモ
温度 $T$K(ケルビン)またはCelsius絶対温度と摂氏の混同に注意。輻射計算では必ず絶対温度を使用
熱伝導率 $k$W/(m·K)鋼: 約50、アルミ: 約237、空気: 約0.026
熱伝達係数 $h$W/(m²·K)自然対流: 5〜25、強制対流: 25〜250、沸騰: 2,500〜25,000
比熱 $c_p$J/(kg·K)定圧比熱と定積比熱の区別(気体で重要)
熱流束 $q$W/m²境界条件としてのNeumann条件

数値例:平板の定常熱伝導(厚み10mm, 鋼k=50W/(m·K), 表面100°C/裏面20°C)

熱流束 q = k×ΔT/L = 50×80/0.01 = 400,000 W/m² 各位置の温度は線形分布

材料別の熱伝導率の比較(数値が大きいほど熱を伝えやすい):

386 W/(m·K)
アルミニウム237 W/(m·K)
50 W/(m·K)
ガラス1.0 W/(m·K)
空気0.026 W/(m·K)

銅は空気の約15,000倍も熱を伝えやすい! ヒートシンクに銅やアルミが使われる理由がこのグラフで一目瞭然です。

簡易計算ツール:伝熱工学基礎

ビオ数 Bi = hL/k を計算し、集中熱容量法の適用可否を判定します(Bi < 0.1で適用可)。

フーリエ数 Fo = αt/L² を計算(α = k/(ρc_p))。非定常熱伝導の無次元時間。

熱解析の境界条件設定は経験と試行錯誤の繰り返し。 — Project NovaSolverは、実務者の知見を活かしやすい解析環境の実現を研究しています。

CAEの未来を、実務者と共に考える

Project NovaSolverは、形状係数における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。

プロジェクトの最新情報を見る →