SUS模型
理论与物理
概述
老师，为什么碳烟（soot）在燃烧CFD中很重要？
碳烟是不完全燃烧产生的碳微粒（粒径10-100 nm），因三个原因而重要：(1) 是废气排放法规的对象（PM：颗粒物），(2) 对辐射传热贡献大（火焰辐射主要由碳烟主导），(3) 危害健康（致癌性）。柴油发动机、航空发动机、工业炉中的碳烟预测是必不可少的课题。
请告诉我碳烟的形成机制。
碳烟的形成是一个四阶段过程。
1. 成核（Nucleation）: PAH（多环芳烃）聚合形成最初的碳烟核。C2H2（乙炔）是PAH生长的主要前驱体
2. 表面生长（Surface Growth）: 通过HACA机制（H-提取-C2H2-加成）在碳烟颗粒表面沉积碳
3. 凝并（Coagulation）: 颗粒之间碰撞、合并而变大
4. 氧化（Oxidation）: 被O2或OH氧化，碳烟燃烧消失
Moss-Brooke模型
请告诉我CFD中使用的碳烟模型的支配方程。
Moss-Brooke 双变量模型求解两个输运方程：碳烟质量分数 $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 diagram）是同时降低碳烟和NOx排放策略的基础。
碳烟模型是化学反应动力学和颗粒力学的融合啊。
是的。因为需要准确描述气相PAH化学和碳烟颗粒动力学两方面，所以它是燃烧建模中难度最高的领域之一。
碳烟是纳米尺寸的碳颗粒——用方程描述直径1nm以下物体生成的挑战
碳烟的生成是直径1~100nm的碳颗粒在燃烧过程中仅数毫秒内形成、生长、凝并的过程。为了在CFD中表达此过程，需要将“成核→表面生长→凝并→氧化”这四个过程公式化。有Fenimore-Jones的双方程模型、Frenklach的详细碳烟模型等不同流派，但特别困难的是颗粒尺寸分布（PSD）的描述。粒径1nm和100nm的行为完全不同，而CFD的网格是毫米量级——也就是说，处理的是物体尺寸仅为计算网格百万分之一的跨维度问题。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：请想象一下打开水龙头的瞬间。最初水流不稳定地喷出，过一会儿就变成稳定的水流了吧？描述这个“正在变化的过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？是只看“经过足够时间流动稳定之后”——也就是将此项设为零。计算成本大幅降低，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着运往下游吧。这就是“对流”——流体运动搬运物体的效果。暖风的暖空气能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——此项包含“速度×速度”，因此是非线性的。也就是说，流速变快此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水。因为蜂蜜的粘性（$\mu$）高，不易流动。粘性越大扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针尖有力地射出对吧？为什么呢？因为活塞侧压力高，针尖侧压力低——这个压力差成为推动流体的力。大坝放水也是同样原理。天气图中等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——冬天房间里开了暖气但暖空气不上升，得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意点·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C，水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判定指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值方法的细节
请告诉我碳烟模型的数值实现。
CFD中使用的碳烟模型大致分为三类。
| 模型 | 精度 | 计算成本 | 特征 |
|---|---|---|---|
| 经验双变量模型 | 低-中 | 低 | Moss-Brooke，输运 $Y_s$ 和 $N$ |
| 矩方法 (MoM) | 中-高 | 中 | 输运粒径分布的矩 |
| 分区法 | 高 | 高 | 通过离散分区解析粒径分布 |
经验双变量模型
请从最简单的模型开始讲。
Moss-Brooke（Fluent标准搭载）用 $Y_s$（碳烟质量分数）和 $N$（数密度）两个变量描述颗粒群。对成核、表面生长、凝并、氧化各过程使用Arrhenius型速率公式。简便但粒径分布信息只有平均值。
矩方法 (MOMIC)
矩方法是什么？
是求解粒径分布函数 $n(v,t)$（$v$ 是颗粒体积）的矩 $M_r = \int_0^\infty v^r n(v) dv$ 的输运方程的方法。$M_0$ 对应数密度，$M_1$ 对应体积分数。Frenklach & Harris 提出的MOMIC（带插值闭合的矩方法）在 Fluent 2020 及以后版本中可用。
分区法
分区法的优点是什么？
将粒径范围分割成离散的分区（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）的反应机理。全局单步机理不包含C2H2，因此无法进行碳烟计算。请使用DRM-19以上的机理。
碳烟与辐射的耦合
碳烟和辐射是如何耦合的？
碳烟颗粒发射和吸收连续光谱的辐射。碳烟的吸收系数由下式近似。
其中 $C_0$ 和 $C_2$ 是光学常数。气体辐射（CO2, H2O的带状辐射）加上碳烟的连续辐射，因此碳烟多的火焰中辐射损失大幅增加。
碳烟模型是燃烧模型＋颗粒模型＋辐射模型的三重耦合啊。
是的。模型复杂度和计算成本的平衡很重要，先用Moss-Brooke双变量模型把握趋势，必要时再转向MOMIC或分区法是实务性的做法。
碳烟模型的数值方法——输运方程法与分区法的选择
CFD中的碳烟分析主要有“双方程模型（Moss-Brookes：碳烟数密度N＋碳烟体积分数f）”和“分区法（用离散区间追踪粒径分布）”。双方程模型计算成本低，是Fluent的标准实现，但无法获得粒径分布的详细信息。分区法可以按粒径解析颗粒的成核、凝并、表面生长，预测光散射、吸收特性的精度高。最近MoM（矩方法）作为两者的折中方案受到关注，在计算成本和建模精度方面平衡性优异。航空发动机NVPM（非挥发性颗粒物）法规对应正在加速分区法的实用化。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM: 自然满足守恒定律。CFD的主流。FEM: 对复杂形状、多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式法: CFL ≤ 1 是稳定条件。隐式法: CFL > 1 也稳定，但影响精度和迭代次数。LES: 推荐 CFL ≈ 1。物理意义：一个时间步内信息传播不超过一个网格。
残差监控
相关主题
なった
詳しく
報告