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 以下,同时监测出口温度、壁面热流和火焰位置是否稳定。

后处理重点看四类结果:温度云图判断火焰形状与高温区位置;壁面热流分布判断热负荷是否集中;出口组分判断燃烧完全程度;流线图判断回流区是否稳定。若火焰贴壁或出口温度异常偏低,通常需要回到网格或边界条件排查。

旋流燃烧器二维轴对称仿真:几何—网格—边界—结果判读 ① 几何与网格(轴对称剖面) axis 轴对称 fuel_inlet air_inlet 旋流 火焰锋面 wall_liner T=900 K, ε=0.7 outlet p=0 Pa 火焰区局部加密 壁面边界层 ② 关键参数与物理判断 当量比 φ ≈ 1.07(接近化学计量) 旋流数 S = 0.8 → 中心回流区 湍流模型:Realizable k-ε + 增强壁面 燃烧模型:非预混平衡 PDF / EDC 辐射模型:DO,壁面发射率 0.7 压力基求解器,SIMPLE 耦合 ③ 收敛与结果判读 能量残差 1e-7 组分残差 1e-6 出口温度 稳定 壁面热流 稳定 判据:残差达标 + 监测量平稳 + 物理合理 异常:火焰贴壁 / 出口温度偏低 → 查网格与边界 几何简化 → 网格加密 → 模型选择 → 边界条件 → 求解收敛 → 后处理判读
建议在正式计算前,先用简化模型(如非预混平衡 PDF + 粗网格)跑通全流程,确认边界条件与量级合理,再切换到 EDC 或详细化学反应机理进行精细计算。这样能显著减少因模型不匹配导致的反复调试。

整条流程中,最容易出问题的环节往往不是求解器设置,而是几何简化与边界条件之间的匹配:例如把旋流叶片简化为入口切向速度后,若旋流数给得过大,中心回流区会过强,火焰被推向下游,出口温度反而偏低。此时应回到边界条件检查旋流数定义,而不是盲目加密网格。理解每一步的物理含义,比记住面板位置更重要。