放射ネットワークシミュレーター 戻る
伝熱シミュレーター

放射ネットワークシミュレーター — 表面抵抗と空間抵抗

2つの灰色表面間の放射熱伝達を、電気回路と同じ抵抗ネットワークで可視化。温度と放射率を変えて、なぜ非黒体表面で熱流が落ちるのかを学べます。

パラメータ設定
表面1温度 T_1
K
表面2温度 T_2
K
放射率 ε_1
放射率 ε_2

A_1 = A_2 = 1 m²、形態係数 F_12 = 1(平行無限平板)を仮定しています。

計算結果
放射熱流 Q_12
黒体限界 Q_BB
灰色効率 η = Q/Q_BB
総抵抗 R_tot
放射抵抗ネットワーク

左ノード=σT_1^4 / 右ノード=σT_2^4 / 中間に表面抵抗 R_1, R_2 と空間抵抗 R_12 が直列接続

理論・主要公式

黒体放射に基づく駆動力 σT^4 と、表面抵抗・空間抵抗を直列接続した電気回路アナロジーで放射熱流を計算します。

表面抵抗(放射率 ε < 1 による反射損失):

$$R_s = \frac{1-\varepsilon}{A\,\varepsilon}$$

空間抵抗(形態係数 F_12 による幾何学的損失):

$$R_{12} = \frac{1}{A_1\,F_{12}}$$

2面間の放射熱流。σ = 5.670×10⁻⁸ W/(m²·K⁴):

$$Q_{12} = \frac{\sigma\,(T_1^4 - T_2^4)}{R_1 + R_{12} + R_2}$$

黒体限界(ε_1 = ε_2 = 1)と灰色効率:

$$Q_{BB} = \sigma A (T_1^4 - T_2^4),\quad \eta = \frac{Q_{12}}{Q_{BB}}$$

平行無限平板(F_12 = 1, A_1 = A_2 = A)では η = 1/(1/ε_1 + 1/ε_2 - 1) と簡略化されます。

放射ネットワークシミュレーターとは

🙋
真空中でも熱が伝わるって聞きました。あれが放射ですよね。でも、なぜわざわざ「ネットワーク」って言うんですか?
🎓
そう、放射は電磁波だから真空でも伝わる。ざっくり言うと、放射伝熱の解析を電気回路と同じ形にすると、ものすごく見通しよく解けるんだ。「電圧」が σT^4、「抵抗」が表面抵抗と空間抵抗、「電流」が熱流 Q——この対応関係そのものを「放射ネットワーク」って呼ぶんだよ。
🙋
「表面抵抗」と「空間抵抗」って何が違うんですか?
🎓
表面抵抗 R_s = (1−ε)/(Aε) は、表面が黒体じゃないせいで起きる「反射損失」。放射率 ε が小さい鏡みたいな面ほど、放射エネルギーが反射されて自分に戻ってきてしまう。空間抵抗 R_12 = 1/(A_1·F_12) は、形態係数 F_12 で決まる「幾何学的な届きにくさ」。シミュレーターで放射率スライダーを下げてみて。表面抵抗が跳ね上がって、熱流が一気に減るのが見える。
🙋
温度が4乗で効くってよく聞きます。スライダーを動かすと、Q_12 がすごい勢いで増えますね。
🎓
それがステファン・ボルツマン則 E_b = σT^4 の威力だ。T_1 を800Kから1600Kに上げると、T^4 は16倍。同じ温度差「800K差」でも、800K vs 0K と 1600K vs 800K では桁違いに放射が強い。だから炉や燃焼器みたいな高温機器では、対流より放射が支配的になるんだよ。
🙋
じゃあ、効率 η が灰色だと66.7%しかないのは、どこに「失われた」んですか?
🎓
どこにも消えてはいない。表面抵抗の段で「自分の表面に反射されて戻ってきた」分が、結果的に有効な熱流に寄与しないんだ。だから魔法瓶やMLI(多層断熱)は逆にこれを利用する。ε を極端に小さくしたアルミ箔を何枚も重ねて、表面抵抗を意図的に大きくして放射熱を遮る。シミュレーターで ε_1 = ε_2 = 0.05 にしてみて。Q が一気に1/10以下になるはずだ。

よくある質問

