一种重建超燃冲压发动机燃烧室火焰温度与发射率二维分布的方法技术

一种基于多光谱辐射强度重建超燃冲压发动机燃烧室火焰温度与发射率二维分布的方法，首先使用多光谱成像设备测量超燃冲压发动机燃烧室火焰二维多光谱辐射强度，为待测火焰温度与发射率赋任意初值并设置阻尼系数，采用优化Levenberg

【技术实现步骤摘要】
一种重建超燃冲压发动机燃烧室火焰温度与发射率二维分布的方法


[0001]本专利技术涉及基于火焰光谱辐射强度重建火焰温度与发射率二维空间分布的测量方法，尤其是涉及重建超燃冲压发动机燃烧室内处于强瞬态的火焰温度与发射率分布。

技术介绍

[0002]超燃冲压发动机燃烧室内燃料流动和燃烧机理极为复杂，由于燃烧室长度有限,气流速度高,燃料在燃烧室内的滞留时间很短。在这极短的时间内要完成燃料的喷射、雾化、与空气掺混、点火和燃烧,难度很大。对燃烧室火焰进行高精度的燃烧检测以保证火焰稳定燃烧和提高燃烧效率是解决该问题的主要方法。众多燃烧检测方法中辐射测温法以其非接触、抗干扰性强的优点广泛应用于发动机火焰燃烧检测领域。
[0003]然而已有辐射测温法中，由于已有发射率模型与迭代算法受迭代初值影响较大并且测量设备无法实现时空同步成像，因此无法对超燃冲压发动机燃烧室内火焰进行精确的燃烧检测。所以发展一种基于能够基于多光谱辐射强度图像同时重建超燃冲压发动机燃烧室火焰温度与发射率分布的方法具有重要意义。

技术实现思路

