烟烟模型
烟烟的理论基础
概述
老师，烟烟（soot）为什么在燃烧CFD中重要？
烟烟是由不完全燃烧产生的碳微粒子（粒径10-100 nm），具有3个重要原因：(1) 排放气体规制对象（PM：颗粒物），(2) 对辐射传热的贡献大（火焰辐射由烟烟主导），(3) 健康危害（致癌性）。柴油发动机、航空发动机、工业炉中的烟烟预测是必须的课题。
请教烟烟的生成机制。
烟烟的生成通过4个阶段过程进行。
1. 核形成（Nucleation）：通过PAH（多环芳香族烃）的聚合形成最初的烟烟核。C2H2（乙炔）是PAH成长的主要前驱体
2. 表面成长（Surface Growth）：通过HACA机制（H-Abstraction-C2H2-Addition）使烟烟粒子表面堆积碳
3. 凝聚（Coagulation）：粒子相互碰撞、合并而变大
4. 氧化（Oxidation）：通过O2和OH使烟烟燃烧消失
Moss-Brooke模型
请教CFD中使用的烟烟模型的控制方程。
Moss-Brooke 2变量模型求解烟烟质量分率 $Y_s$ 和烟烟数密度 $N$（粒子数/kg）的两个输运方程。
烟烟体积分率 $f_v$ 由下式求得。
其中 $\rho_s \approx 1800$ kg/m$^3$ 是烟烟密度，$d_p$ 是平均粒径。
烟烟生成的主要参数
什么条件下烟烟会大量产生？
总结烟烟生成的主要因子。
| 因子 | 烟烟增加方向 | 原因 |
|---|---|---|
| 当量比 | 过浓（$\phi > 1$） | 氧气不足导致不完全燃烧 |
| 温度 | 1500-1800 K | 核形成的最适温度域 |
| 压力 | 高压 | 碰撞频率增加 |
| 燃料结构 | 芳香族 > 直链 | 易产生PAH前驱体 |
| 停留时间 | 长 | 保证烟烟成长时间 |
原来1500-1800 K是烟烟生成的温度窗口。
对的。低于此温度时核形成速度缓慢，高于此温度时OH氧化占优势使烟烟燃尽。这个"烟烟生成窗"通过$\phi$-T 图谱可视化。柴油燃烧的$\phi$-T 图谱（Dec图）是烟烟和NOx同时低减策略的基础。
烟烟模型是化学反应速率论和粒子力学的融合。
对的。必须准确描述气相PAH化学和烟烟粒子动力学两方面，所以这是燃烧建模中难度最高的领域之一。
烟烟是纳米级碳粒子——用方程式描述直径1nm以下物体生成的挑战
烟烟生成是直径1～100nm的碳粒子在燃烧中短短数毫秒内形成、成长、凝聚的过程。要在CFD中表现这个过程，需要将"核形成→表面成长→凝聚→氧化"这四个过程数学化。Fenimore-Jones的2方程模型和Frenklach的详细烟烟模型等流派各有风格，特别困难的是粒子尺寸分布（PSD）的描述。粒径1nm和100nm的行为完全不同，而CFD网格的单位是毫米级——也就是说物体的尺寸比计算网格小100万倍，这是异次元的问题。
烟烟的数值计算方法
数值方法的详细说明
请教烟烟模型的数值实现。
CFD中使用的烟烟模型大体分为3类。
| 模型 | 精度 | 计算成本 | 特点 |
|---|---|---|---|
| 经验2变量模型 | 低-中 | 低 | Moss-Brooke、输运$Y_s$和$N$ |
| 方法的矩（MoM） | 中-高 | 中 | 输运粒径分布的矩 |
| Sectional法 | 高 | 高 | 用离散截面解析粒径分布 |
经验2变量模型
请从最简单的模型开始教。
Moss-Brooke（Fluent标准搭载）用$Y_s$（烟烟质量分率）和$N$（数密度）这2个变量描述粒子集团。核形成、表面成长、凝聚、氧化各过程使用Arrhenius型速率式。简便但粒径分布信息只有平均值。
方法的矩 (MOMIC)
矩法是什么？
求解粒径分布函数 $n(v,t)$（$v$为粒子体积）的矩 $M_r = \int_0^\infty v^r n(v) dv$ 的输运方程。$M_0$ 对应数密度，$M_1$ 对应体积分率。Frenklach & Harris提出的MOMIC（Method of Moments with Interpolative Closure）在Fluent 2020后可用。
Sectional法
Sectional法有什么优点？
将粒径范围离散分成多个截面（bin），各截面的数密度单独输运。粒径分布形状可任意表现，精度最高，但需要20-30个额外标量输运，计算成本很大。STAR-CCM+和CONVERGE可用。
Fluent中的设置
请教Fluent中烟烟模型的设置步骤。
1. Models > Species > Species Transport（前提是燃烧模型已设置）
2. Models > Soot > Moss-Brooke（简易）或MOMIC（推荐）
3. PAH前驱体化学种设置（C2H2、C6H6等）-- 必须包含在反应机制中
4. 与辐射模型的耦合 -- 将烟烟吸收系数传给辐射模型
重要提示：烟烟模型需要包含PAH前驱体（至少C2H2）的反应机制。全局1阶段机制不含C2H2，无法进行烟烟计算。要用DRM-19以上的机制。
烟烟与辐射的耦合
烟烟和辐射怎样耦合？
烟烟粒子放射和吸收连续谱辐射。烟烟的吸收系数近似为下式。
其中 $C_0$ 和 $C_2$ 是光学常数。气体辐射（CO2、H2O的带状辐射）加上烟烟的连续辐射，所以烟烟多的火焰辐射损失会大幅增加。
烟烟模型是燃烧模型＋粒子模型＋辐射模型的三重耦合。
对的。模型的复杂性和计算成本的平衡很重要。先用Moss-Brooke 2变量模型掌握趋势，必要时再进行MOMIC或Sectional法，这是实务方法。
烟烟模型的数值方法——输运方程法和sectional法的选择
CFD烟烟分析主要有"2方程模型（Moss-Brookes：烟烟数密度N＋烟烟体积分率f）"和"Sectional法（粒子尺寸分布用离散bin追踪）"。2方程模型计算成本低、Fluent标准搭载，但得不到粒径分布的详细信息。Sectional法可按粒径解析核形成、凝聚、表面成长，光散射吸收特性预测精度高。最近MoM（方法的矩）作为两者的折中而受关注，计算成本和建模精度的平衡优异。航空发动机的NVPM（非挥发性粒子物质）规制对应中Sectional法的实用化在加速。
烟烟的实务应用
实践指南
请教烟烟分析的实务步骤。
柴油发动机或燃气轮机烟烟预测的流程。
分析流程
1. 燃烧模型的确立 -- 首先准确计算温度场和化学种分布（特别是C2H2、OH）
2. 烟烟模型的加入 -- 燃烧模型稳定后打开烟烟模型
3. 与辐射模型的耦合 -- 将烟烟吸收系数传给辐射（DO + WSGGM + 烟烟吸收）
4. 后处理 -- 评估烟烟体积分率 $f_v$、粒径分布、辐射热流
柴油喷雾火焰中的烟烟分析
柴油发动机烟烟分析有什么特殊注意点？
柴油烟烟分析的关键是与喷雾模型的耦合。
| 要素 | 推荐设置 | 备注 |
|---|---|---|
| 喷雾模型 | KH-RT（Kelvin-Helmholtz / Rayleigh-Taylor） | 液滴分裂标准模型 |
| 蒸发模型 | Frossling | 多成分蒸发 |
| 燃烧模型 | SAGE（CONVERGE）或EDC（Fluent） | 需要详细化学反应 |
| 烟烟模型 | MOMIC或Sectional | 粒径分布重要 |
| 反应机制 | n-heptane骨架（44种/112反应） | 柴油代理燃料 |
n-正庚烷是柴油的代理。
对的。实际柴油燃料是数百种成分的混合物，但为了再现着火特性和烟烟生成倾向，用n-正庚烷或n-十二烷作代理燃料。ECN（Engine Combustion Network）的Spray A基准是柴油喷雾火焰的标准验证工况。
$\phi$-T 图谱的烟烟/NOx分析
请教$\phi$-T 图谱的使用方法。
将CFD结果各单元按当量比 $\phi$ 和温度 $T$ 散点分布。烟烟生成区（$\phi > 2$、$T = 1500-1800$ K）和NOx生成区（$\phi \approx 1$、$T > 2200$ K）重叠，可视化烟烟-NOx平衡。
常见失败与对策
| 症状 | 原因 | 对策 |
|---|---|---|
| 完全不产生烟烟 | 反应机制不含C2H2 | 换用DRM-19以上机制 |
| 烟烟过多 | 氧化模型弱 | 确认OH氧化路径 |
| 烟烟空间分布不符 | 喷雾附近网格分辨率不足 | 细化喷雾附近网格 |
| 辐射损失太大 | 烟烟吸收系数过估 | 确认烟烟光学常数 |
烟烟分析涉及燃烧、喷雾、辐射的全耦合，分阶段构建是铁则。
对的。冷喷雾 → 着火、燃烧 → 烟烟 → 辐射的顺序逐步积累。一次全打开往往收敛不了。
柴油发动机的"NOx-PM平衡"——烟烟模型支撑规制对应实务的原因
柴油发动机排气规制对应中，工程师直面的"NOx-PM平衡"矛盾：提高燃烧温度，氮氧化物（NOx）增加，降低温度、降低氧气浓度，反而烟烟（PM）增加。Euro 6、国内重型车规制对应需要同时优化EGR（排气再循环）流量、喷射时序、共轨压力这3个变量。装备烟烟模型的CFD（Fenimore-Jones模型或Two-Equation Soot模型），能将多变量优化的试验次数减至台架试验的1/10以下。"先用CFD缩小方向，最后仅在发动机台架确认"的开发流程现已成为汽车制造商的标准。
烟烟的软件比较
商用工具比较
| 工具 | 经验模型 | MOMIC | Sectional | PAH详细 | 最佳用途 |
|---|---|---|---|---|---|
| Ansys Fluent | Moss-Brooke | 2020R1后 | 无 | CHEMKIN联动 | 通用燃烧器 |
| STAR-CCM+ | 2变量 | 搭载 | 搭载 | DARS联动 | 燃气轮机 |
| CONVERGE | Hiroyasu | 搭载 | 搭载 | SAGE | 柴油发动机 |
| OpenFOAM | 社区 | 社区 | 研究实现 | Cantera联动 | 研究应用 |
柴油发动机烟烟分析还是要用CONVERGE？
CONVERGE的喷雾+燃烧+烟烟全耦合最成熟。AMR自动追踪喷雾、火焰、烟烟生成区域，无需手动网格调整，这是大优点。
Fluent固有功能
Fluent长期采用Moss-Brooke作为标准，但2020R1后加入了MOMIC。PAH前驱体可处理到C6H6（苯）甚至A4（芘）。烟烟-辐射耦合用DO/P-1模型自动设置。
STAR-CCM+固有功能
STAR-CCM+强项是标配Sectional法。20截面可解析粒径1-500 nm分布，还能与DPF（柴油微粒滤清器）耦合分析。
选择指南
最后怎样选？
烟烟模型选择由所求的输出（总量 vs 粒径分布）决定。
对的。仅与排放规制值比较，2变量模型足够，但需要粒子数规制（PN）或粒径分布评估，就非要用Sectional法不可。
烟烟CFD工具——Ansys Fluent vs OpenFOAM sootFoam
烟烟模型CFD工具以Ansys Fluent（Moss-Brookes/Enhanced Moss-Brookes）和OpenFOAM的sootFoam（sectional法）为代表。Fluent在工业燃烧器、燃气轮机设计中实绩丰富，详细化学反应与烟烟模型耦合完善。OpenFOAM的sootFoam可追踪粒径分布，多用于航空发动机排气粒子（NVPM）研究。Numeca提供燃烧专用工具FINE/Open，具备高精度的乱流燃烧-烟烟耦合分析。所有工具都能导入通用反应机制文件（CHEMKIN格式），燃料变更对应容易。
烟烟的先端研究
先端课题和研究动向
请教烟烟模型的最前沿研究。
烟烟建模是燃烧CFD中最活跃的研究领域。有3个方向。
PAH化学的精细化
PAH（多环芳香族烃）化学怎样进展？
烟烟核形成机制，除传统的HACA（H-Abstraction-C2H2-Addition）外，还有PAH间的范德华力结合（Physical Nucleation）和自由基-自由基结合（CHRCR机制）的新建议。
| 核形成机制 | 研究团队 | 特点 |
|---|---|---|
| HACA（化学成长） | Frenklach（伯克利） | 经典、高温有效 |
| 物理核形成 | Kraft（剑桥） | 范德华力结合 |
| CHRCR | Wang（斯坦福） | 自由基-自由基结合 |
| 共鸣稳定化 | Mebel（FIU） | 共鸣稳定化自由基路径 |
混合方法的矩
矩法怎样进展？
HMOM（混合方法的矩）由Mueller & Pitsch（2012）提出，克服MOMIC粒径分布假设的限制。能处理双峰（bimodal）粒径分布，核形成和表面成长同时进行的区域精度上升。
LES与烟烟模型的耦合
LES能预测烟烟？
LES的网格以下尺度烟烟生成速度非线性变动，滤波反应速度处理是课题。条件矩闭包（CMC）或输运概率密度函数与Sectional法的组合LES-CMC-Sectional是最高精度手法。计算成本极大，但Sandia flame D和航空发动机模型燃烧器的验证在进行。
机器学习的应用
烟烟研究是分子级化学和粒子物理学的最前沿。
对的。烟烟核形成机制60年来仍在讨论。CFD模型精度上升需要基础化学的进步。
烟烟模型的最前线——PAH（多环芳香族烃）成长的直接数值模拟
烟烟生成的核心是"PAH（多环芳香族烃）的成长、凝聚、表面反应"，但完全再现的计算成本巨大。2010年代后DNS（直接数值模拟）追踪了PAH形成的素反应，定量化了从乙炔（C₂H₂）经苯（C₆H₆）到芘（C₁₆H₁₀）的成长路径反应速度常数。现在的LES/RANS级烟烟模型（Moss-Brookes、semi-empirical、HACA）由DNS得到的参数验证。面向2030年代航空发动机设计，非化石燃料（SAF）的精密烟烟生成特性预测是CAE开发的最前线课题。
烟烟故障排除
故障排除
请教烟烟计算的故障。
1. 完全不产生烟烟
原因检查清单：
- 确认反应机制含C2H2。全局1阶段机制不含烟烟前驱体
- 确认烟烟模型正确启用
- 火焰温度是否在烟烟生成窗（1500-1800 K）内
- 是否有局部过浓区（$\phi > 1$）
完全预混合稀薄燃烧条件下完全没有烟烟是正常？
正常。$\phi < 0.8$ 均质稀薄燃烧下烟烟生成几乎为零。烟烟是过浓条件（柴油喷雾、富燃域）下的问题。
2. 烟烟比实验大10-100倍
原因与对策：
- Moss-Brooke模型的核形成常数过大。缺省值对乙烯火焰标定，对甲烷火焰易过估
- OH氧化弱。OH烟烟氧化由Fenimore模型描述，OH质量分率不正确会导致氧化不足
- 乱流-烟烟相互作用建模不足。RANS从平均场计算烟烟生成，未考虑局部温度、浓度变动
3. 烟烟与辐射的耦合问题
| 症状 | 原因 | 对策 |
|---|---|---|
| 火焰温度急速下降 | 烟烟辐射过大 | 确认烟烟光学常数，$f_v$是否过大 |
| 辐射反馈振荡 | 烟烟→辐射→温度→烟烟非线性循环 | 增加Under-Relaxation、分段耦合 |
| 壁面热流与测不符 | 烟烟辐射的壁面吸收模型不足 | 增加DO模型的角度分割 |
4. 计算速度慢（Sectional法）
Sectional法加20-30个截面会大幅增加化学种数，内存和计算时间激增。
对策：
- 减少截面数（10个截面也能掌握趋势）
- 设烟烟模型的Update Frequency（每5-10次迭代更新1次，而非每次）
- 用AMR仅在烟烟生成区提高分辨率，其他区降低（CONVERGE）
调试步骤
1. 首先只运行燃烧模型，确认温度、C2H2、OH分布
2. 加入烟烟模型，确认烟烟体积分率$f_v$的分布
3. 加辐射模型，确认烟烟-辐射耦合影响
4. 各阶段与实验数据比较
烟烟计算必须分阶段构建和实验数据对照。
对的。烟烟模型参数对燃料和工况依赖大，没有对象工况接近实验数据的验证，结果不能信任。
"烟烟预测与实测差异大"——反应机制和乱流模型的相互作用
CFD烟烟计算与光学计测（LII法、消光法）结果往往相差1个数量级。原因多是①乱流-化学反应相互作用模型的不正确，②烟烟前驱体（PAH）的不确定反应速度常数，③放射冷却对烟烟形成的影响。特别是乱流波动使烟烟局部形成速度变调的"乱流烟烟相互作用（TSI）"的精确建模仍是研究课题。实测与差异分解为"成分（核形成/凝聚/氧化）"的敏感度分析和CFD-实验联动方法正成为工业燃烧器设计精度上升的标准手法。
相关课题
有用
详情
报告