4. 传热基础与Fluent实现

燃烧仿真的本质，是流动、组分混合与传热三者耦合的结果。火焰向壁面的散热决定了壁温，壁温又反过来影响辐射与对流强度，而温度场又直接改变反应速率。因此在 Fluent 中做燃烧问题，传热设置不是“附加项”，而是与化学反应同等重要的物理约束。本节先把三种传热模式的物理含义与 Fluent 中的实现路径讲清楚，再落到能量方程、壁面热条件与热源项的具体设置上。

一、三种传热模式的物理图像与Fluent实现

导热是分子尺度上的能量传递，由傅里叶定律描述：热流密度与温度梯度成正比，比例系数为导热系数 k。在 Fluent 中，导热不需要单独“开启”，只要能量方程打开，固体域与流体域内部的导热项就自动参与求解。真正需要用户干预的是导热系数的来源：固体壁面材料（如钢、耐火砖）的 k 通常取常数或随温度变化的多项式；流体混合物的 k 则由组分与温度共同决定，燃烧问题中建议使用 kinetic-theory 或 ideal-gas-mixing-law，避免用常数导致高温区导热严重失真。

对流是流体宏观运动携带的能量，由牛顿冷却定律描述：q = h(T_w − T_f)。Fluent 并不直接求解 h，而是求解近壁面的动量与能量边界层，h 是后处理量。因此对流换热的“设置”实质上是边界层网格与湍流模型的设置：壁面附近 y+ 应落在所选壁面函数的适用区间（标准壁面函数约 30~300，增强壁面处理要求 y+≈1），否则算出的壁面热流不可信。对于燃烧室壁面，若使用可缩放壁面函数，建议首层网格 y+ 控制在 30~100 之间。

辐射是高温燃烧中占比最大的传热模式，在 1500 K 以上的火焰中，辐射热流常占总热流的 30%~60%。Fluent 提供多种辐射模型，选择依据如下：光学厚度 aL > 1 的燃烧室优先用 DO（离散坐标）模型，它支持参与性介质、半透明壁面与局部热源；aL < 1 且以气体辐射为主时可用 P-1，计算量小但精度有限；DTRM 适合光学薄介质但计算昂贵；S2S 只适用于无参与性介质的表面辐射，燃烧问题中很少单独使用。燃烧仿真中推荐 DO 模型，并配合灰气体加权和模型（WSGGM）计算气体吸收系数。

二、能量方程的开启与源项设置

在 Fluent 中，能量方程通过 Models → Energy → On 开启。对于燃烧问题，能量方程必须开启，否则组分输运与温度场无法耦合。能量方程中的源项 S_h 包含三部分：化学反应放热、辐射源项、以及用户自定义热源。化学反应放热由组分输运模型自动计算，不需要手动输入；辐射源项在启用辐射模型后自动加入；只有用户自定义热源（如电加热、微波加热、局部点火源）需要手动设置。

设置热源项时，在 Cell Zone Conditions 中选中对应流体域，勾选 Source Terms，在 Energy 一栏填入体积热源值，单位为 W/m³。例如在 0.01 m³ 的燃烧室内模拟 5 kW 的点火加热，热源值应填 5×10⁵ W/m³。若热源随温度或位置变化，可用 UDF 的 DEFINE_SOURCE 宏实现。注意：热源项是体积源，不是面源；若在壁面上加热，应使用壁面热流边界条件而非体积源。

三、壁面热条件的四种类型与选择

Fluent 壁面热条件在 Boundary Conditions → Wall → Thermal 中设置，共四类：

• 定温（Temperature）：直接给定壁面温度 T_w。适用于壁面外有强冷却、壁温基本不变的场景，如带水冷套的燃烧室。示例：T_w = 600 K。
• 定热流（Heat Flux）：给定通过壁面的热流密度 q。适用于已知加热功率的场景，如电加热壁面。示例：q = 50 kW/m²，正值表示热量从外部流入流体域。
• 对流（Convection）：给定外部对流换热系数 h 与外部自由流温度 T_∞。适用于壁面外侧有空气自然对流或强制风冷的场景。示例：h = 15 W/(m²·K)，T_∞ = 300 K。
• 辐射+对流混合（Radiation + Convection）：同时考虑外部辐射与对流。适用于壁面暴露在高温环境中的场景，如炉膛外壁。需设置外部发射率 ε 与外部辐射温度 T_∞。

选择壁面热条件时，一个实用判断是：若壁面外侧有明确冷却介质且换热系数已知，用对流或混合条件；若壁面温度由实验测得且变化不大，用定温；若已知加热功率，用定热流。对于燃烧室壁面，推荐使用“对流+辐射”混合条件，因为它同时考虑了外部冷却与火焰辐射的耦合。

四、壁面热流的分区与耦合设置

实际燃烧室壁面往往不是单一材料，而是由耐火层、保温层、金属外壳组成。Fluent 中可用 Wall 的 Layer 功能设置多层壁面，每层给定厚度与导热系数，Fluent 会自动计算等效热阻。例如：耐火层 0.1 m、k = 1.5 W/(m·K)，保温层 0.05 m、k = 0.1 W/(m·K)，金属壳 0.01 m、k = 45 W/(m·K)。总热阻 R = 0.1/1.5 + 0.05/0.1 + 0.01/45 ≈ 0.567 K·m²/W，等效换热系数约 1.76 W/(m²·K)。这个值可直接用于对流边界条件，也可在 Layer 中逐层设置。

五、一个可复现的示例设置

以直径 0.3 m、长 0.5 m 的圆柱燃烧室为例，甲烷与空气预混燃烧，当量比 1.0，入口速度 5 m/s，入口温度 300 K。传热设置如下：能量方程开启；辐射模型选 DO，吸收系数用 WSGGM；壁面热条件选“对流+辐射”，外部对流系数 h = 20 W/(m²·K)，外部自由流温度 300 K，壁面发射率 ε = 0.8，外部辐射温度 300 K；壁面材料为钢，厚度 0.01 m，k = 45 W/(m·K)。点火阶段在燃烧室中心 0.001 m³ 区域设置体积热源 1×10⁶ W/m³，持续 0.5 s 后关闭。迭代 2000 步后，可观察壁面平均热流约 35 kW/m²，壁面平均温度约 620 K，出口平均温度约 1750 K。这些数值为示例量级，实际结果取决于网格、湍流模型与化学反应机理，需以收敛后的监测曲线为准。

判断传热设置是否合理，可检查三项：壁面热流是否与输入热源量级匹配；能量方程残差是否降至 10⁻⁶ 以下；出口温度是否在绝热火焰温度与入口温度之间。若壁面热流远大于输入热源，说明辐射或对流设置过强；若出口温度接近入口温度，说明能量方程未正确耦合或热源未生效。