ステファン・ボルツマンの法則 — 輻射伝熱の基礎とCAE実装

理論と物理

概要 — なぜ4乗なのか

プランクの分光放射輝度は

これを全波長 $\lambda = 0 \sim \infty$ で積分して全放射エネルギー $E_b$ を求める。変数変換 $x = hc/(\lambda k_B T)$ を使うと、積分が

になる。この定積分は $\Gamma(4)\,\zeta(4) = 6 \times \pi^4/90 = \pi^4/15$ という確定値。だから $E_b \propto T^4$ が出てくる。4乗は「光子の状態密度が3次元空間で $\nu^2$ に比例する」ことと「ボーズ・アインシュタイン分布のエネルギー平均」から来ていて、空間の次元と量子統計の帰結なんだ。

2019年のSI再定義以降、$k_B$, $h$, $c$ はすべて定義値になったから、$\sigma$ も厳密に確定した値だ。CAEの教科書では $5.67 \times 10^{-8}$ と丸めて使うことが多いけど、知っておくと精度議論のときに役立つよ。

支配方程式

ステファン・ボルツマンの法則の基本形は以下の通りである。

実際の物体は黒体ではないため、放射率（emissivity） $\varepsilon$ を導入する：

代表的な放射率を以下に示す。

灰色体の放射交換

これが灰色体の基本式だ。ここで重要なのは $T^4$ の差であって、$(T_1 - T_2)^4$ ではない点。この非線形性がCAEを厄介にする。

例えば、鋳片（$T_1 = 1200°C = 1473\,\text{K}$）と周囲（$T_2 = 30°C = 303\,\text{K}$）で計算してみよう。

一方、同じ条件での自然対流は $h \approx 10$ W/(m²K) として $q_{\text{conv}} = 10 \times 1170 \approx 12$ kW/m²。つまり輻射は対流の約18倍。1200°Cの世界では輻射が圧倒的に支配的なんだ。

多面間の放射交換と形態係数

閉じた空間（エンクロージャ）内での多面間放射交換は、形態係数（View Factor） $F_{ij}$ を用いて記述される。

灰色体・拡散面の場合はラジオシティ法（Radiosity Method）を用いる：

ここで $J_i$ は面 $i$ のラジオシティ（放射＋反射エネルギーの合計）。これはN元連立方程式となり、CAEソルバーでは行列演算で解く。形態係数の性質として以下が重要：

• 相反則： $A_i F_{ij} = A_j F_{ji}$
• 総和則： $\sum_{j=1}^{N} F_{ij} = 1$（閉じたエンクロージャの場合）
• 凸面： $F_{ii} = 0$（凸面は自分自身を見ない）

輻射支配条件 — いつ輻射が主役になるか

ただし、真空中では対流がゼロだから、宇宙機の熱設計では常温でも輻射が唯一の排熱手段になる。ISS（国際宇宙ステーション）のラジエーターは約280K前後で動作していて、$\sigma T^4$ だけで排熱量を決める。温度が低いぶん面積を大きくする必要があるわけだ。

数値解法と実装

T⁴項の線形化

ここで放射熱伝達係数 $h_{\text{rad}}$ は

$h_{\text{rad}}$ は前回の反復での温度を使って計算し、Newton-Raphson法で更新していく。対流の $h_{\text{conv}}$ と足し合わせられるから、既存の熱伝達係数の枠組みに乗せやすい。

これをNR法の接線剛性行列に入れる。Abaqusの*RADIATIONやAnsysのSRDOPTは内部的にこの手法を使っている。収束は速いが、初期温度推定が悪いと発散しやすい。

有限要素法における輻射境界条件

熱伝導方程式の弱形式に輻射境界条件を組み込む。エネルギー方程式は：

輻射面 $\Gamma_r$ 上の境界条件：

Galerkin法による弱形式への組み込み後、要素レベルで以下の非線形残差ベクトルが生じる：

ここで $\mathbf{N}$ は形状関数ベクトル。接線マトリクスは：

連成解法のストラテジー

1. 弱連成（Sequential）: CFDで対流係数 → 熱伝導+輻射をFEMで解く → 温度をCFDにフィードバック。実装は楽だが収束が遅い。
2. 強連成（Monolithic）: 流体・固体・輻射を1つの方程式系にまとめて同時に解く。Ansys FluentのS2S(Surface-to-Surface)モデルはこの方式。収束は速いがメモリ消費大。
3. DO/MCモデル: 関与媒体（燃焼ガス等）がある場合、Discrete Ordinates法やモンテカルロ法で輻射輸送方程式(RTE)を解く。S-Bの法則は壁面境界条件として使われる。

実務では「まず輻射なしで流れ場を収束させ、その後に輻射モデルを有効化」が安定する。いきなり全部オンにすると発散のリスクがある。

実践ガイド

解析フロー

1. 放射面の特定： 高温面と、それが「見える」面を洗い出す。遮蔽の有無も確認。
2. 放射率の設定： 温度依存性があるか確認。酸化状態で大きく変わるから文献値を鵜呑みにしない。
3. 形態係数の計算： ソルバーの自動計算（ヘミキューブ法等）を使うか、解析解があれば検証に使う。
4. メッシュ検討： 輻射面はT⁴の非線形性があるため、温度勾配の急な部分（エッジ、角部）を細かくする。
5. 初期温度の設定： 実際の運転温度に近い値を初期条件にする。300Kスタートだと1500Kまで行くのに反復回数が爆発する。
6. 段階的な解法： 対流→対流+輻射の順で有効化。いきなり全部オンにしない。

事例1: 連続鋳造の放射冷却