材質と表面状態で決まります。研磨されたアルミは0.04〜0.07、酸化したアルミは0.20〜0.30、酸化した鋼は0.60〜0.85、白色塗料は0.80〜0.95、黒体に近い黒色塗料は0.95以上です。同じ材料でも温度や酸化状態で大きく変わるため、設計時は実測値か信頼できるハンドブックの値を使います。波長依存性が強い場合は分光放射率を考える必要があります。
放射率が波長によらず一定と見なす近似で、解析を非常に簡単にします。多くの工業材料では中赤外域(典型的な工業温度域 300〜2000K の放射ピーク)でほぼ成立します。ただしガラスやセラミックスのように分光特性が強い材料、あるいは選択吸収膜(太陽集熱器の選択面)では、波長帯ごとに分けて扱う必要があります。
温度域で大きく異なります。室温付近では自然対流の熱伝達係数が5〜10 W/(m²·K) 程度、放射の等価熱伝達係数も同程度のため両者は競合します。500K以上では放射が急速に支配的になり、1000Kを超えると対流の数倍〜十数倍に達します。逆に低温では対流・伝導が支配的で、放射は無視されることもあります。
主要な3手法があります。S2S(surface-to-surface)法は形態係数行列を事前計算し本ツールと同じネットワークで解く高速手法、放射部分が透明媒質のときに最適です。DO(discrete ordinates)法は媒質吸収・散乱を扱える有限体積系の手法、燃焼室など参加性媒質に使います。モンテカルロ法は光線追跡で複雑形状や鏡面反射に強く、検証用や高精度計算で使われます。

実世界での応用

炉・燃焼器の設計:製鉄炉・ガラス溶融炉・ボイラ火炉では、放射が伝熱の60〜90%を占めます。バーナーから壁へ、壁から被加熱物へという放射ネットワークを精密に解くことが、燃料消費・温度均一性・耐火物寿命を直接左右します。CFDコードのS2Sやモンテカルロ放射モデルは、まさにこの放射ネットワークを大規模化したものです。

宇宙機の熱設計:真空中の人工衛星では対流が存在せず、ほぼ全ての熱排出が放射です。MLI(多層断熱、放射率0.05のアルミ蒸着フィルム十数層)で機器を覆って表面抵抗を巨大化し、ラジエータ面(高放射率)で宇宙背景へ熱を捨てる設計は、放射ネットワーク解析そのものです。表面処理(白色塗装・銀メッキ・OSR)の選択で熱バランスを制御します。

建物の窓と断熱:Low-E(低放射率)ガラスは、ガラス内面に酸化錫等の薄膜をコーティングして放射率を0.85→0.10程度まで下げ、室内外の放射熱交換を大幅に減らします。複層ガラスでは隙間のアルゴンガスが対流を抑え、Low-Eコーティングが放射を抑えるため、熱貫流率(U値)が単板の1/3以下にまで下がります。

赤外線サーモグラフィの校正:サーモグラフィは表面から放出される放射エネルギーを温度に換算しますが、ε に大きく依存します。同じ温度でも黒色塗料面(ε≈0.95)と研磨ステンレス面(ε≈0.10)では検出強度が10倍近く違うため、測定対象に放射率テープを貼ったり、既知放射率塗料を吹いて校正します。本ツールの ε スライダーを動かしたときの変化が、ちょうどこの誤差源に対応します。

よくある誤解と注意点

最も多い誤解は、「黒い物ほど熱を吸収するから熱くなる」と単純化しすぎることです。確かに可視光域では「黒色=高吸収率」が成り立ちますが、伝熱で重要なのは中赤外域(典型的な工業温度の放射ピーク)の放射率であって、可視色とは一致しません。例えば白色塗料は可視光は反射しますが、中赤外域では ε ≈ 0.9 と「黒体に近い」挙動をします。シミュレーターの ε スライダーは可視色ではなく「対象温度域での放射率」を入れる必要があると覚えてください。

次に多いのが、形態係数 F_12 を「相手が見えていれば1」と勘違いすることです。F_12 は「面1から放出された放射のうち面2に到達する割合」であって、平行無限平板や同心円筒の内側→外側のような完全包囲の場合だけ1になります。一般の凸形状では F_12 < 1 で、残り (1−F_12) は他の面や周囲環境へ向かいます。本ツールでは平行無限平板(F_12 = 1)に簡略化しているので、この場合の単純なネットワーク解が得られますが、実機では形態係数の正確な評価が解析精度を左右します。

最後に、表面抵抗を温度差や絶対温度で割ってはいけない点に注意してください。放射ネットワークは「σT^4 を電位、Q を電流」とするオーム則の類似物で、電位差は σ(T_1^4 − T_2^4) であって T_1 − T_2 ではありません。対流伝熱で慣れた「熱伝達係数 × 温度差」の感覚で R を温度差ベースに置き換えると、低温域では合っているように見えても高温域で大きくずれます。放射の本質は4乗則であり、ネットワークは厳密にこの4乗の差を駆動力とした線形回路として閉じている、という点を押さえましょう。