[0004]为了实现对流畅环境复杂，火焰传播与释热过程呈现强不稳定的超燃冲压发动机燃烧室火焰进行高精度燃烧检测，本专利技术提供了一种基于多光谱辐射强度图像的同时重建火焰温度与发射率二维空间分布的方法，该方法能够基于多光谱二维辐射强度图像给出时空同步的超燃冲压发动机燃烧室火焰温度与发射率二维空间分布。
[0005]本专利技术的技术方案如下：
[0006]一种基于二维多光谱辐射强度图像，同时重建超然冲压燃烧室火焰温度与发射率二维空间分布的方法，其特征在于本方法包含以下步骤：
[0007]1)利用多光谱成像设备测量超燃冲压发动机燃烧室火焰二维多光谱辐射强度图像；
[0008]2)基于HOTTEL发射率模型，超燃冲压发动机燃烧室火焰光谱辐射强度可表示为：
[0009][0010]式中，I为火焰光谱辐射强度，单位为W/m3/sr，c1＝3.742
×
10
‑
16
W
·
m2，c2＝1.4388
×
10
‑2m
·
K，κL为火焰吸收能力，λ为波长，单位为nm，T为温度，单位为K；
[0011]3)设置任意值为火焰吸收能力κL0与火焰温度T0作为迭代初值；
[0012]4)设置初始阻尼系数μ与缩放系数η；
[0013]5)根据拟海森矩阵H构造增量方程并计算迭代步长，更新阻尼系数，求解κL
k+1
，T
k+1
，k为迭代阶数；将上一阶迭代结果作为迭代初值求解火焰吸收能力与火焰温度，直至迭代结果不随迭代阶数的增加而增加，视为迭代收敛，此时的火焰吸收能力与火焰温度视为火
焰真实吸收能力与真实温度。
[0014]本专利技术所述的上述技术方案中，步骤4)中，初始阻尼系数μ设置为0.01，缩放系数η设置为10，以适用于超燃冲压发动机燃烧室火焰的温度与发射率重建。
[0015]本专利技术所述的上述技术方案中，步骤5)中采用优化Levenberg
‑
Marquarelt算法迭代求解超燃冲压发动机燃烧室火焰温度与发射率，其包含以下子步骤：
[0016]a.根据本次迭代的火焰吸收能力、火焰温度与阻尼系数输入值计算拟海森矩阵；
[0017]b.根据步骤a计算所得拟海森矩阵构造增量方程计算火焰吸收能力与火焰温度的迭代步长；
[0018]c.得到迭代步长后，可基于下式对火焰吸收能力与火焰温度进行修正：
[0019][κL
k+1
,T
k+1
]＝[κL
k
,T
k
]‑
[ΔκL
k
,ΔT
k
]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0020]式中，k为迭代阶数，Δ为修正值；
[0021]d.依据迭代前后的多光谱辐射强度残差平方和对阻尼系数进行迭代修正；
[0022]e.根据更新后的火焰吸收能力与火焰温度的变化趋势，判断迭代是否收敛，若收敛则执行步骤f，若未收敛则返回步骤a，将迭代修正后的火焰吸收能力、火焰温度与阻尼系数作为下次迭代的输入值进行上述过程；
[0023]f.输出迭代收敛值。
[0024]本专利技术提供的计算方法能够通过多光谱成像设备测量所得多光谱辐射强度图像同时重建超燃冲压发动机燃烧室火焰温度与发射率空间分布，不依赖迭代初值的精度与任何先验条件，结果可靠，精度较高。
附图说明
[0025]图1为本专利技术所述方法的总体流程图；
[0026]图2基于Levenberg
‑
Marquarelt算法重建超燃冲压发动机燃烧室火焰多光谱辐射强度二维分布；
[0027]图3为不同阶数重建温度与多光谱辐射强度残差平方和；
[0028]图4为重建光谱辐射强度与光谱发射率；
[0029]图5为基于Levenberg
‑
Marquarelt算法重建超燃冲压发动机燃烧室火焰温度二维分布；
[0030]图6为基于Levenberg
‑
Marquarelt算法重建超燃冲压发动机燃烧室火焰发射率二维分布。
具体实施方式
[0031]下面结合附图对本专利技术的具体实施方法进行说明。本专利技术的具体实施流程见图1。利用多光谱成像设备测量超然冲压发动机燃烧室火焰多光谱辐射强度图像，将二维多光谱辐射强度数据作为输入数据并基于HOTTEL发射率模型与LM优化算法计算火焰吸收能力与火焰温度的迭代步长，根据迭代步长修正火焰吸收能力与火焰温度后，依据迭代前后火焰多光谱辐射强度的残差平方和对阻尼系数进行放缩调整，将迭代修正后的火焰吸收能力、火焰温度与阻尼系数输入下一次迭代，不断重复上述步骤直至迭代收敛，即为火焰真实吸收能力与温度。
[0032]一、具体实施过程
[0033]利用多光谱成像设备测量超燃冲压发动机燃烧室火焰多光谱辐射强度图像，输出多光谱成像设备测量波段范围内16张空间分辨率为512(H)
×
272(V)的光谱辐射强度图像，如图2所示。
[0034]基于HOTTEL发射率模型，火焰光谱辐射强度可表示为：
[0035][0036]式中，c1＝3.742
×
10
‑
16
W
·
m2,c2＝1.4388
×
10
‑2m
·
K，κL为火焰吸收能力，λ为波长，单位为nm，T为温度，单位为K；
[0037]发射率与火焰吸收能力κL有如下关系：
[0038][0039]为火焰吸收能力与火焰温度设置任意值为迭代初值，设置初始阻尼系数μ＝0.01，阻尼系数的缩放系数η＝10。依据目标函数式(1)本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于多光谱辐射强度图像重建超燃冲压发动机燃烧室火焰温度与发射率二维空间分布的方法，其特征在于，所述方法包括如下步骤：1)利用多光谱成像设备测量超燃冲压发动机燃烧室火焰二维多光谱辐射强度图像；2)基于HOTTEL发射率模型，超燃冲压发动机燃烧室火焰光谱辐射强度可表示为：式中，I为火焰光谱辐射强度，单位为W/m3/sr，c1＝3.742
×
10
‑
16
W
·
m2，c2＝1.4388
×
10
‑2m
·
K，κL为火焰吸收能力，λ为波长，单位为nm，T为温度，单位为K；3)设置任意值为火焰吸收能力κL0与火焰温度T0作为迭代初值；4)设置初始阻尼系数μ与缩放系数η；5)根据拟海森矩阵H，计算迭代步长并更新阻尼系数，求解κL
k+1
，T
k+1
，k为迭代阶数；将上一阶迭代结果作为迭代初值求解火焰吸收能力与火焰温度，直至迭代结果不随迭代阶数的增加而增加，视为迭代收敛，此时的火焰吸收能力与火焰温度视为火焰真实吸收能力与真实温度。2.根据权利要求1所述的基于火焰多光谱辐射强度图像计算超燃冲压发动机燃烧室火焰温度与发射率二维空间分布的方法，其特征在于：在步骤4)中，初始阻尼系数μ设置为0.01...

【专利技术属性】
技术研发人员：孙明波，郑树，朱家健，李庆，才伟光，王岩，张明轩，陆强，
申请(专利权)人：华北电力大学，
类型：发明
国别省市：