参与媒质的辐射传递（Participating Media Radiation）
参与媒质的辐射传递（Participating Media Radiation）的理论基础
什么是参与媒质
老师，我第一次听到"参与媒质"这个词。普通空气也会对辐射产生影响吗？
好问题。氮气 (N₂) 和氧气 (O₂) 这样的对称二原子分子实际上几乎不吸收也不放出红外线。所以常温空气对辐射来说基本上是"透明的"。
但是 CO₂ 和 H₂O 这样的多原子分子就不一样了。它们在分子振动的固有振动频率对应的特定波长带（2.7μm、4.3μm、15μm 等）中强烈吸收红外线，同时也放出。这种"与辐射能相互作用的媒体"叫做参与媒质（participating media）。
但在日常生活中，我们感觉不到空气在吸收红外线啊？什么时候这才变得重要呢？
关键在于温度和尺度。常温下相隔1米的气体辐射基本可以忽略。但在以下情况，情况就完全不同了：
- 燃气轮机燃烧室（1500°C以上）：燃烧气中的CO₂和H₂O辐射占壁面热流通量的30～50%
- 工业炉、锅炉（1000～1600°C）：炉内高温气体到加热物体的主要传热方式变成辐射
- 大规模火灾：火焰的辐射热是周围蔓延的原因。煤烟粒子的贡献很大
- 大气辐射：地球温室效应本身就是CO₂·H₂O的红外线吸收现象
简单说，高温气体存在的系统中忽略辐射会导致严重的设计错误。
哇，壁面热流通量的30～50%来自气体辐射！这不能忽略。普通的表面间视角因子就不够了吗？
完全正确。表面间辐射用视角因子 $F_{ij}$ 就可以只计算表面之间的交互，但参与媒质时媒体本身既是发光体又是吸收体。光线通过媒体时每次都有能量变化。所以必须求解辐射传递方程（RTE：Radiative Transfer Equation）这个积分微分方程。
辐射传递方程（RTE）
听到RTE就觉得很难…是什么方程呢？
RTE是"光线沿某个方向传播时能量强度如何变化"的方程。定常状态下的形式如下：
来整理一下各项的意义：
- $I(\mathbf{r}, \hat{\mathbf{s}})$：位置 $\mathbf{r}$、方向 $\hat{\mathbf{s}}$ 的辐射强度 [W/(m²·sr)]
- $\kappa$：吸收系数 [1/m] — 媒体吸收辐射的程度
- $\sigma_s$：散射系数 [1/m] — 媒体散射辐射的程度
- $I_b = \frac{\sigma T^4}{\pi}$：黑体辐射强度（普朗克函数的方向积分）
- $\Phi(\hat{\mathbf{s}}' \to \hat{\mathbf{s}})$：散射相位函数 — 从方向 $\hat{\mathbf{s}}'$ 散射到 $\hat{\mathbf{s}}$ 的概率分布
左边是"光线传播距离 $ds$ 时强度的变化"，右边的4项对应"放出增加""吸收减少""向其他方向散射减少""从其他方向散射到达增加"。
4项各自代表增减。那这个方程为什么比普通微分方程难解呢？
RTE难解有3个原因：
- 7维问题：位置3维 × 方向2维 × 波长1维 × 时间1维。加上时间和全波谱，计算量巨大
- 积分微分方程：散射项含有全立体角 $4\pi$ 的积分，各个方向的解互相耦合
- 非线性耦合：通过放出项 $I_b \propto T^4$ 与能量方程非线性耦合
所以RTE的计算成本与Navier-Stokes方程相当，有时甚至更大。
从RTE导出的重要关系是辐射热源项（加入能量方程的项）：
这里 $G = \int_{4\pi} I \, d\Omega$ 是入射辐射，$\sigma$ 是Stefan-Boltzmann常数（$5.67 \times 10^{-8}$ W/(m²·K⁴)）。通过这一项，辐射与传导、对流耦合。
光学性质与光学厚度
"光学厚度"这个词我经常看到，具体是什么意思呢？
光学厚度（optical thickness） $\tau$ 表示媒体衰减辐射的程度，是无量纲数：
其中 $\beta = \kappa + \sigma_s$ 是消光系数（extinction coefficient），$L$ 是特征长度。
光学厚度的值决定了媒体行为变化很大：
| 光学厚度 $\tau$ | 分类 | 物理行为 | 解法指导 |
|---|---|---|---|
| $\tau < 0.1$ | 光学上薄 | 辐射基本透过媒体。吸收、放出贡献小 | 表面间辐射＋薄气体修正可应对 |
| $0.1 < \tau < 3$ | 中间区域 | 吸收、放出、散射都重要。最难分析 | 需要DOM（S4以上）或蒙特卡罗法 |
| $\tau > 3$ | 光学上厚 | 辐射被短距离吸收和再放出。扩散性行为 | P1近似或Rosseland扩散近似有效 |
例如，燃气轮机燃烧器内的CO₂丰富燃烧气 $\tau \approx 5\sim10$，属光学上厚。而大气中1米距离 $\tau \ll 0.1$，属光学上薄。
光学厚度不同就能用的近似法不同啊。中间区域最复杂，我能理解。
气体辐射物性模型
老师，CO₂和H₂O的吸收系数随波长变化很大对吧？怎么处理呢？
这是气体辐射中最棘手的问题。实际的CO₂、H₂O吸收谱由数十万条谱线组成，结构极其复杂。波长每隔0.01cm⁻¹吸收系数就变化好几个数量级，全谱直接计算的Line-by-Line (LBL) 法需要数万～数十万个波长格子，工程计算不现实。
所以实务上用以下带模型或全局模型：
WSGGM（加权灰气体模型）
最广泛使用的模型。用 $N$ 个虚拟灰气体（gray gas）的加权和近似实际气体的辐射特性：
其中 $a_{\varepsilon,i}(T)$ 是温度依赖的加权系数，$\kappa_i$ 是第 $i$ 个灰气体的吸收系数，$pL$ 是压力×光路长。$i=0$ 对应"透明窗"，令 $\kappa_0 = 0$。
WSGGM的优点是，通常3～5个灰气体就能获得工程精度（与LBL相比误差5%以内）。代表性的加权系数是Smith et al. (1982) 提出的，几乎所有CFD求解器都默认搭载。
但也有注意事项：
- 基于均一温度、均一组成的气体柱假设 → 非均一场应用需要修正
- CO₂/H₂O摩尔比变化的系统（煤燃烧 vs. 气燃烧）需要换系数
- 高压下（10 atm以上）需要考虑压力增宽效应的系数
WSGGM以外还有其他方法吗？
当然有。根据精度和计算成本的权衡来选择：
| 模型 | 精度 | 计算成本 | 适用场景 |
|---|---|---|---|
| Line-by-Line (LBL) | 最高（参考解） | 极高 | 基准验证用 |
| SNB（统计窄带） | 高 | 高 | 研究用详细分析 |
| WSGGM | 工程上充分 | 低 | 工业CFD标准 |
| 灰气体（Gray Gas） | 低 | 最低 | 概念设计阶段 |
煤烟粒子的辐射
燃烧分析中经常听到"煤烟（soot）"的话题。有煤烟时辐射会怎么变呢？
从辐射的角度看，煤烟是非常特殊的存在。气体分子仅在特定波长带吸收、放出，但煤烟粒子在连续波谱全域都与辐射相互作用。也就是说，煤烟表现为基本灰色的参与媒质。
煤烟的吸收系数近似为：
其中 $f_v$ 是煤烟体积分数（典型 $10^{-7} \sim 10^{-5}$），$T$ 是温度 [K]，$C$ 是常数（约 1860 m⁻¹K⁻¹）。总吸收系数为 $\kappa_{total} = \kappa_{gas} + \kappa_{soot}$ 相加。
例如体积分数 $f_v = 10^{-6}$、温度1500K的煤烟，$\kappa_{soot} \approx 2.8$ m⁻¹，与CO₂的吸收系数相当。少量煤烟对辐射的影响也很大。柴油发动机燃烧室或池火中，煤烟辐射往往占主导。
Hottel的气体辐射图表 — 辐射工学的起点
参与媒质辐射研究始于1920年代MIT的Horace Hottel。他通过实验将CO₂、H₂O的全放射率（total emissivity）整理为温度和压力×光路长（pL）的函数图表。数字计算机出现之前，工程师只靠这张图表和手算就设计工业炉。1982年Smith、Shen、Friedman提出的WSGGM是对Hottel图表的数式拟合，使CFD求解器可以调用，是划时代的工作。至今WSGGM的精度仍通过与Hottel图表的一致性来验证。
参与媒质的辐射传递（Participating Media Radiation）的数值计算方法
离散坐标法（DOM / Discrete Ordinates Method）
RTE实际在计算机上怎样求解呢？方向积分看起来很复杂啊…
最广泛使用的方法是离散坐标法（DOM：Discrete Ordinates Method）。思路很简单，把连续的方向 $\hat{\mathbf{s}}$ 用有限个离散方向 $\hat{\mathbf{s}}_m$（$m=1, \ldots, M$）代表，对每个方向分别求解RTE。
DOM中，RTE对每个离散方向 $\hat{\mathbf{s}}_m$ 近似为：
其中 $I_m = I(\mathbf{r}, \hat{\mathbf{s}}_m)$，$w_{m'}$ 是方向 $m'$ 的求积权重。方向选择用S$_N$ 近似广泛应用，$N$ 值决定了精度：
| S$_N$ 次数 | 方向数（3D） | 精度 | 计算成本 | 应用指导 |
|---|---|---|---|---|
| S2 | 8 | 低 | 低 | 预备计算、趋势把握 |
| S4 | 24 | 中 | 中 | 一般工业计算最低标准 |
| S6 | 48 | 高 | 高 | 标准详细分析 |
| S8 | 80 | 非常高 | 非常高 | 光学上薄区域精密分析 |
方向增多精度上升，但计算成本也上升。实务上S4是标准吗？
完全对。实务上S4～S6最常见。S2的方向分辨率太低，容易出现射线效应（ray effect）的人工条纹。反之S8以上计算成本跳跃式增加，只有特别需要精度时才用。
DOM的优点包括：
- 对光学厚度无依赖（从薄到厚都能用）
- 散射强的系统也能应对
- 与非结构网格兼容性好
- 几乎全部CFD求解器都搭载
缺点是"方向数 × 空间网格数"的联立方程组，内存消耗大。
P1近似（球面谐函数法）
DOM以外还有什么方法吗？听说过P1近似这个说法。
P1近似是球面谐函数展开（P$_N$展开）在1阶截断的方法。将辐射强度 $I$ 的方向依赖性展开：
应用这个近似，RTE变成入射辐射 $G$ 的椭圆型偏微分方程：
看这个式子，形状完全像热传导方程对吧？所以P1近似又叫"辐射扩散近似"。已有的FEM或FVM框架就能求解，实现容易。
P1近似的特点总结为：
- 优点：计算成本极低。比DOM少1/10～1/50。只需求一个标量方程
- 缺点：方向依赖性仅在1阶截断，辐射强度异向性强的系统（光学上薄区域、局部热源）精度下降
- 适用范围：$\tau > 1$（光学上厚区域）效果好。光学上薄系统误差大
光学上厚区域能用，薄区域不行。燃烧炉内部厚，出口附近薄，所以要注意…对吧？
完全正确！实务中常见的做法是，先用P1做预备计算摸清全貌，最终精密计算时换成DOM。Ansys Fluent的GUI上很容易切换。
离散传递法（DTM）与蒙特卡罗法
还有其他解法吗？
还有2个重要方法：
离散传递法（DTM：Discrete Transfer Method）
从每个壁面单元发射代表光线，光线穿过各个单元时分析地积分RTE。计算量由壁面单元数 × 光线数决定。比DOM的内存消耗小，但不适合散射强的系统。STAR-CCM+在"Surface-to-Surface"模式中提供DTM。
蒙特卡罗法（MCM）
大量随机发射光子束，统计再现辐射场的方法。理论上RTE的严密解，可信度高。但计算成本巨大（需数百万光子），工业计算中很少用。OpenFOAM的fvDOM或OpenRad等可用。
解法选择指南
结局就是，应该选哪个方法？判断基准教我一下。
实务的选择流程是这样的：
- 光学上厚（$\tau > 3$）＋无散射 → P1近似或Rosseland扩散近似
- 中间～薄（$\tau < 3$） → DOM（S4～S6）
- 散射重要（煤烟、粉尘） → DOM（含散射相位函数）
- 超复杂形状＋高精度要求 → 蒙特卡罗法
- 预备计算、快速掌握趋势 → P1近似
现场最常见是DOM的S4。这对大多数工业燃烧问题都够用。
与能量方程的耦合
RTE求了之后，怎么反馈给温度场啊？
在能量方程里加辐射热源项 $-\nabla \cdot \mathbf{q}_r$：
$\nabla \cdot \mathbf{q}_r = \kappa(4\sigma T^4 - G)$ 所以，$G > 4\sigma T^4$ 的单元（外来辐射多）温度升高，反之降低。$T^4$ 的非线性在这里出现，温度微小变化也会大幅改变辐射热源，因此需要迭代收敛。
实务上，CFD求解器的每个迭代步：
- 从现在的温度场计算 $I_b = \sigma T^4/\pi$
- 求解RTE得到辐射场（$G$ 或 $\mathbf{q}_r$）
- 将 $-\nabla \cdot \mathbf{q}_r$ 代入能量方程更新温度
- 反复1～3直到收敛
Fluent等求解器默认10步迭代仅解一次辐射，降低成本。但强辐射耦合的系统（如燃气轮机燃烧器）每步都解会更稳定。
参与媒质的辐射传递（Participating Media Radiation）的实务应用
分析流程
从零开始做包含参与媒质辐射的CFD分析，怎样的步骤呢？
标准的工作流如下：
- 问题评估：辐射重要吗？（温度1000°C以上？CO₂/H₂O/煤烟存在？）
- 光学厚度估算：$\tau = \kappa L$ 概估 → 选解法（DOM/P1）
- 气体辐射物性模型选择：WSGGM（标准）、灰气体（简易）、SNB（高精度）
- 网格生成：壁面附近充分细化（与壁面热流通量精度直接相关）
- 边界条件设定：壁面放射率 $\varepsilon_w$、温度、气体组成
- 无辐射模型预备计算：先稳定流场和温度
- 追加辐射模型重算：P1粗→DOM精
- 网格感度、方向感度确认：S4→S6结果不变？
- 后处理：壁面热流通量（对流+辐射）、温度分布确认
网格设计指南
辐射分析中网格有特殊要求吗？
几个重要注意事项：
- 壁面附近分辨率：辐射热流通量在壁面单元评估。需要与边界层相当的精度
- 与光学平均自由程 $\ell = 1/\beta$ 的关系：网格大小超过 $\ell$ 就无法适当解辐射场。光学上厚系统（$\ell$ 小）网格要求宽松
- DOM射线效应对策：低阶S$_N$ 近似用粗网格会出数值条纹。提高网格密度或增加S$_N$ 次数
- 火焰面分辨：温度梯度陡峭的火焰面，$T^4$ 非线性会放大网格粗度导致的温度误差
边界条件的设置
辐射的边界条件，与普通热分析有什么不同吗？
壁面的辐射边界条件如下：
重要参数：
- 壁面放射率 $\varepsilon_w$：氧化金属0.7～0.9、耐火砖0.6～0.8、抛光金属0.05～0.2。现场要注意经年氧化层增厚导致数值变化
- 入出口辐射边界：来自外部的放射。燃气轮机高温燃烧气入口放射强度设置错误会导致壁面热流低估
- 半透明壁面（玻璃、石英）：某些波长透射。需要波长依赖透射率模型
现场最常见错误是壁面放射率随意设为1.0。实际金属壁约0.3～0.8，与1.0的差异能导致壁面热流通量20～40%的差异。
实务技巧
实务上要特别注意的地方教我！
实战工程师常常卡壳的地方及其对策整理如下：
| 常见问题 | 原因 | 对策 |
|---|---|---|
| 壁面热流仅为实验值的一半 | 忽略气体辐射或WSGGM组成比用默认值 | 激活辐射模型，正确设CO₂/H₂O摩尔比 |
| 应该有局部温度极值的地方反而平坦 | P1近似的适用极限（光学上薄区域） | 换DOM。最少S4 |
| DOM壁面出现条纹 | 射线效应（低阶S$_N$ + 粗网格） | 提高S$_N$ 次数或细分网格 |
| 带辐射的计算发散 | $T^4$ 非线性导致的发散 | lower-relaxation系数0.8→0.5。先无辐射求稳定解 |
| 含煤烟火焰温度过高 | 未启用煤烟模型 | 输入煤烟体积分数，追加 $\kappa_{soot}$ |
壁面放射率的设置错误会导致20～40%的差异，很可怕。没有实测数据的时候怎么办？
文献值（Incropera等传热工学教科书）引用典型值后做灵敏度分析是铁则。放射率±0.1地变化，看结果变化多大。如果壁面热流通量对放射率敏感，就要要求实测或采用保守数值。
参与媒质的辐射传递（Participating Media Radiation）的软件比较
Ansys Fluent
首先Ansys Fluent中的设置方法教一下。
Fluent是辐射模型支持最充实的求解器之一：
- 可用模型：DO (S2～S8以上)、P1、DTM (DTRM)、S2S (Surface-to-Surface)、蒙特卡罗
- 气体物性：WSGGM (Smith et al. 默认)、Gray、EWBM (Exponential Wide Band Model)
- 设置路径：Models → Radiation → Discrete Ordinates → Angular Discretization (Theta/Phi Divisions)
- 煤烟模型：1-step soot model、2-step Tesner-Magnussen model
- 实务Tips：Theta Divisions = Phi Divisions = 2 相当于S4。
solve/set/discretization-scheme/do-intensity 1切换到2阶精度
Simcenter STAR-CCM+
STAR-CCM+怎么样？
STAR-CCM+的Radiation模型有以下：
- DOM：S2～S8（用Quadrature Level指定）
- S2S (Surface-to-Surface)：无参与媒质时高速
- 气体物性：WSGGM (Coppalle and Vervisch模型也支持)
- 设置路径：Physics → Radiation → Participating Media Radiation → DOM
- 多面体网格相容性：DOM的空间离散对多面体网格优化。比四面体网格射线效应少
OpenFOAM
开源的OpenFOAM中怎么做呢？
OpenFOAM也有辐射模型实装：
- fvDOM：有限体积版DOM。
constant/radiationProperties配置 - P1：
radiationModel P1;激活 - 气体物性：
greyMeanAbsorptionEmission（灰色）、wideBandAbsorptionEmission（宽带）。WSGGM无默认，但用户库可追加 - 注意事项：
Qr（辐射热流通量）字段要在能量方程中参照。否则辐射不反映到温度
OpenFOAM设置自由度大但参数组合易出错。商用求解器的GUI相比学习成本高，但自定义模型实装的强度是优势。
其他工具
| 软件 | 辐射模型 | 气体物性 | 特色 |
|---|---|---|---|
| COMSOL Multiphysics | DOM, P1 | WSGGM, Gray | 多物理场耦合容易。Heat Transfer Module中激活辐射 |
| Ansys CFX | DTM, 蒙特卡罗 | WSGGM | Fluent的DOM与不同，基于DTM方法 |
| CONVERGE | DOM, Gray | WSGGM | 发动机燃烧专化。与AMR（自适应网格细分）联动 |
参与媒质的辐射传递（Participating Media Radiation）的前沿研究
机器学习加速辐射模型
最近的研究中怎样的进展？
最热的话题是机器学习（ML）代替辐射模型（代理模型）。RTE计算占全计算时间的30～70%，这部分高速化影响很大。
- 神经网络代替WSGGM：温度、压力、组成→吸收系数、加权系数的映射用DNN学习。用WSGGM成本实现LBL精度
- 直接预测辐射热源项：温度场→$\nabla \cdot \mathbf{q}_r$ 用CNN直接预测。无须求RTE
- 物理约束神经网络 (PINNs)：RTE本身作损失函数组入的物理约束神经网络
但学习数据范围外（外推）的信赖性是问题，2025年现在还处研究阶段。实务采用有限。
光子蒙特卡罗法的进展
GPU性能向上，GPU并列光子蒙特卡罗法渐成实用。従来CPU版需数日的计算，NVIDIA A100等GPU数小时内完成的报告有。这样蒙特卡罗法从"基准专用"变成"实务计算可用"。
多尺度辐射分析
另一个重要方向是多尺度方法。从分子级光谱数据库（HITRAN/HITEMP）直接计算吸收系数，再传给CFD辐射模型。离开传统的WSGGM这样的经验拟合，实现第一原理的辐射分析变可能。特别是H₂燃烧（零碳燃烧）这样新燃烧技术，既有WSGGM系数不可用，这方法不可欠。
水素燃焼CO2出既存的WSGGM使的盲点！ 脱炭素時代放射模型変。
参与媒质的辐射传递（Participating Media Radiation）的故障排除
收敛问题与对策
加入辐射后计算根本不收敛…？
辐射收敛问题是实务中常见的。原因与对策整理：
1. 初期化问题
温度场为零或非现实值时 $T^4$ 会暴走。先无辐射稳定流场和温度，然后加辐射是铁则。
2. Under-Relaxation调整
辐射under-relaxation（松弛系数）从默认1.0降至0.5～0.8。Fluent中 solve/set/under-relaxation/do-intensity 0.7 设置。
3. 辐射更新频度
每步求RTE成本高。Fluent默认10步求一次，但强辐射耦合时5步求一次会改善收敛。
4. $T^4$ 线性化
放出项 $\kappa I_b \propto T^4$ 用Taylor展开线性化有效：
这样线性化能强化能量方程系数矩阵的对角项，稳定性上升。大多数CFD求解器自动处理，但定制代码需要手工实装。
常见错误
初心者特別是容易陷入的陷阱是什么？
看过多次的典型错误列举：
- 温度单位混淆：辐射与 $T^4$ 成正，摄氏度(°C)和绝对温度(K)混淆差异大。300°C→573K vs 300K辐射差12倍
- "空气透明"的刻板印象：高温燃烧气的CO₂、H₂O贡献不能忽视。温度1000°C以上必考虑参与媒质
- WSGGM默认系数不变：Smith et al.系数以CH₄/空气燃烧为前提。煤燃烧、富氧燃烧的CO₂/H₂O摩尔比不同，需换系数
- 壁面放射率随意设1.0：实际工业材料0.3～0.9。特别不锈钢抛光面 $\varepsilon \approx 0.15$
- 方向分辨率不足：S2"凑合能跑"后本计算就用它→射线效应结果非物理
- 散射常设为零：含煤烟、飞灰系统散射改变辐射场。煤燃烧需考虑flyash散射
验证基准
自己的分析结果检证用的基准问题有吗？
辐射分析常用的基准问题介绍：
| 基准名称 | 概要 | 特点 | 参考文献 |
|---|---|---|---|
| 1D平板间辐射 | 2个无限平板间灰气体 | 有分析解。光学厚度影响检证 | Modest (2013) 教科书 |
| 正方形围壳体 | 4壁面围参与媒质 | DOM/P1/MC精度比较常用 | Kim & Baek (1991) |
| Sandia Flame D | 甲烷/空气拡散火焔 | TRI（乱流-辐射相互作用）检証 | TNF Workshop |
| IFRF 1MW炉 | 工業用燃烧炉実測数据 | 有温度、热流通量詳細计测 | IFRF Database |
1D平板（有分析解可比对）始自模型设定正确认後，複雑問題進良手順。
参与媒质的辐射传递，真的非常深奥啊…。从RTE各项物理意义到实务中的设置方法，整体脉络都清楚了。首先从1D平板基准开始试试！
好的，好的方法。辐射不是"加上去精度就上升"的东西，正确的模型和适当的参数才有价值。光学厚度估算→解法选择→气体物性模型→壁面放射率，这4步要始终意识到。加油！
更详细
错误