臨界断熱半径
臨界断熱半径の理論基礎
臨界断熱半径とは
先生、断熱材を巻けば必ず放熱が減ると思ってたんですが、逆に増えることがあるって本当ですか？
本当だ。円筒や球の外面に断熱材を巻くと、伝導の熱抵抗は増えるが外表面積も増大する。対流の熱抵抗は面積に反比例するから、ある半径までは対流抵抗の減少が伝導抵抗の増加を上回り、全体の放熱量がかえって増える。
直感に反しますね。
放熱量が最大になる半径を臨界断熱半径 $r_{cr}$ という。円筒の場合、全熱抵抗は
これを $r$ で微分してゼロとおくと
球の場合は $r_{cr} = 2k/h$ になる。
$k$ が小さくて $h$ が大きいほど臨界半径が小さいんですね。
そう。自然対流（$h \approx 5$ W/(m$^2$ K)）でエポキシ断熱材（$k \approx 0.2$ W/(m K)）なら $r_{cr} = 0.04$ m = 40 mm だ。元の配管外径が40mm未満なら、断熱材を巻いても逆効果になる領域が存在する。
物理的解釈
全熱抵抗を分解して考えると理解しやすい。
| 半径 $r$ | 伝導抵抗 | 対流抵抗 | 全抵抗 | 放熱量 |
|---|---|---|---|---|
| $r < r_{cr}$ | 増加（小） | 減少（大） | 減少 | 増大 |
| $r = r_{cr}$ | — | — | 最小 | 最大 |
| $r > r_{cr}$ | 増加（大） | 減少（小） | 増加 | 減少 |
対流抵抗の減少が伝導抵抗の増加に勝つかどうかの勝負なんですね。
その通り。実用的には、電線の被覆やパイプ断熱で問題になる。設計時に $r_i > r_{cr}$ であることを確認してから断熱材厚を決定する。
rcrit = λ/h の直感的意味
円筒断熱の臨界半径rcrit = λins/hoは「断熱材の熱伝導抵抗増加速度」と「外表面の対流熱伝達増加速度」が釣り合う点だ。典型的な屋外条件（ho = 10 W/m²·K）とグラスウール（λ=0.04 W/m·K）ではrcrit = 4 mm。多くの産業配管の外半径はこれをはるかに超えるため実問題では出現しないが、電線被覆や医療チューブの設計では重要になる。
臨界断熱半径の数値計算手法
解析的解法
臨界断熱半径は解析解がきれいに出る問題なんですね。
そう。1次元円筒座標の定常熱伝導として厳密解が得られる。温度分布は
放熱量 $q$ は
$r_o$ を変数として $dq/dr_o = 0$ から臨界半径が導かれる。
数値解析による検証
FEMでも同じ結果が出ますか？
当然出る。むしろ理論解との比較はソルバー検証のベンチマークとして使われる。Ansys MechanicalでSOLID70の円筒メッシュを作り、内面に温度拘束、外面に対流条件を設定すれば良い。
| パラメータ | 値 |
|---|---|
| 内径 $r_i$ | 5 mm |
| 断熱材 $k$ | 0.2 W/(m K) |
| 外部 $h$ | 10 W/(m$^2$ K) |
| $r_{cr}$ 理論値 | 20 mm |
断熱材厚を5mmから50mmまで変化させてパラメトリック解析すると、放熱量が $r_o = 20$ mm で最大になることが確認できる。理論値との誤差は0.1%以下だ。
検証用にはうってつけの問題ですね。
APDLマクロでDOループを組めば全ケースを自動で回せる。Abaqusでも同様にPython scriptingで $r_o$ をパラメータにしたスタディが可能だ。
温度依存性を含む場合
断熱材の $k$ が温度依存する場合、$r_{cr} = k(T)/h$ の $k$ をどの温度で評価するかが問題になる。Newton-Raphson反復で自己無撞着解を求める必要がある。
現実の断熱材は温度が上がると $k$ も上がることが多いですよね。
グラスウールは200℃で常温の約2倍になる。高温配管の断熱設計では温度依存性を無視すると過小設計になりかねない。
熱損失曲線のプロットで視覚確認
臨界半径の存在を実感するには、断熱外径を横軸に熱損失Qを縦軸にプロットする。Q(r)=2πLΔT/[ln(r/ri)/λ + 1/(h·r)]を微分してdQ/dr=0から臨界点を求められる。1970年代のInciropera & DeWittの熱伝達教科書（現在第7版まで改訂）は、このグラフを標準例題として世界中の熱工学教育に普及させた。
臨界断熱半径の実務適用
設計への適用
臨界断熱半径を実務でどう活用するんですか？
主に2つの場面がある。
1. 断熱設計: 配管やダクトの断熱材厚を決めるとき、$r_i > r_{cr}$ であることを確認してから厚み最適化
2. 放熱設計: 電線の被覆厚を決めるとき、$r_i < r_{cr}$ なら被覆を厚くしても放熱が改善される
放熱設計に逆利用するのは面白いですね。
電線の被覆が典型例だ。AWG24の銅線（外径0.56mm）にPVC被覆（$k = 0.16$ W/(m K)）を巻く場合、自然対流 $h = 10$ W/(m$^2$ K) として $r_{cr} = 16$ mm。被覆外径が32mmに達するまでは被覆を厚くするほど放熱が良くなる。
代表的な断熱材の物性
| 断熱材 | $k$ [W/(m K)] | 使用温度範囲 | $r_{cr}$（$h$=10） |
|---|---|---|---|
| グラスウール | 0.04 | 〜450℃ | 4 mm |
| ロックウール | 0.04 | 〜700℃ | 4 mm |
| 発泡ウレタン | 0.02 | 〜100℃ | 2 mm |
| シリカエアロゲル | 0.015 | 〜650℃ | 1.5 mm |
| セラミックファイバー | 0.08 | 〜1200℃ | 8 mm |
高性能断熱材ほど $r_{cr}$ が小さいから、細い配管でも安心して使えるんですね。
エアロゲルは $r_{cr} = 1.5$ mm なので、外径3mm以上の配管ならほぼ問題にならない。コストは高いが宇宙・LNG分野で実績がある。
結果の検証
臨界断熱半径の解析結果検証ポイントは以下だ。
- 理論値との比較: $r_{cr} = k/h$（円筒）または $2k/h$（球）と一致するか
- 放熱量カーブ: $r = r_{cr}$ でピーク、その前後で単調変化するか
- エネルギー収支: 内面の入熱と外面の放熱が一致するか
FEMの結果がきれいなカーブにならない場合は何が原因ですか？
メッシュ不足、特に薄い断熱層に要素が1層しかない場合に精度が落ちる。径方向に最低3要素以上を確保すること。
医療用カテーテルの断熱設計
血管内を通る加温カテーテル（外径2〜4 mm）は臨界半径の影響を受ける典型例だ。医療用シリコーン（λ≈0.2 W/m·K）の被覆を厚くしすぎると体液との対流熱伝達面積増大で逆に熱損失が増える。米国FDA提出の510(k)申請資料では断熱設計根拠にFourier熱伝導解析の添付が求められるようになった。
臨界断熱半径のソフトウェア比較
商用ツールでのパラメトリック解析
臨界断熱半径の確認を商用ツールでやるにはどうしますか？
断熱材厚をパラメータにしたスイープ解析が直接的だ。
| ツール | パラメトリック手法 | 所要時間目安 |
|---|---|---|
| Ansys Mechanical | APDLマクロでDOループ | 10ケース×数秒 |
| Ansys Workbench | Design Explorerでパラメータスイープ | 自動メッシュ更新込み |
| COMSOL | Parametric Sweepスタディ | GUIで断熱厚を変数定義 |
| Abaqus | Python scripting + DCAX4要素 | parametric studyで自動化 |
Ansys Workbenchのパラメトリック解析は便利ですね。
WorkbenchではCADパラメータ（SpaceClaimの断熱材厚寸法）をそのまま設計変数にできる。メッシュ更新から求解まで全自動で、結果をResponse Surfaceとして可視化できる。
APDL実装例
最小限のAPDLパラメトリック解析コードはこうだ。
```
*DO,t_ins,1,50,1 ! 断熱材厚1〜50mm
/PREP7
ET,1,PLANE55,,,1 ! 軸対称
MP,KXX,1,0.2 ! 断熱材
CYL4,0,0,5,0,5+t_ins,90
ESIZE,0.5
AMESH,ALL
/SOL
D,NODE(5,0,0),,100 ! 内面100℃
SFL,LINE(外面),CONV,10,25
SOLVE
*GET,Q_total,FSUM,,HEAT
*ENDDO
```
1mmきざみで50ケース回すんですね。
この程度の問題なら全体で数分で完了する。結果をプロットすれば $r_{cr}$ の位置が一目瞭然だ。
ツール選定の指針
| シナリオ | 推奨 |
|---|---|
| 単純円筒の概算 | 手計算（$r_{cr}=k/h$）で十分 |
| 温度依存性・放射を考慮 | Ansys Mechanical, COMSOL |
| 3D配管の断熱最適化 | Workbench + Design Explorer |
| 配管系全体の熱損失評価 | Flownex, COMSOL Pipe Flow |
手計算で済む問題はわざわざFEMを回さなくていいと。
そう。工学的判断力とは問題に対して適切な道具を選ぶことだ。
Wolfram Mathematicaで解析解を即時可視化
Wolfram Mathematica 12（2019年）のManipulate関数を使うと、rcrit = λ/hのパラメータをスライダーで動かしながら熱損失曲線をリアルタイムで描画できる。大学の熱工学講義でこのノートブックを配布する教員が増えており、MITのOCWでも2021年から関連Mathematicaノートを公開している。
臨界断熱半径の先端研究
放射を含む臨界断熱半径
高温環境だと放射の影響も出てきますよね？
そうだ。外表面からの放射を含めると全熱抵抗の解析が複雑になる。放射の線形化熱伝達係数は
全外面熱伝達係数 $h_{\text{eff}} = h_{\text{conv}} + h_r$ として臨界半径を求めると
$h_r$ は温度依存なので反復計算が必要だが、対流だけの場合より $r_{cr}$ が小さくなる。
放射を含めると臨界半径が小さくなるんですね。
ステファン・ボルツマンの法則 $q = \varepsilon \sigma A T^4$ で放射は面積に比例する。表面積増大による放熱増加がさらに大きくなるからだ。300℃以上の高温配管では放射の寄与が対流と同程度になり、無視できない。
多層断熱の最適化
複数の断熱材を重ねる場合、材料の配置順序で性能が変わる。
高温側に耐熱性の高い材料（セラミックファイバー）、低温側に低熱伝導率材料（エアロゲル）を配置するのが一般的だ。
層構成の最適化は手計算だと大変そうですね。
ExcelのソルバーやPythonのSciPy.optimizeで十分対応できる。3層以上ならCOMSOLのOptimization Moduleで自動最適化するのも効率的だ。
真空断熱への発展
真空断熱パネル（VIP）は $k = 0.004$ W/(m K) と極端に低く $r_{cr} = 0.4$ mm。ほぼ全ての形状で臨界断熱厚は無問題だ。LNG輸送タンクや宇宙機の多層断熱（MLI）では有効熱伝導率が $10^{-4}$ W/(m K) オーダーになる。
宇宙機の断熱設計は次元が違いますね。
真空環境では対流がないため、放射遮蔽が主な断熱機構になる。地上の臨界断熱半径の概念は適用外だが、放射シールド枚数の最適化という類似の問題が生じる。
球殻断熱には臨界半径が存在しない
円筒断熱とは異なり、球殻断熱では外半径を増やしても熱損失は単調減少する（極値なし）。これは球の表面積が半径の2乗で増えるのに対し、断熱抵抗も2乗成分で増えるためバランスが崩れないからだ。Cengel & Ghajarの熱伝達教科書（2020年第6版）は球殻と円筒を比較した図解でこの直感を整理している。
臨界断熱半径のトラブル対応
よくあるトラブルと対策
臨界断熱半径まわりで実務上困ることってありますか？
典型的な問題を整理しよう。
1. 断熱材を追加したのに温度が下がった
原因: 元の管径が $r_{cr}$ 未満で、断熱材の追加が放熱を増大させた。
対策: $r_{cr} = k/h$ を計算して現在の管径と比較。$r_i < r_{cr}$ なら、より低い $k$ の断熱材を選定するか、外面に遮風カバーを付けて $h$ を下げる。
2. 断熱材が結露する
低温配管の断熱で結露が問題になることがありますよね。
断熱材内部の温度が露点を下回ると内部結露が生じる。繊維系断熱材は含水で $k$ が急増し、さらに結露が進行する悪循環に陥る。
対策: 防湿層（PE/ALフィルム）を断熱材の外面（温度が高い側）に設ける。低温配管では必須だ。
3. 解析結果が理論値と合わない
FEM結果との不一致チェックリスト。
| チェック項目 | よくある原因 |
|---|---|
| 軸対称設定 | 2D要素でaxisymmetric設定を忘れている |
| 対流面 | 内面に対流条件を間違えて設定 |
| 単位系 | mm系で $k$ の値を変換し忘れ |
| メッシュ | 薄い断熱層に要素が1層しかない |
| 端面条件 | 端面が断熱（デフォルト）になっているか |
軸対称の設定忘れはありがちですね。
2D平面のデフォルトで計算すると面積変化の効果が入らないので、臨界断熱半径の現象が再現できない。Ansys PLANE55ならKEYOPT(3)=1、Abaqus DCAX4なら軸対称要素を明示的に選択すること。
細い電線の過熱トラブルに注意
直径3 mm以下の電線に塩化ビニール絶縁被覆（λ≈0.17 W/m·K）を巻く場合、被覆厚さを増やすと発熱体（電線）の温度が上がるケースがある。UL規格では絶縁被覆の最適厚さを電流容量ごとに規定しているが、カスタム設計品で規格外の組み合わせを選んだ際に過熱トラブルが生じた事例がUL Whitepaperに記録されている。
なった
詳しく
報告