燃烧不稳定仿真从机理到实操：热声振荡判定、路线选择与排错指南

干这行的人应该都有体会——“燃烧不稳定”这五个字，比“难收敛”更让人头大。模型燃烧室跑了几万步，压力信号突然开始规律性振荡，幅度越来越大，火焰像被什么东西“推着”一呼一吸，燃烧室壁面的脉动载荷直线上升。这不是数值炸了，而是你撞上了真实的物理：热声振荡。这个主题在燃烧仿真系列里常被单独拿出来讲，就是因为它的物理机制、数值方法和普通燃烧流场仿真完全是两套逻辑。本文我会从机理、判定准则、仿真路线选择、实操流程到调试排错，完整梳理一遍，特别适合做燃气轮机、航空发动机加力室、火箭推力室以及工业燃烧器仿真的工程师和研究生参考。

1. 先搞懂“燃烧不稳定”到底在闹什么

1.1 一场压力与火焰的“合谋”

燃烧不稳定的本质不是单纯的压力脉动，而是压力振荡与火焰放热振荡之间形成了正反馈。可以把燃烧室想象成一个声学封闭腔体，内部存在一系列声学模态，比如一阶纵向驻波、切向驻波。火焰在某个位置燃烧，如果它的放热速率恰好以和当地压力振动相同的节奏波动，也就是“放热脉动”与“压力脉动”同相，那么火焰就像一只不停给秋千施加推力的手，每一轮振荡都会从燃烧释放的能量里抽取一部分，喂给声场。声场变强了，压力脉动变大，反过来又影响流场和混合过程，促使火焰放热脉动进一步增强。这样循环下来，小扰动被不断放大，最终形成大幅值的自激振荡。

这里有个很关键的概念：时滞。火焰对上游扰动（比如速度波动、当量比波动）的响应不是瞬时的，从扰动产生到放热响应的延迟时间就是时滞。声波周期是固定的，如果时滞恰好处在让放热波动与压力波动相位配合的区间，正反馈条件就成立；如果相位不对，火焰反而会成为声波的阻尼器。所以你会发现，同一个燃烧室，改一改喷嘴结构，或者换一种燃料，稳定性边界可能完全不一样——本质上改的是火焰响应特性和时滞。

1.2 为什么工程上对它闻风丧胆

普通燃烧仿真关心的是温度场、速度场、污染物排放，压力脉动可以当成噪声指标来看。但燃烧不稳定一旦发生，问题往往是结构级的。我是见过真实工程案例的：某型样机在做耐久试验时，燃烧室压力脉动幅度达到平均压力的10%以上，下游涡轮导向器边缘出现了明显的疲劳裂纹，测点支架直接被振断。这种高周疲劳破坏在热端部件上特别致命，因为燃烧室的声学频率通常在数百赫兹到数千赫兹，恰好落在金属材料的高周疲劳敏感区间，几小时内就能积累大量的交变应力循环。

除了结构破坏，振荡还会带来一系列连锁反应：大幅压力脉动导致热释放率剧烈波动，局部热负荷峰值明显升高，热障涂层加速剥落；火焰瞬时回退进入喷嘴区域造成烧蚀；燃气侧的传热规律完全偏离设计值，冷却设计失效。在某些工况下，振荡甚至会造成熄火或者回火，直接中断燃烧。这就是为什么发动机研制流程中，燃烧不稳定仿真是“一票否决”级别的检查项——你可以不计较NOx差几个ppm，但你不可以无视压力振荡放大趋势。

1.3 判定方向的物理标尺：Rayleigh准则

判断热声不稳定的经典判据是瑞利准则（Rayleigh Criterion）。它的本质其实很朴素：在一个声学周期内，如果放热脉动对应的体积分与压力脉动的乘积的时均值大于零，则局部区域在向声场注入能量，也就是“激励”；反之则是在“吸收能量”，也就是“阻尼”。用数学表达就是瑞利积分：

R = ∭∫ p'(x,t) · q'''_'(x,t) dt dV