1. スプレー冷却帯： 水スプレーによる強制対流（$h = 500\sim2000$ W/(m²K)）が支配的。ここではS-B輻射は全体の10〜20%程度。
2. 放射冷却帯： スプレーが終わった後、鋳片温度はまだ900〜1100°Cある。この区間では$q_{\text{rad}} = \varepsilon\sigma(T_s^4 - T_\infty^4)$ が冷却の80%以上を担う。
3. ロール接触帯： ガイドロールとの接触伝導。接触面積は小さいが局所的に重要。

CAEモデルでは、鋳造方向に沿ったラグランジュ的追跡（スライス法）が一般的だ。鋳片の断面2Dモデルを切り出して、引き抜き速度に合わせて各冷却帯を通過させる。凝固潜熱は $L = 270$ kJ/kg程度をエンタルピー法で処理する。

事例2: 宇宙機の熱設計

宇宙空間では対流がゼロのため、排熱はすべてS-B法則による輻射に依存する。設計上の基本式は：

宇宙空間の輻射温度 $T_{\text{space}} \approx 3\,\text{K}$（宇宙マイクロ波背景放射）は無視できるため、実質的に $\varepsilon\sigma A T_{\text{rad}}^4$ だけで設計する。ISSのラジエーターは $\varepsilon \approx 0.9$、$T_{\text{rad}} \approx 280\,\text{K}$ で運用されている。$280^4 = 6.15 \times 10^9$ に $\sigma$ を掛けると約 $350\,\text{W/m}^2$ — これがラジエーター1m²あたりの排熱能力だ。

太陽入熱（$\alpha_s \cdot S_0 \approx 0.3 \times 1361 = 408\,\text{W/m}^2$）との熱収支から、ラジエーター面積と運用温度が決まる。ここで $\alpha_s / \varepsilon$ 比（太陽光吸収率と赤外放射率の比）が熱設計の最重要パラメータとなる。

事例3: 工業炉の炉壁設計

特にガス炉では燃焼ガス（CO₂, H₂O）自体が放射に関与するため、S-B法則だけでは不十分。Weighted-Sum-of-Gray-Gases (WSGG) モデルで気体の放射率を計算し、輻射輸送方程式を解く必要がある。Ansys FluentのDOモデルやSTAR-CCM+のDO/S2Sモデルがこれに対応している。

ソフトウェア比較

主要CAEツールの輻射機能

オープンソースでの実装

radiationModel  viewFactor;
viewFactorCoeffs
{
    nBands    1;          // 灰色体=1バンド
    smoothing true;       // 形態係数の平滑化
}
emissivityModel lookup;   // テーブル参照

ElmerFEMでも Heat Equation ソルバーに Radiation = Diffuse Gray を指定するだけで灰色体の面間輻射が使える。ただしOSSでは大規模View Factor計算（100万面超）の並列効率に課題があり、商用ツールに比べて10〜50倍遅いことがある。

先端技術

分光放射モデル

灰色体近似（波長に依存しない一定の $\varepsilon$）は多くの実用問題で十分だが、以下の場合には分光（スペクトル）放射モデルが必要になる：

• 選択的放射体： ガラス、セラミックコーティング、半導体ウェハーなど波長による放射率変化が大きい材料
• 太陽光＋赤外輻射の共存： 太陽光（ピーク~0.5μm）と物体の赤外放射（ピーク~10μm）で放射率が大きく異なる場合
• 薄膜干渉： 光学コーティングの熱設計

帯域モデル（Band Model）では波長域を数バンドに分割し、各バンドで異なる $\varepsilon_i$ を使う。Ansys Fluentでは Non-Gray DO モデル、COMSOLでは Multi-band radiation として実装されている。

関与媒体の輻射

高温ガス（CO₂, H₂O, CO, CH₄）は特定の波長帯で輻射の吸収・放射を行う関与媒体（Participating Media）である。この場合、S-B法則は壁面の境界条件としてのみ使用され、ガス内部の輻射輸送は輻射輸送方程式（RTE）で記述される：

ここで $I$ は放射強度、$\kappa$ は吸収係数、$\sigma_s$ は散乱係数、$\Phi$ は散乱の位相関数。壁面での $I_b = \sigma T^4/\pi$ がS-B法則からの寄与である。

GPU加速モンテカルロ法

学術分野では機械学習（PINN: Physics-Informed Neural Networks）で輻射輸送方程式を近似する研究も進んでいる。学習済みネットワークでRTEを高速に解き、燃焼シミュレーションのターンアラウンドを短縮する試みだ。ただし、実務への適用はまだ限定的だね。

トラブルシューティング

よくあるエラーと対策

収束困難への対処

1. 緩和係数（Under-Relaxation）を導入： $T^{n+1} = \omega\,T^{n+1}_{\text{calc}} + (1-\omega)\,T^n$ で $\omega = 0.3\sim0.7$ から始める。
2. 温度クランプ： 1反復あたりの温度更新を最大100K以下に制限する。特に初期数反復で有効。
3. 段階的な有効化： まず $\varepsilon = 0.1$ で解を収束させ、段階的に $\varepsilon$ を実際の値まで上げる（ラムビング）。
4. 初期温度の改善： 定常解析でも、予想される温度分布に近い初期場を与えると反復回数が半減することがある。
5. メッシュ品質： 輻射面のアスペクト比が10を超える要素は温度振動の原因。リメッシュする。

Abaqusでは *STEP, NLGEOM の *CONTROLS で FIELD=TEMPERATURE の許容増分を制限できる。Ansysでは NROPT,UNSYM（非対称NR法）にすると、輻射のヤコビアンが非対称になる問題に対処できるよ。