分光辐射（光谱辐射）分析

分光辐射理论基础

为什么需要按波长分析

• 太阳电池设计：硅太阳电池在带隙能1.1eV处对应波长约1.1μm。波长大于此值的光无法被激发电子，无法进行能量转换。太阳谱中多少百分比的能量在此波长以下必须通过分光辐射计算。
• 红外传感器·热像仪：中红外（3～5μm）和远红外（8～14μm）大气透过窗口不同。根据目标温度范围选择传感器感应波长段，需要了解对象物的分光辐射特性。
• 选择性吸收面设计：太阳热集热器使用"选择吸收面"，在短波长（太阳光）处高吸收率≈0.95，在长波长（自身热辐射）处低辐射率≈0.05。灰体近似无法设计这种结构。

普朗克辐射定律

其中：

• $C_1 = 2\pi hc^2 = 3.742 \times 10^{8} \ \text{W} \cdot \mu\text{m}^4/\text{m}^2$（第一辐射常数）
• $C_2 = hc/k_B = 1.439 \times 10^{4} \ \mu\text{m} \cdot \text{K}$（第二辐射常数）
• $h$：普朗克常数，$c$：光速，$k_B$：玻尔兹曼常数

维恩位移定律

白炽灯灯丝（2800K）峰值在1μm附近的近红外。可见光范围（0.4～0.7μm）只占全能量的约10%。这就是白炽灯电能的90%被"浪费"成热的原因——LED的优势正在于此。

波段辐射率 F(0→λT)

$F$ 的取值范围是0到1，作为 $\lambda T$ 的函数是一一对应的。波长 $\lambda_1$ 到 $\lambda_2$ 范围内的辐射能比例为：

太阳光（5778K）的可见光范围（0.4～0.7μm）能量占比为 $F(0 \to 0.7 \times 5778) - F(0 \to 0.4 \times 5778) = F(0 \to 4045) - F(0 \to 2311)$ ≈ 0.491 − 0.124 = 0.367（约37%）。

分光辐射率与实际表面

全波长的全辐射率 $\varepsilon(T)$ 是分光辐射率对普朗克分布的加权平均：

一个例子：玻璃在可见光(0.4～0.7μm)处 $\varepsilon_\lambda \approx 0.05$（几乎透明），但在红外(5μm以上)处 $\varepsilon_\lambda \approx 0.95$。这正是温室效应的物理根源——太阳短波光通过玻璃，但室内物体的长波红外被玻璃吸收再辐射回室内。

分光辐射数值计算方法

波长波段分割法（Band Model）

分割的要点：

• 辐射率变化剧烈的波段要细分
• 普朗克分布能量集中的波段（峰值附近）要细分
• 能量极少的波段可粗分或忽略

加权和灰气体模型（WSGGM）

每个灰气体 $j$ 有吸收系数 $\kappa_j$ 和温度相关权重 $a_j(T)$：

这里 $p$ 是气体分压，$s$ 是光路长。$j=0$ 对应透明窗（$\kappa_0=0$）。

蒙特卡洛法结合

优点是避免了波段分割的离散化误差。缺点是需要大量光子（通常10⁶～10⁸）才能降低统计噪声，计算成本高。多用于空间飞行器热设计（如Thermal Desktop）和炉内辐射的参考解。

数值积分实现

1. 级数求和：用前面的无穷级数，求和到 $n=10$ 左右。当 $\xi = C_2/(\lambda T) > 2$ 时收敛快，最容易实现。
2. 数值积分：台形法或辛普生法直接对普朗克式积分。用对数尺度采样波长效率更高。
3. 表查补插：预制 $\lambda T$ 对应 $F$ 的表格，线性或3次样条补插。教科书的放射函数表就是这类。

实务中通常用方法1，精度能达 $10^{-8}$ 相对误差。

分光辐射工程应用

分光辐射率数据获取

1. 文献数据库：NASA技术报告的"Thermal Radiation Properties Survey"，Touloukian等的《Thermophysical Properties of Matter》系列（全13卷）是经典。几千种金属、陶瓷、涂料的分光辐射率都在里面。
2. FTIR实测：傅里叶红外分光（FTIR）可测2.5～25μm波段，分辨率0.1cm⁻¹。自己的样品、不同表面处理、温度影响等细节必须自测。
3. 厂商数据：涂料厂（如Acktar超黑涂料）或光学镀膜厂有时公布分光反射率数据，用 $\varepsilon_\lambda = 1 - \rho_\lambda$（不透明体）转换。