对这个积分更工程化的理解，是看每个空间位置的“瑞利指数”——压力脉动和放热脉动的相位差。相位差在±90°以内时，燃烧对声场贡献正功，振荡被放大；相位差超过90°时，燃烧消耗声能，振荡被抑制。所以仿真分析燃烧不稳定性的时候，最核心的一步不是看压力云图有多漂亮，而是提取压力脉动与放热脉动的相位关系，计算瑞利积分的空间分布。

但要注意，瑞利准则只是必要条件而非充分条件。一个系统究竟稳不稳定，要看激励和阻尼的综合平衡。阻尼包括声波向上下游出口的辐射、壁面的声吸收、液滴蒸发和湍流耗散、喷嘴的声阻抗等。即使瑞利积分为正，如果阻尼比激励更强，系统依然是稳定的。所以定量评估时，需要把线性增长率算出来：增长率 = 激励贡献 - 阻尼贡献。正增长率对应不稳定，负增长率对应稳定。

2. 仿真路线怎么选：三条技术路径的取舍

2.1 低阶声学网络模型：设计初期的快速预筛

燃烧不稳定性仿真和普通流场仿真在方法论上最大的区别，是“你不可能一上来就做高保真三维非定常计算”——结构还没定、参数还在扫，每一轮都跑LES是不现实的。所以行业里成熟的流程是多级递进：先低阶粗判，再中阶提取火焰响应，最后高阶精确验证。

最低阶的方法就是声学网络模型，也叫集总参数声学模型。把燃烧室拆成声学元件：直管段、变截面、突扩腔、喷嘴、容腔，每个元件用传递矩阵或者阻抗边界来描述，整个系统串成一个电路一样的网络；火焰通过一个火焰传递函数元件接进去，描述放热脉动对上游速度脉动的响应。然后求解这个网络的特征方程，得到各阶声学模态的复频率。复频率的实部是振荡频率，虚部的正负就代表线性稳定还是不稳定，虚部越大说明失稳趋势越强。

这个路线的优点非常突出：计算量几乎可以忽略，一张设计参数表扫几千个工况也就几秒钟；而且物理图景清晰，你可以直接看到是哪个模态最容易激发、哪些几何参数对阻尼影响最大。缺点也明显：火焰传递函数本身通常来自经验公式或半经验模型，几何做过多简化之后，绝对精度有限。所以它的正确用途是“趋势筛选”和“机制解释”，不能作为最终交付的稳定性结论。

2.2 火焰描述函数（FDF/FTF）：连接火焰与声场的桥梁

如果你问做燃烧不稳定性仿真的人，项目中最值钱的数据是什么，十有八九会回答：火焰响应函数。火焰传递函数FTF（Flame Transfer Function）定义很简单，就是放热脉动与上游速度脉动在某个频率下的复数比值；如果把增益随扰动幅度的变化也考虑进去，就升级成火焰描述函数FDF（Flame Describing Function）。

为什么要花大力气去测这个函数？因为声学网络模型里最不确定的输入、LES非线性分析里最关键的机理信息都在这里。火焰响应包含了时滞、饱和效应、剪切层动力学、当量比扰动耦合等一系列物理信息。在工程操作中，这个函数可以通过高保真LES加受迫扰动来提取：在入口叠加单频正弦速度扰动，记录下游火焰面的体积放热脉动，做傅里叶变换得到该频率下的增益和相位。扫完一系列频率，就得到完整的FTF曲线。相位曲线里隐藏着时滞信息，增益曲线里能看到哪些频段最容易放大扰动。

实验测量也做类似的事，用扬声器阵激励空气流动，用OH·化学发光信号代表放热脉动。但实验成本高、测点受限，仿真提取FTF已经成为主流做法。需要注意，提取FTF时的扰动幅值一定要控制在“线性区”，否则增益曲线会混入非线性饱和效应，拿去给线性稳定分析用就会失真。

2.3 高保真LES仿真：终极验证

