放射冷却 — CAE用語解説
放射冷却
輻射による冷却現象とCAE応用
放射冷却って夜に霜が降りるやつですよね？工業CAEで関係ありますか？
理論と物理
放射冷却の基本概念
「放射冷却」って、夜に地面が冷える現象のことですか？ それとも、もっと広い工学の概念なんですか？
気象現象としての放射冷却はその一例です。工学、特にCAEでは、物体が周囲空間に電磁波（主に赤外線）を放射することで熱を失う「放射伝熱」の一種として捉えます。例えば、宇宙空間に露出した人工衛星のパネルは、太陽からの輻射加熱と、深宇宙（約3K）への放射冷却のバランスで温度が決まります。
輻射と対流や伝導とどう違うんですか？ 数式的にはどう表すんですか？
根本的な違いは、熱の移動に媒体を必要としない点です。真空でも熱は移動します。基本的な式はステファン-ボルツマンの法則で、単位面積あたりの放射エネルギー
放射率εって何で決まるんですか？ 材料によって全然違う？
表面の材質、粗さ、酸化状態、温度、波長に強く依存します。例えば、研磨されたアルミニウムは約0.05と低く、黒体塗装（アクリル系ブラック）は0.95以上です。CAEでは、材料ライブラリにこれらの値が登録されていますが、実測値を使うのが理想です。ASTM E408やJIS R 3106などで測定法が規格化されています。
数値解法と実装
放射伝熱の離散化と計算負荷
放射の計算は、普通の熱伝導のFEMと比べて何が難しいんですか？
最大の難点は「非局所性」です。伝導は隣接要素にしか影響しませんが、放射は離れた面同士が直接「見え合って」熱を交換します。この見え合い関係を表す「形態係数（View Factor）」の計算が、幾何学的に複雑だと非常に重くなります。N個の面がある場合、理論上はN²の組み合わせを計算する必要があります。
N²って…現実的なモデルでは計算できないのでは？
その通りです。そこで「レイトレーシング法」や「ヘミキューブ法」などの近似アルゴリズムが使われます。例えば、Ansys Mechanicalの「Radiation」機能では、形態係数を事前計算するオプションがあり、メッシュが変わらない限り再利用できます。COMSOLでは「放射伝熱」インターフェースで、面と面の可視判定に高速な幾何学アルゴリズムを採用しています。
放射と伝導・対流が組み合わさった問題は、どう解くんですか？ 非線形が強そうです。
連立して解きます。支配方程式は、例えば節点iの熱平衡で書くと、
実践ガイド
放射冷却解析のワークフロー
実際に放射冷却を含む解析をやる時、最初に何を確認すべきですか？
まず、放射が支配的かどうかの見極めです。対流熱伝達率が10 W/m²K以上ある環境では、放射の影響は相対的に小さくなりがちです。しかし、真空や低圧力環境、あるいは高温（400°C以上）では無視できません。LED照明のヒートシンク解析でも、自然対流と放射冷却は同等の寄与を持つことがあります。
モデルを簡素化するコツはありますか？ 全部の面に放射を定義するのは大変そうです。
重要なのは「放射ネットワーク」を考えることです。熱的に重要な面（高温部、放熱フィン表面）と、それらが見ている対象（筐体内部、外部環境）だけを放射面として定義します。背面や隠れた面は無視できます。また、複雑な幾何形状は「包絡面」で置き換えて形態係数を計算すると効率的です。
境界条件で「環境温度」を設定する時、放射用の環境温度と対流用の流体温度は同じでいいんですか？
必ずしも同じではありません。対流は周囲の空気温度、放射は物体が「見ている」面の平均温度です。屋内で壁が室温なら同じですが、屋外では放射のシンク温度は実効天空温度（大気温度より数度低い）を使うことがあります。Ansysの「Radiation to Ambient」では、この区別を設定できます。
ソフトウェア比較
主要CAEソフトにおける放射機能
Ansys、Abaqus、COMSOLで放射の扱い方に違いはありますか？
大きく異なります。Ansys Mechanicalは「Surface to Surface Radiation」と「Radiation to Ambient」を提供。前者は閉空間内の面間放射に強く、形態係数計算に専用ソルバーを使用します。Abaqus/Standardでは「*RADIATION」というキーワードで定義され、主にキャビティ放射に用いられます。COMSOLは「熱放射」が伝熱モジュールに統合され、面間放射と環境への放射を柔軟に組み合わせられ、非等温流体との連成も直感的です。
専用の熱輻射解析ソフトはあるんですか？
あります。ThermoAnalyticsの「RadTherm」や、シーメンスの「STAR-CCM+」に含まれる放射伝熱モジュールは、自動車のサーマルマネジメントなど、複雑な放射・対流連成問題に特化しています。特に波長依存性や散乱を考慮した放射（例えば太陽光のスペクトル）を扱う場合、これらの専用ツールの精度が高いです。
無償・低価格ソフトではどうでしょう？
OpenFOAMの「radiativeHeatTransfer」ソルバーや、Elmer FEMの「Radiation Solver」で基本的な面間放射は計算可能です。ただし、形態係数の計算アルゴリズムやユーザーインターフェースは商用ソフトに比べて劣り、設定が煩雑です。学習用や単純形状の検証には有用ですが、実務での複雑な3Dモデルには不向きな場合が多いです。
トラブルシューティング
放射解析のよくあるエラーと対策
放射を入れたら計算が全然収束しなくなりました。まず疑うべき点は？
まず、放射率εが物理的にあり得ない値（>1や<0）になっていないか確認。次に、形態係数の合計が1から大きく外れていないかチェック。これは「見落とし」がある証拠で、エネルギー保存則を破り、発散の原因になります。ソフトウェアの「View Factor Check」機能を使い、0.95〜1.05の範囲に収まっているか確認しましょう。
「放射による熱流束が大きすぎる」という警告が出ました。なぜですか？
主に二つの原因。第一に、温度の初期値が現実離れしている（例：室温20°Cの部屋で、初期温度を1000°Cに設定）。T^4に比例するので、初期段階で莫大な熱が放射され、数値的に不安定になります。第二に、放射を定義した面同士が極端に近接している（メッシュサイズ以下）場合、形態係数が異常値になることがあります。面間距離を確認してください。
真空チャンバー内の装置の温度を解析しています。放射だけなのに、定常状態で温度が下がり続けます。
それは「閉じた系」になっていない典型的な症状です。装置から放射された熱は、チャンバー内壁に吸収され、内壁の温度が上昇します。その内壁からの再放射（バックヒート）を考慮していないと、熱が無限に宇宙に消えていくようにモデル化され、永遠に冷却されます。内壁を放射面として含め、内壁の熱容量や外部への放熱（冷却水など）も正しくモデル化する必要があります。
放射を考慮した計算時間が想定より10倍かかっています。高速化の方法は？
1. **形態係数の再計算をオフにする**: メッシュや幾何が変わらない定常解析では、一度計算した形態係数をファイルから読み込む。2. **放射面のメッシュを粗くする**: 熱伝導には細かいメッシュが必要でも、放射計算用の面メッシュは粗くできる。3. **対称性を利用する**: 幾何学的に対称なら、放射計算も対称セクターで行う。4. **ソルバー設定**: 放射マトリックス用のソルバーを「反復法」から「直接法」に変える（メモリは食うが収束が早い場合あり）。Ansysでは「Radiation Matrix Method」の選択肢を確認しましょう。
関連トピック
なった
詳しく
報告