1. 喷雾燃烧仿真技术概述
在工业燃烧器和航空发动机领域，喷雾燃烧是最常见的燃烧形式之一。当液体燃料通过喷嘴雾化成微小液滴并与氧化剂混合燃烧时，会形成复杂的多相湍流燃烧过程。这种燃烧方式在燃气轮机、柴油机和火箭发动机中广泛应用，其燃烧效率直接影响设备性能和排放水平。
我从事燃烧仿真工作已有八年，处理过上百个喷雾燃烧案例。实际工程中最头疼的就是液滴破碎、蒸发和燃烧的耦合计算，这些过程发生在毫秒级时间内，却需要同时考虑流体力学、热力学和化学反应动力学。下面我将分享喷雾燃烧建模的核心技术要点。
2. 喷雾燃烧模型框架解析
2.1 基础控制方程体系
喷雾燃烧仿真需要求解三大类控制方程：
- 气相Navier-Stokes方程（包含湍流模型）
- 液滴运动与传热传质方程
- 燃烧化学反应方程
以典型的RANS方法为例，时均化的连续方程和动量方程为：
\frac{\partial \bar{\rho}}{\partial t} + \nabla \cdot (\bar{\rho}\tilde{\mathbf{u}}) = S_m\frac{\partial}{\partial t}(\bar{\rho}\tilde{\mathbf{u}}) + \nabla \cdot (\bar{\rho}\tilde{\mathbf{u}}\otimes\tilde{\mathbf{u}}) = -\nabla \bar{p} + \nabla \cdot \bar{\mathbf{\tau}} + \bar{\rho}\mathbf{g} + \mathbf{S}_u其中源项Sm和Su来自液滴与气体的质量、动量交换。
2.2 液滴追踪方法对比
主流液滴追踪方法有三种：
- 拉格朗日粒子法：逐个计算液滴运动轨迹
- 优点：物理清晰，可跟踪单个液滴
- 缺点：计算量大，需要统计平均
- 欧拉多相流法：将液滴视为连续相
- 优点：计算效率高
- 缺点：难以处理大尺寸分布
- 混合方法：大液滴用拉格朗日，小液滴用欧拉
经验提示：对于工业燃烧器，推荐使用拉格朗日法结合随机轨道模型，在保证精度的同时控制计算量。
3. 关键子模型实现细节
3.1 液滴破碎模型选型
常用破碎模型对比：
| 模型类型 | 适用场景 | 计算复杂度 | 典型参数设置 |
|---|---|---|---|
| TAB模型 | 低韦伯数(<12) | 低 | C3=1.732, C4=0.5 |
| KH-RT模型 | 航空发动机 | 中 | B0=0.61, B1=10 |
| 动态拖曳模型 | 高湍流强度环境 | 高 | Cd系数动态调整 |
在柴油机仿真中，KH-RT模型表现最好。其核心公式：
\tau_{KH} = \frac{3.788B_1 r}{Ω\sqrt{We}}\tau_{RT} = \frac{C_{RT}}{√(g|a_t|/ρ_l)}其中We为韦伯数，Ω为液滴振荡频率。
3.2 蒸发模型参数设置
采用Spalding蒸发模型时，关键参数包括：
- 蒸发常数：0.6-1.2（轻质燃油取低值）
- 努塞尔数修正系数：1.5-2.0（考虑强制对流）
- 斯蒂芬流修正：必须开启
实测数据表明，对于柴油喷雾，将蒸发常数设为0.85时，液滴寿命预测误差可控制在15%以内。
4. 燃烧耦合策略
4.1 有限速率/涡耗散模型
推荐使用EDC模型处理湍流燃烧：
EDC { modelVersion 2; volumeFractionCoeffs { Cgamma 2.1377; Ctau 0.4082; } }注意调整Cgamma和Ctau参数：
- 增大Cgamma：增强湍流混合效应
- 减小Ctau：加快化学反应速率
4.2 辐射换热处理
必须考虑辐射对液滴蒸发的影响：
- 采用DO辐射模型
- 设置液滴辐射吸收系数
- 灰体辐射假设（ε=0.9）
典型设置：
radiationModel discreteOrdinate; radiation { solverFreq 1; absorptionEmissionModel greyMean; scatterModel none; sootModel none; }5. 典型问题排查指南
5.1 液滴堆积问题
现象：液滴在燃烧室壁面堆积不蒸发解决方案：
- 检查壁面温度设置（应>燃料沸点）
- 验证壁面液滴边界条件（建议设为reflect）
- 增加壁面热流耦合计算
5.2 燃烧不充分问题
诊断步骤：
- 检查当量比分布（φ=0.8-1.2为佳）
- 验证湍流混合时间尺度（应<化学反应时间）
- 监测OH自由基分布（判断火焰位置）
调整策略：
- 增大喷嘴锥角（15°→30°）
- 优化旋流器角度（45°±5°）
- 调整二次风比例（20%-30%）
6. 计算加速技巧
6.1 并行计算优化
推荐分区策略：
- 气相场：几何分解
- 液滴群：动态负载平衡
- 化学反应：CPU-GPU异构计算
实测案例：使用128核并行时，采用混合分区策略可提升效率40%。
6.2 时间步长控制
采用自适应步长策略：
Δt = min(Δt_{CFL}, 0.1τ_{evap}, 0.01τ_{chem})其中：
- τ_evap：典型液滴蒸发时间
- τ_chem：最短化学反应时间
在FLUENT中可通过UDF实现自动调整。
7. 实验验证方法
7.1 PDPA测量对比
相位多普勒粒子分析仪(PDPA)验证要点：
- 测量位置选择：
- 喷嘴下游1D、3D、5D处
- 径向间隔5mm布点
- 对比参数：
- SMD（索特平均直径）
- 轴向速度分布
- 液滴数密度
7.2 燃烧诊断技术
推荐组合方案：
- OH-PLIF：定位火焰锋面
- 高速摄影：捕捉点火过程
- 排气分析：测量CO/NOx排放
误差允许范围：
- 火焰长度：±10%
- 温度场：±8%
- 排放值：±15%
8. 工程应用案例
某型燃气轮机燃烧室改造项目：
- 原始问题：出口温度分布系数>0.25
- 仿真发现：燃油喷雾集中在头部
- 改进措施：
- 优化旋流器叶片角度（40°→35°）
- 增加预混段长度（150mm→180mm）
- 效果：
- 温度分布系数降至0.18
- NOx排放降低22%
关键仿真设置：
Physical Model: - Lagrangian Multiphase - Realizable k-epsilon - EDC Combustion - Skeletal Mechanism (29 species)9. 未来改进方向
根据最新研究趋势，建议关注：
- 机器学习辅助的破碎模型
- 用DNN预测次级液滴分布
- 可减少50%模型误差
- 高保真LES模拟
- 需配合自适应网格
- 计算量增加但精度提升
- 碳烟生成预测
- 耦合详细碳烟模型
- 需要验证光学测量数据
在最近参与的某预研项目中，采用神经网络修正的KH-RT模型使液滴贯穿距离预测误差从18%降至7%。