到了方案冻结、需要给设计报告“盖章”的阶段，就得靠三维高保真LES来兜底。和前面两条路线的最大区别是，LES不做线性化假设，也不预设火焰响应函数，它直接求解非稳态湍流流动、混合和化学反应过程，理论上可以自己“涌现”出振荡。

LES能捕捉到极限环振荡——线性分析只能预测“会不会失稳”，但无法预测“失稳之后振荡幅度多大”，因为非线性饱和会让振幅稳定在一个平台。极限环振幅对工程非常有价值：它决定了脉动载荷的大小，直接影响结构强度评估。此外LES还能给出振荡的空间结构、火焰在振荡过程中与壁面的相互作用、局部热流脉动等细节，这些是低阶模型完全给不了的。

当然代价是巨大的：网格量动辄几千万，时间步长受声学和化学时间尺度限制，一个工况计算量可能是一般燃烧流场仿真的10倍以上。而且数值格式必须低耗散，否则声波在传播过程中被数值阻尼吃掉，得到“假稳定”的结果，比算不明白更危险。还有一点很容易被忽略：LES因为截断了小尺度湍流，对阻尼的模拟能力存在先天缺陷，尤其是壁面边界层、小尺度耗散这类机制，所以哪怕高保真计算显示稳定，工程上也会保守处理，给足余量。

我把三条路线的定位总结在下面这张表里：

3. 一个典型的燃烧不稳定性仿真项目实操

3.1 几何建模与网格设计的两个原则

这套流程在实际项目中怎么落地？我以一个简化的某型燃气轮机模型燃烧室作为模拟项目X来展开。几何上不追求完整复刻真实产品，保留关键声学特征即可：上游供气段、喷嘴出口、突扩台阶、主燃烧区、收缩出口。声学仿真对几何的一个特殊要求是“变截面”不能被随便抹平，因为截面突变会产生声波反射，直接影响模态频率和振型；同样长度的燃烧室，有没有突扩段，一阶纵向模态频率可能差出百分之十几。

网格设计要遵守两个原则。第一，火焰区分辨率必须细化到能解析大尺度湍流结构和火焰面褶皱。对中等压力、中等温升的实验室尺度燃烧室，近喷嘴区网格尺寸在0.2mm到0.5mm量级；主燃区适当放稀到1mm左右。总网格量随燃烧室尺寸不同差异很大，实验室尺度一般在2000万到4000万之间。第二，声波传播方向上的网格过渡要平缓，避免突然的网格疏密跳变引起声波数值反射。这个道理跟光学里的阻抗失配类似，网格密度变化相当于声阻抗突跃，部分声波会在界面被“反射”回来，造成假模态。

3.2 边界条件、燃烧模型与求解参数

边界条件在热声仿真里比普通流场仿真讲究得多。入口我给的是质量流量和总温，并在后面加入受迫扰动的入口脉动段，这样不会因为直接在边界上叠加脉动而污染入口流场；出口必须使用非反射或者部分反射边界。标准压力出口会把压力波“原样弹回”计算域，形成虚假的高反射腔。我是用特征边界条件配合出口阻抗模型处理：先估计出口导向器对声波的反射系数，把反射系数折进边界条件里。

燃烧模型的选取是热声仿真里最容易被低估的一环。简单涡耗散模型虽然在小火焰模拟中很流行，但它在构造放热脉动相位时过于粗放，反应速率由混合时间尺度直接控制，放热响应与真实化学动力学脱节，拿去做稳定性分析可信度很低。稳定的做法是用火焰面模型或者有限速率化学机制。如果目标燃尽率和NOx预测精度要求不高，但放热脉动相位必须可靠，可以选简化三步或四步总包机理，配合火焰面密度输运模型，兼顾计算效率和相位保真度。