铝的抛光面(ε=0.04)和阳极氧化面(ε=0.80)相差20倍！说"铝的辐射率"而不指定表面处理，会导致设计错误。

分析流程

1. 把握对象温度范围：用维恩定律估算关键波长域。300K～1500K的场景，2～10μm带占主导。
2. 获取分光辐射率数据：所有表面的 $\varepsilon_\lambda(\lambda)$，从文献或实测入手。波长覆盖范围至少从普朗克分布 $E_{b\lambda}$ 达到1%的短波到衰减到1%的长波。
3. 确定波段分割：按辐射率变化点、大气透过窗、传感器灵敏波段等确定分割。
4. 各波段灰体辐射伝热计算：每个波段作灰体处理，积分所有波段结果。
5. 结果验证：用全辐射率的灰体解析解对比，与实测红外像仪数据交叉核对。

常见失误与对策

分光辐射软件对比

商业工具的分光辐射支持

• 太阳谱带(0.2～3μm)：直射日光+地球反照率。用太阳吸收率 $\alpha_s$
• 星球红外带(3～50μm)：地球或月球红外放射。用星球红外反照率和放射温度表征
• 深宇宙带(全波）：宇宙背景辐射(2.7K)。基本可忽略，极低温部件才考虑

太阳电池面用 $\alpha_s / \varepsilon_{\text{IR}}$ 比（太阳吸收率/红外辐射率）来表示特征。这种二波段设计正是分光辐射应用的典型。

开源方案

• OpenFOAM（fvDOM模型）：最多4波段波段分割。波段数少，但可自定义扩展。
• Py4CAtS / RADIS：Python红外辐射传递代码。和HITRAN数据库联动，Line-by-Line精确计算气体分光吸收放射。可与CFD耦合。
• libRadtran：大气辐射传递代码。特化于太阳和地球辐射谱计算。太阳电池日射谱推算有用。

分光辐射前沿研究

机器学习预测分光特性

1. 分光辐射率预测：用神经网络从材料组成、微观结构参数推算 $\varepsilon_\lambda(\lambda)$。以前每个样品要FTIR测，现在可快速探索材料设计空间。
2. WSGGM高精化：深度学习对Line-by-Line结果作代理模型，温度、压力、气体组成变化时实时输出最优WSGGM系数。比传统固定多项式拟合覆盖范围广。
3. 反演问题（推断分光辐射率）：从多波长辐射计测量数据同时推断温度和 $\varepsilon_\lambda$。贝叶斯推论和PINN（物理信息神经网络）有应用。

超材料与选择性辐射面

2014年斯坦福Raman团队用光子晶体结构实现白天辐射冷却突破。目前研发方向包括：

• 纳米薄膜：多层膜和超材料精密设计 $\varepsilon_\lambda(\lambda)$
• 涂料型：普渡大学钡硫酸涂料（可见光反射率98.1%）商用化进行
• 纺织品：可穿戴的辐射冷却纤维

这些全依赖 $\varepsilon_\lambda(\lambda)$ 的精心设计，分光辐射分析是核心技术。

分光辐射分析故障排除

分光辐射分析典型问题

1. 波段辐射率合计与全辐射率对不上

• 现象：$\sum \varepsilon_i \cdot \Delta F_i \neq \varepsilon_{\text{总}}$（文献值）
• 原因：波段分割过粗，或波段边界与辐射率变化点不对齐
• 对策：增加波段数看收敛。叠加 $\varepsilon_\lambda(\lambda)$ 图和波段边界，确认变化点被波段内部包含

2. 气体WSGGM温度范围超限错误

• 现象：高温区(>2500K)辐射过高或过低
• 原因：WSGGM系数的拟合温度范围被超过
• 对策：用宽温度范围系数（如Kang等2025版400～3000K）或改用CKD/SLW法

3. 结果与灰体解析相差很大

• 现象：分光和灰体温度差>50K
• 原因：不一定是bug。选择性面或多温度源混合时，这偏差是物理正常
• 对策：先做单温单面简单验算（$q = \varepsilon(T) \sigma T^4$），确认代码逻辑无误再进复杂问题

4. 太阳照射条件的光谱输入不对

• 现象：太阳照射下温度过高或过低
• 原因：AM0(太空1361 W/m²)和AM1.5(地表1000 W/m²)混淆，或忽视大气散射
• 对策：用ASTM G173-03的AM1.5标准谱。区分直达光和散射光

5. 计算时间爆炸

• 现象：波段数×空间网格×角度离散=计算时间难以承受
• 对策：最小化波段（只保留物理关键波段）。DO法的角离散从θ×φ=3×3开始，逐步增加做收敛性检查