3. 辐射换热基础与Fluent辐射模型选型

自然对流问题中，当温度水平较高（例如超过300℃）或参与介质具有明显辐射吸收能力时，辐射换热量往往与对流换热量处于同一量级，忽略辐射会导致壁面温度、热流密度和流场结构出现系统性偏差。本节从辐射换热的基本物理量出发，说明Fluent中三种常用辐射模型（DOM、DTRM、S2S）的适用条件与选型依据。

一、辐射换热的基本物理量

黑体辐射是辐射换热的基准。黑体在任意温度下对所有波长、所有方向的辐射能力达到最大值，其光谱辐射力由普朗克定律描述，全波长积分后得到斯特藩–玻尔兹曼定律：

E_b = σ·T^4      σ = 5.67e-8 W/(m²·K⁴)

式中T为绝对温度（K）。注意温度以四次方出现，因此高温工况下辐射项增长极快：温度从300 K升到600 K，黑体辐射力增大16倍，而对流换热系数变化通常远小于此。

发射率ε（黑度）是实际表面辐射力与同温度黑体辐射力之比，0<ε≤1。它同时决定表面的发射能力和吸收能力（基尔霍夫定律：光谱定向发射率等于同条件下的光谱定向吸收率）。工程中常见取值：抛光铝约0.05，氧化铝约0.3～0.5，耐火砖约0.8～0.9，高温烟气中的炭黑颗粒接近0.9以上。Fluent中发射率在壁面边界条件（Thermal选项卡）中设置，是辐射计算最敏感的参数之一。

吸收系数a描述半透明介质（如高温烟气、玻璃熔体、参与性气体）对辐射强度的衰减能力，单位1/m。它与散射系数σ_s共同构成消光系数β=a+σ_s。吸收系数越大，辐射在介质中穿透距离越短，介质越接近“光学厚”。判断光学厚度用τ=β·L（L为特征长度）：τ>>1为光学厚介质，τ<<1为光学薄介质。

散射相函数Φ(Ω'→Ω)描述辐射从方向Ω'散射到方向Ω的概率分布。各向同性散射时Φ=1；前向散射（如大颗粒）Φ在入射方向附近取大值；后向散射则相反。Fluent中可用各向同性、线性各向异性或用户自定义相函数。对于大多数工程燃烧与换热问题，若颗粒浓度不高，常先按各向同性处理，再评估敏感性。

二、三种模型的适用条件与选型逻辑

DOM是通用性最强的模型。它把4π立体角离散为若干方向，对每个方向求解辐射传递方程，因此能处理吸收、发射、散射以及半透明壁面。自然对流中若介质为空气且温度不高，DOM并非必需；但若介质为高温烟气、水蒸气或CO₂，或存在玻璃、石英等半透明窗口，DOM是首选。方向离散数（Theta/Phi Divisions）需权衡精度与成本：3×3适合快速估算，5×5适合大多数工程问题，对强各向异性散射或光学薄介质可提高到8×8以上。

DTRM通过追踪射线积分RTE，概念直观，早期版本中常用于光学厚介质。其缺点是射线数随网格和方向数增加而迅速膨胀，且对散射处理能力弱。在Fluent现代版本中，DTRM已较少作为首选，更多被DOM替代。若必须使用，建议射线数控制在合理范围，并确认介质散射可忽略。

S2S只考虑表面之间的辐射，假设介质不参与辐射（透明）。它通过计算角系数F_ij建立表面间的辐射交换关系，适合封闭腔内的表面辐射问题，例如电子设备机箱内壁面之间的辐射、真空炉内工件与炉壁的换热。S2S不计算介质吸收和散射，因此不能用于烟气、玻璃等参与介质。角系数矩阵规模随表面数平方增长，表面数超过数千时内存和计算时间会显著增加，此时可考虑用DOM替代。

三、Fluent中的关键设置与操作

在Fluent中启用辐射模型：Models → Radiation → Edit，选择模型后需设置以下参数。

对于DOM：Radiation Model选Discrete Ordinates (DO)，设置Theta Divisions和Phi Divisions（默认2×2，建议3×3起步），Energy Iterations per Radiation Iteration通常取1，Radiation Iteration Frequency取1～10。若介质参与辐射，需在Materials中设置Absorption Coefficient和Scattering Coefficient。

对于S2S：选择Surface to Surface (S2S)，需指定参与辐射的壁面，并设置Radiation Boundary Condition。角系数计算在初始化时完成，可通过Compute/Write查看角系数矩阵规模。

壁面发射率设置：Boundary Conditions → Wall → Thermal → Radiation，输入Internal Emissivity。注意区分Internal Emissivity（参与介质侧）和External Emissivity（外部环境侧）。对于不透明壁面，Diffuse Fraction通常取1。

一个典型教学示例：一个1 m×1 m×1 m的封闭方腔，左侧壁面800 K，右侧壁面300 K，其余壁面绝热，内部为空气。若忽略辐射，左侧壁面附近空气温升仅由对流决定；开启S2S后，左侧壁面向右侧壁面的辐射换热使右侧壁面温度明显升高，方腔内温度分布更均匀。若将内部介质改为吸收系数a=1 m⁻¹的烟气，则需改用DOM，此时介质吸收和再发射会进一步改变温度场。该示例参数为合理假设，用于说明模型差异，实际结果需运行后读取。

四、常见错误与判断依据

常见错误包括：在参与介质中使用S2S导致辐射被完全忽略；DOM方向数过低导致辐射热流出现方向性偏差；发射率设置过大或过小使壁面热流偏离实际；散射相函数选择不当使介质内辐射强度分布失真。

判断依据：开启辐射后，若壁面热流密度或温度分布与纯对流结果差异显著（例如超过5%～10%），说明辐射不可忽略；若差异极小，可评估是否关闭辐射以节省计算资源。对于DOM，可通过监测Radiation Heat Flux和Incident Radiation随迭代的变化判断收敛；对于S2S，可检查角系数矩阵是否满足互易性和封闭性。

选型不是一次性的，建议在初步计算后用不同模型或不同方向数做一次对比，确认结果对模型选择不敏感后再固定设置。这样既能保证精度，也能避免不必要的计算开销。