时间步长的设定有个实用公式：先按流动CFL条件估算，Δt流动 ≈ 0.5Δx / U_max，如果局部网格0.5mm，最高速度30m/s，Δt就是微秒量级。然后再检查声学CFL：声速约1000m/s，声学CFL ≈ cΔt/Δx = 2.5，看起来很大，但压力基求解器配合全隐式时间推进，声学波的解析能力更多取决于每周期内的步数，一般建议每个声学周期至少50步以上。500Hz对应的周期是2ms，按2微秒步长算，一个周期约1000步，分辨率充足。

3.3 在入口叠加扰动：扫频的具体操作方法

提取火焰响应函数或者稳定性信息，核心操作是“受迫扰动扫描”。我在模拟项目X里是这样做的：先跑一个纯冷态或者稳态燃烧LES，等待流场充分统计收敛，然后开始在每个入口通道叠加速度脉动：

u(t) = U0 · (1 + A·sin(2πft))

扰动幅值A取3%到5%之间。为什么取这个值？太小了信噪比不够，扰动信号淹没在湍流脉动里；太大了又会触发非线性饱和，提取出来的增益和相位不代表线性响应。3%到5%正好能保证扰动信号在FFT中清晰地凸显出来，又不至于偏离线性区。频率扫描从150Hz开始，以25Hz为间隔扫到900Hz，覆盖该尺度燃烧室最可能激发的前几阶纵向模态。值得一提的经验是，频率一定要扫描过去，不要只算预测模态频率附近——火焰响应的峰值频段和声学模态频率不一定重合，往往“失稳模态”恰恰落在响应增益较高而声学模态阻尼较弱的某个组合频段上，只算预测频点等于碰运气。

每个频点要跑多久？这里有个质量指标：至少覆盖5个声学周期以上再做统计。但只跑5个周期是不够的，统计收敛需要更多样本。我的习惯是每个频点跑50ms物理时间，约合250个声学周期，压力探针和放热统计曲线确认进入周期性稳态后，再额外记录10到20ms的时序数据用于后处理。这个时间尺度下，一个频点的LES大概需要数万步迭代，整条扫频曲线的计算量确实不小，但这是拿高保真数据的必要代价。

3.4 数据处理：从原始信号到稳定判据

计算结束后，真正的分析工作才开始。我通常在燃烧室内沿轴向布置7个压力探针，火焰区放一组积分探头用来统计放热脉动。数据处理有几个要点直接决定结论准不准。

第一，瞬态段必须剔除。入口开始强迫扰动后的前几个周期，流场处于暂态调整过程，这段信号不能进入FFT统计，否则会引入非物理的频谱成分。我用的是每段信号先去均值、去线性趋势，然后乘以Hann窗，重叠50%分段平均。第二，热释放信号不能简单用某几个网格点的温度脉动代替。对预混火焰，放热率脉动要通过对整个火焰区的化学反应放热项做体积积分得到；对非预混火焰，还应考虑当量比脉动带来的影响，必要时用简化火焰面参数来重构放热信号。第三，相位和增益的计算要在频域做。对每个激励频率，取该频率处的压力脉动复振幅p̃和放热脉动复振幅Q̃，计算火焰传递函数F(f) = (|Q̃|/Q̄)/(|ũ|/Ū)，相位就是arg(Q̃/ũ)。

瑞利积分最终可以化成一个频域形式的判据：按频率把压力脉动和放热脉动的相位差算出来，如果相位差落在±90°之内，这个频率就是“激励频率”候选；接着把激励位置对应的声学模态阻尼估算出来，如果激励能量大于阻尼，系统线化预测是不稳定的。实际操作里，我会做一张瑞利指数空间分布云图，看看是火焰根部还是火焰中段在“推”声场，这对指导改结构非常有帮助——挪一挪喷嘴位置或者改变掺混孔，本质上就是在改变瑞利指数的空间分布。

4. 常见问题与排查技巧实录

4.1 频率对不上？先查声学边界和网格

