[发明专利]一种燃烧不稳定线性增长率不确定度的量化方法

说明书

本发明提供了一种燃烧不稳定线性增长率不确定度的量化方法，包括：步骤一，燃烧系统动力学过程模型构建；步骤二，线性增长率识别；步骤三，噪声强度识别；步骤四，不确定度量化。本发明基于线性随机动力学Fokker‑Planck方程，从极限环压力振荡数据中提取噪声强度。利用蒙特卡罗模拟方法数值求解非线性随机微分方程，生成100组以上的从小扰动到有限幅值增长过程时间序列。选取不同时间区间，比较识别结果的平均值和标准差，选择偏差和方差较小的区间作为线性增长阶段。识别各组线性增长阶段的增长率并取平均值和标准差，来量化该工况下的不确定度。

技术领域

本发明属于航空航天技术领域，涉及燃烧系统，具体涉及一种燃烧不稳定线性增长率不确定度的量化方法。

背景技术

航空航天发动机、燃气轮机等燃烧系统研制过程中经常会遇到燃烧不稳定问题，通常表现为燃烧室声模态频率附近的高幅值压力振荡。燃烧不稳定会加剧燃烧室壁面热传递和结构振动，进而导致系统性能下降和组件失效等严重后果。燃烧不稳定涉及燃烧、流动和声场间复杂的耦合过程，其产生机理、主被动控制、裕度评估和预测方法一直被学术界和工业界所关注。燃烧系统中一旦出现燃烧不稳定，通常有减少激励能量和增加燃烧室声学阻尼两种解决措施。由于燃烧现象的复杂性，减少激励能量方法的难度较大。相比较而言，增加燃烧室声学阻尼的设计要容易一些。

线性增长率表征小扰动下燃烧室振荡能量的变化速率。对于出现燃烧不稳定的燃烧室，新设计阻尼装置的衰减率应大于燃烧室的增长率并保留一定的裕度，级需要燃烧室不稳定声模态的增长率作为设计输入参数。学者们提出了多种方法来识别增长率。主要分为外部激励方法(输入-输出识别)和无外部激励方法(仅输出识别)。外部激励方法中，可以利用分叉前的噪声相干共振特性，识别系统的非线性燃烧响应，外推出系统分叉点的位置和类型、增长率、极限环振荡幅值等。实际发动机试验中加入特定形式的外部激励比较困难，因此外部激励方法在工程应用中不方便。

无外部激励的方法中，考虑自发燃烧不稳定从小扰动增加到有限幅值(极限环)的过程，振荡幅值初始增长阶段(满足线性假设)以指数形式增长，因此可以对该阶段的压力振荡时间序列瞿包络，用指数函数去拟合，或将包络取对数后利用线性函数去拟合，从而识别增长率。

现有无外部激励识别方法存在以下不足：实验数据存在噪声时，如何量化识别的燃烧不稳定线性增长率的不确定度，如何确定振荡幅值从小扰动增加到有限幅值过程中的线性增长阶段，这些问题在现有技术中还没有得到很好地解决。

发明内容

针对现有技术存在的不足，本发明的目的在于，提供了一种燃烧不稳定线性增长率不确定度的量化方法，以解决上述背景技术中提出的实验数据存在噪声时，燃烧不稳定线性增长率识别中的不确定度量化问题，以及振荡幅值从小扰动增加到有限幅值过程中确定线性增长阶段的问题。

为了解决上述技术问题，本发明采用如下技术方案予以实现：

一种燃烧不稳定线性增长率不确定度的量化方法，该方法包括以下步骤：

步骤一，燃烧系统动力学过程模型构建：

燃烧系统的燃烧室的压力震荡为单阶声模态振荡，燃烧室的压力振荡用非线性振动方程描述：

式中：

表示η的二阶导数；

表示η的一阶导数；

η表示脉动压力，单位为Pa；

ω₀表示声模态角频率，单位为rad/s；

ν表示线性增长率，单位为rad/s；

κ表示非线性饱和系数，单位为1；

ξ表示燃烧噪声，单位为Pa；