灰色体間の放射交換

灰色体放射交換の理論基礎

前提となる3つの仮定

面間の放射交換計算の標準モデル「灰色拡散面のエンクロージャ理論」は、①灰色（放射特性が波長によらない：適用限界はキルヒホッフの法則参照）、②拡散（放射・反射が方向によらない）、③各面が等温・一様放射、の3仮定の上に建っています。この3つを飲めば、任意の面数の放射交換が線形代数に落ちる——工学的にきわめて費用対効果の高い理論です。逆に言えば、結果が実測と合わないときに疑うべき仮定リストでもあります。

形態係数という幾何の言葉

面 \( i \) から出た放射のうち面 \( j \) に直接届く割合が形態係数 \( F_{ij} \) です。純幾何量であり、2つの恒等式が常に成り立ちます。

$$ A_i F_{ij} = A_j F_{ji} \;\;(\text{相反則}), \qquad \sum_j F_{ij} = 1 \;\;(\text{総和則・閉空間}) $$

この2式は計算した形態係数の検算式としてそのまま使えます。総和則が満たされない形態係数セットは、エネルギーを漏らす（または創る）計算につながるため、放射解析の検証はまずここから始めます。

放射抵抗ネットワーク——2面系の完全解

灰色面の交換は電気回路のアナロジーで見通せます。各面に表面抵抗 \( (1-\varepsilon)/(\varepsilon A) \)、面間に空間抵抗 \( 1/(A_1 F_{12}) \) を置くと、2面エンクロージャの正味熱交換は

$$ Q_{12} = \frac{\sigma (T_1^4 - T_2^4)}{\dfrac{1-\varepsilon_1}{\varepsilon_1 A_1} + \dfrac{1}{A_1 F_{12}} + \dfrac{1-\varepsilon_2}{\varepsilon_2 A_2}} $$

です。平行大平板（\( F_{12}=1, A_1=A_2 \)）なら \( q = \sigma(T_1^4 - T_2^4)/(1/\varepsilon_1 + 1/\varepsilon_2 - 1) \)、大空間中の小物体なら包絡側の抵抗が消えて \( Q = \varepsilon_1 A_1 \sigma (T_1^4 - T_2^4) \)——実務頻出の公式群は全てこのネットワークの特殊形です。「表面抵抗はεが小さいほど大きい」という電気回路的直感が、低放射率面（研磨金属・断熱シールド）の効き方をそのまま説明します。

多面系の計算手法

ラジオシティ法——N面系の標準解法

面数が増えたら、各面のラジオシティ \( J_i \)（放射+反射の合計射出）について連立方程式を解きます。

$$ J_i = \varepsilon_i \sigma T_i^4 + (1-\varepsilon_i) \sum_j F_{ij} J_j $$

温度規定面と熱流束規定面（断熱・再放射面は \( q=0 \) の特殊形）が混在してもそのまま解け、解いた後に各面の正味熱流束が得られます。CFD/FEMツールの「S2S（surface-to-surface）」モデルはこの実装です。数値上の主役は方程式ではなく形態係数行列の計算で、面数Nに対して \( N^2 \) 個の幾何積分（遮蔽判定込み）が必要——ここが計算コストと精度の両方を支配します。

形態係数の計算法と検証

どの方法でも、出てきた行列に総和則・相反則の充足度チェックをかけます。ツールは通常、総和則を満たすよう正規化するオプションを持ちますが、正規化前の誤差が大きい（例：総和が0.95未満）ならレイ数・解像度不足のサインで、正規化は誤差を隠すだけです。

面分割の設計——「等温面」仮定を守る

理論の仮定③（各面等温）は、温度勾配の大きい面を1枚として扱わないという分割指針に翻訳されます。長い放熱フィンを1面にすると根元と先端の温度差が潰れ、放射量を誤ります。逆に細かすぎる分割は形態係数計算を爆発させるため、実務では「予想温度差が数十K以内になる単位で分割」が目安です。S2Sのクラスタリング（複数メッシュ面を1放射面に束ねる）設定は、まさにこのトレードオフの調整弁です。

実務適用の指針

手計算例——数字の相場観を持つ

放射シールドの設計

2面の間に低εのシールドを1枚挟むと、放射抵抗が直列に増えて熱流はほぼ半減します（同εなら \( 1/(N+1) \)：N枚で \( 1/(N+1) \) 倍）。多層断熱材（MLI）はこの原理の極致で、宇宙機では数十層で実効放射率0.01級を実現します。設計上の注意は2つ：①シールド効果はεが小さいほど劇的（表面抵抗支配）、②実物ではスペーサ・縫い目の接触伝導と端部リークが理論値を大きく損なう——MLIの実効性能が理論の数分の一になるのは通例で、実測係数（実効放射率）での設計が実務です。

放射の線形化——熱回路・過渡解析への組み込み

温度差が小さい場合、放射は線形化して対流と同じ形で扱えます：\( h_r = 4\varepsilon\sigma T_m^3 \)（\( T_m \)は平均温度）。室温付近では \( h_r \approx 5\varepsilon \) W/m²K——自然対流と同オーダーという相場が即座に出ます。熱回路網・1D解析に放射を安価に組み込む定石ですが、温度差が大きい（比で1.5倍超）場面では線形化誤差が無視できないため、4乗のまま解くか反復更新にします。

ツールでの取り扱い

ツール別の設定要点

開放系のモデル化——「仮想面」を忘れない

エンクロージャ理論は閉空間が前提です。開いた筐体・屋外機器では、開口部に周囲温度の仮想黒体面を張って閉空間化するのが定石で、これを忘れると総和則が破れて熱が消えます。ツールによっては「environment temperature」設定が自動でこの役を担いますが、その周囲温度が何を代表しているか（空・地面・室内壁）は解析者が物理として決める事項です。屋外なら空の実効温度（晴天夜間で外気−10〜−20 K相当）の扱いが結果を左右します。

先端動向

鏡面・スペクトルへの拡張

灰色拡散の先には、部分鏡面反射（拡散+鏡面の分解、さらにBRDFによる一般化）、波長帯域分割（非灰色）、偏光まで含む精密モデルがあります。レイトレーシングのGPU化により、かつて「灰色拡散で我慢」だった大規模モデルでもモンテカルロ分光計算が現実的になりつつあり、ヘッドランプ・光学系・高温炉の設計で採用が進んでいます。まず灰色拡散で全体を掴み、感度の高い面だけ精密化する階層戦略が実務の定石です。

参与媒体（ガス・粒子）との複合

燃焼ガス（CO₂・H₂O）・煤・噴霧が介在する系では、面間交換に媒体の吸収・放射・散乱が加わり、DO/DOM・P1・モンテカルロ等の放射輸送解法とガスの非灰色モデル（WSGG・SLW・k分布）の組み合わせ選定が主題になります。面のS2S理論は「媒体が透明」という近似の明示——炉・エンジン筒内でこの近似が破れることの認識が、モデル選択の出発点です。

近接場放射——理論の外側

面間距離が熱放射の波長（〜10 μm）に近づくと、エバネッセント波のトンネリングで黒体限界を桁で超える熱輸送が起きます（近接場放射）。MEMS・熱光起電力・ナノギャップデバイスで実測が進み、従来の形態係数理論が適用できないスケールの設計論が形成されつつあります。マクロ設計者にとっては「μmギャップでは別の物理」という境界線の知識として押さえておく段階です。

トラブル対応

症状別の原因と対策

放射解析の検証の型

関連記事：キルヒホッフの法則、放射伝熱の記事一覧、NAFEMS T4（境界条件検証の型）。