做热声仿真最常遇见的问题，就是数值预测的自激振荡频率和实验或者理论解析对不上。如果偏了超过10%，第一步不是怪火焰模型，而是检查声学边界条件。出口反射系数设置的偏差对模态频率的影响是系统性的，尤其对半波和四分之一波模态，出口阻抗直接改变声学“管长”的有效长度。排查办法是：先用一个不点火的冷态算例，只送空气，在入口加一个宽频扰动，统计出声学模态频率，和解析公式对比。冷态声模态算准了，再考虑燃烧温度场对声速的影响是否被正确捕捉。最后才检查网格分辨率——声波在粗网格区域的数值色散会造成相速度偏差，这种偏差通常是高频模态更明显。

4.2 压力振荡发散？不一定是物理发散

还有一类问题是扰动后压力振荡持续增长，甚至出现数值爆炸。这个情况要冷静区分：到底是物理不稳定正在发展，还是数值格式在“造假”。我归纳了一个快速分类法：先把燃烧关掉，做冷态受迫扫描，如果冷态算例加扰动后压力不能衰减，那一定是你边界条件或者数值格式有问题——真实的等温腔体声音是会衰减的。冷态正常，点火后压力振荡确实增长，再检查增长速率是否和线性稳定性理论的预期吻合。还有一个常见的数值元凶是时间步长太大：在非定常压力基求解器里，内迭代次数不足时，跨声速压力波每步都来不及完全传播，积累起来会造成压力“爬升”假象。遇到这种情况，先把步长缩小一半试算，再对比增长趋势是否变化，如果趋势变化大，基本就是时间分辨率不足。

4.3 FTF相位测不准？信号处理细节

提取火焰传递函数时相位抖动大，是大家都踩过的坑。这里我吃过一次大亏：入口强迫幅值取到了15%，结果火焰增益曲线在低频段出现明显的凹陷，相位曲线也出现“非物理”的跳变。后来把幅值降到4%，曲线立刻变得平滑。原因就是扰动太大以后，火焰面的非线性饱和效应让响应不再是正弦形式，高阶谐波能量变强，FFT主频幅值被“分流”了。另外，相位测量最怕频谱泄漏，尤其是当激励频率恰好不在FFT分辨率整数倍附近的时候，因此信号长度要保证覆盖整数个扰动周期，或者用Hann窗配合足够的章数平均。如果做临时修正，也可以把强迫信号本身拿来做锁相参考，对所有统计量按强迫相位进行同步平均，能显著压低湍流噪声。

4.4 避坑速查表

我把这几年在热声仿真上积累的经验教训整理成一张表，算是我自己的“现场手册”：

这份表里最后一条最要命。仿真显示线性增长率是正的，但实验装置稳定运行，通常意味着阻尼被低估了。真实燃烧室的声波衰减来源很复杂：表面冷却孔排布的声学吸声效应、冷却射流对声场的散射、出口导向叶片的反射特性，都很难在简化模型里准确还原。这种情况下，不要急着把责任推给实验，先用强迫振荡实验数据校验你的阻尼项，有时候一个合理的阻尼修正就能把“假不稳定”变成稳定结论。

另一条容易被忽视的经验是：压力探针的位置不要放在声学模态节点上。我见过有人把唯一的压力测点正好布在驻波压力节点附近，结果脉动幅值接近于零，差点得出“没有任何不稳定风险”的错误结论。多布置几个测点、沿轴向跨模态节点分布，是花小钱办大事的做法。

5. 一点个人感受收尾

做了这么多年燃烧不稳定性仿真，最大的体会是：不要被模拟结果本身牵着走。LES给出一个“不稳定”结论，你要追问这是不是数值格式耗散不足造成的假增长；LES给出“稳定”，你也要追问阻尼项是不是给得太充裕了，把这个系统推到了不合理的稳定位置。每次交报告之前，我都会拿低阶声学模型和LES的结果交叉对一遍，确认模态频率、增长趋势、瑞利积分的空间分布三者自洽。这个过程很费时间，但能挡掉大多数“看起来很美”的假结论。燃烧不稳定性是仿真和物理结合最紧密的领域之一，参数、边界、数值、后处理，每一个环节都在考验你对机理的理解有多深。搞定了它，你对燃烧流的认知会上一个台阶。