2. Fluent燃烧传热仿真工作流程
燃烧与传热仿真的难点,往往不在某一个面板怎么点,而在于整条工作链是否前后自洽:几何是否简化得当、网格是否匹配火焰与边界层尺度、求解器与模型是否与流动状态一致、边界条件是否给到了物理上正确的量、后处理是否读到了真正需要判断的指标。本节以一个典型的轴对称旋流燃烧器为例,走一遍从几何到结果判读的完整流程。示例中的尺寸、参数和结果均为教学假设,用于说明操作逻辑,不代表任何实际设备的实测数据。
一、几何处理:先决定“保留什么、简化什么”
示例对象为一个小型天然气旋流燃烧器:中心为燃料喷口,外环为旋流空气通道,下游接一段圆柱形燃烧室,出口为自由出流。几何处理阶段的核心判断有三点。
- 降维判断:燃烧器与燃烧室均为回转体,且旋流器叶片周向均匀分布,可简化为二维轴对称模型。这样网格量可下降一到两个数量级,适合先跑通模型、标定参数。若需要观察火焰偏斜、单喷嘴不对称或旋流进动,则必须回到三维。
- 简化取舍:忽略法兰、螺栓、测压孔、点火器等不影响主流的细节;旋流叶片若只关心下游流场,可用“等效旋流数”代替真实叶片几何,即在入口给定切向速度分量。
- 计算域边界:入口取在喷口上游一段直管处,保证来流均匀;出口取在燃烧室下游足够远处,避免回流区被出口边界截断。
几何处理完成后,需在 Fluent 的 Meshing 或外部网格工具中明确命名边界:fuel_inlet、air_inlet、wall_liner、outlet、axis。命名规范直接影响后续边界条件设置效率,也避免把壁面误设为对称面。
二、网格划分:把分辨率放在梯度大的地方
燃烧与传热对网格的要求集中在三处:火焰锋面附近、壁面边界层、旋流剪切层。示例采用二维轴对称结构化四边形网格,轴向约 400 个节点,径向约 120 个节点,总网格量约 4.8 万。
- 燃料喷口与空气通道出口处加密,径向最小尺寸约 0.2 mm,用于捕捉射流剪切与混合。
- 燃烧室壁面设置 15 层边界层网格,首层高度按 y+ 目标估算。若采用增强壁面处理,目标 y+ 约 1;若采用标准壁面函数,目标 y+ 控制在 30~100。
- 火焰可能驻留的区域(喷口下游 0~0.3 m)做局部加密,远离火焰的下游可适当放粗。
网格质量检查关注三项:最小正交质量、最大长宽比、负体积数量。示例中最小正交质量约 0.62,最大长宽比约 85,无负体积。若正交质量低于 0.1,通常需要回到几何阶段处理尖角或小间隙。
三、求解器与模型选择:先定物理,再定数值
示例为低速、不可压、稳态燃烧问题,选择压力基求解器。若马赫数超过 0.3 或存在强激波,则应改用密度基。湍流模型选择 Realizable k-ε 配合增强壁面处理,兼顾旋流与壁面传热。燃烧模型选择涡耗散概念(EDC)或非预混平衡 PDF:若只关心温度场与主要组分,非预混模型更稳健;若关心 CO、NOx 等慢反应组分,EDC 更合适。辐射模型选择离散坐标模型(DO),因为燃烧室壁面温度较高,辐射换热不可忽略。
四、边界条件设置:把物理量给到正确的位置
示例边界条件如下表所示。质量流量入口比速度入口更适合可压缩或变密度问题,因为燃烧后密度变化大,速度入口难以保证质量守恒。
| 边界 | 类型 | 关键参数 |
|---|---|---|
| fuel_inlet | 质量流量入口 | CH₄ 0.002 kg/s,温度 300 K |
| air_inlet | 质量流量入口 | 空气 0.034 kg/s,温度 600 K,旋流数 0.8 |
| wall_liner | 无滑移壁面 | 温度 900 K,发射率 0.7 |
| outlet | 压力出口 | 表压 0 Pa,回流温度 1200 K |
| axis | 轴对称 | — |
当量比由燃料与空气流量估算:空气质量流量 0.034 kg/s 对应约 0.00117 kmol/s 空气,CH₄ 0.002 kg/s 对应 0.000125 kmol/s,化学计量比下所需空气约 0.00125 kmol/s,因此当量比约为 1.07,属于接近化学计量的贫燃偏富工况。这一判断决定了后续火焰温度预期与 NOx 水平。
五、求解与后处理:从残差到物理判据
求解采用 SIMPLE 压力-速度耦合,压力用二阶格式,动量与能量用二阶迎风,组分用二阶迎风。先以较小松弛因子(压力 0.3、动量 0.5、能量 0.8)跑 500 步,再逐步提高。收敛判据不能只看残差:能量残差降到 1e-6 以下、组分残差降到 1e-5 以下,同时监测出口温度、壁面热流和火焰位置是否稳定。
后处理重点看四类结果:温度云图判断火焰形状与高温区位置;壁面热流分布判断热负荷是否集中;出口组分判断燃烧完全程度;流线图判断回流区是否稳定。若火焰贴壁或出口温度异常偏低,通常需要回到网格或边界条件排查。
整条流程中,最容易出问题的环节往往不是求解器设置,而是几何简化与边界条件之间的匹配:例如把旋流叶片简化为入口切向速度后,若旋流数给得过大,中心回流区会过强,火焰被推向下游,出口温度反而偏低。此时应回到边界条件检查旋流数定义,而不是盲目加密网格。理解每一步的物理含义,比记住面板位置更重要。