CN114894491B - 一种重建rbcc发动机燃烧室火焰温度二维分布的方法 - Google Patents
一种重建rbcc发动机燃烧室火焰温度二维分布的方法 Download PDFInfo
- Publication number
- CN114894491B CN114894491B CN202210454349.7A CN202210454349A CN114894491B CN 114894491 B CN114894491 B CN 114894491B CN 202210454349 A CN202210454349 A CN 202210454349A CN 114894491 B CN114894491 B CN 114894491B
- Authority
- CN
- China
- Prior art keywords
- flame
- iteration
- flame temperature
- absorption capacity
- combustion chamber
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 28
- 108091053398 TRIM/RBCC family Proteins 0.000 title claims abstract 9
- 102000011408 Tripartite Motif Proteins Human genes 0.000 title claims abstract 9
- 238000010521 absorption reaction Methods 0.000 claims abstract description 32
- 230000005855 radiation Effects 0.000 claims abstract description 30
- 230000003595 spectral effect Effects 0.000 claims abstract description 17
- 238000000701 chemical imaging Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 2
- 239000000446 fuel Substances 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004861 thermometry Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000003350 kerosene Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0014—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
- G01J5/0018—Flames, plasma or welding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/02—Details or accessories of testing apparatus
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明实施例提供了一种基于多光谱辐射强度重建RBCC发动机燃烧室火焰温度二维分布的方法。该方法首先使用多光谱成像设备测量RBCC发动机燃烧室火焰二维多光谱辐射强度,采用双色法确定火焰温度与火焰吸收能力的迭代初值,基于牛顿迭代法与最小二乘算法迭代求解火焰温度与火焰吸收能力的迭代步长,利用迭代步长修正火焰温度与火焰吸收能力直至火焰多光谱辐射强度相对误差小于允许误差,此时火焰温度迭代结果视为火焰真实温度。本发明提供的火焰温度迭代求解算法可根据火焰光谱辐射强度图像获得RBCC发动机燃烧室火焰温度二维空间分布,迭代速度快,迭代结果精度高且不受先验条件影响。
Description
技术领域
本发明涉及基于火焰光谱辐射强度重建火焰温度二维空间分布的测量方法,尤其是涉及火箭基组合循环(Rocket Based Combined Cycle,RBCC)发动机燃烧室内释热条件极其不稳定的火焰温度分布。
背景技术
RBCC发动机在不同飞行马赫数下燃烧室内存在复杂的亚、跨及超声速流动共存的流场结构特征,同时火箭支板、燃料支板等结构的存在形成强剪切流动,发动机内火焰燃烧过程表现出多重时间尺度和空间尺度耦合的显著特征。使得发动机在宽来流马赫数条件下的点火、火焰稳定与高效燃烧变得十分困难,因此对RBCC发动机燃烧室内火焰进行高精度的燃烧检测是指导RBCC发动机燃烧室燃料燃烧组织模式修正与维持火焰稳定的重要环节。然而已有辐射测温法在重建强瞬态条件下的火焰温度场时精度较低且受先验条件影响较大,测量设备也无法实现时空同步成像,因此发展一种能够基于多光谱辐射强度重建RBCC发动机燃烧室火焰温度场的方法具有重要意义。
发明内容
为了实现对时空耦合特征显著,点火、火焰稳定与高效燃烧十分困难的RBCC发动机燃烧室火焰进行高精度燃烧检测,本发明提供了一种重建RBCC燃烧室火焰温度二维分布的方法,该方法能够基于多光谱二维辐射强度图像给出时空同步的RBCC燃烧室火焰温度二维分布。
本发明的技术方案如下:
一种重建RBCC发动机燃烧室火焰温度二维分布的方法,其特征在于,包括如下步骤:
步骤1:利用多光谱成像设备测量RBCC发动机燃烧室火焰二维多光谱辐射强度图像;
步骤2:RBCC发动机燃烧室火焰光谱辐射强度表示为以下公式1:
上式中,I为火焰光谱辐射强度,单位为W/m3/sr;c1=3.742×10-16W·m2,c2=1.4388×10-2m·K;κL为火焰吸收能力,λ为波长,单位为nm;T为温度,单位为K;
步骤3:通过双色法计算火焰温度与火焰吸收能力的迭代初值T0,κL0;
步骤4:构造多光谱辐射强度偏差方程组;
步骤5:采用牛顿迭代法与最小二乘算法迭代求解κLk+1、Tk+1,k为迭代阶数,将上一次的迭代结果作为迭代初值求解火焰温度与火焰吸收能力,直至迭代结果不随迭代阶数的变化而变化,视为迭代收敛,此时的火焰温度视为火焰真实温度。
作为上述步骤4的一种优选,建立多光谱辐射强度偏差方程,所述偏差方程为以下公式2:
上式中,I为测量光谱辐射强度,单位为W/m3/sr。
作为上述步骤5的优选,所述采用牛顿迭代法与最小二乘算法迭代求解κLk+1、Tk+1,对火焰温度与火焰吸收能力进行迭代,包括:
a.根据所述公式1和所述公式2构造多光谱辐射强度偏差方程组;
b.构造所述偏差方程对火焰温度与火焰吸收能力的偏导数,求解迭代步长;
c.得到所述迭代步长后,可基于以下公式3对火焰吸收能力与火焰温度进行修正:
[kLk+1,Tk+1]=[kLk,Tk]-[ΔKLk,ΔTk]
上式中,k为迭代阶数,Δ为修正值;
d.依据修正后的火焰温度与火焰吸收能力的变化趋势判断迭代过程是否收敛:若修正后的火焰温度与火焰吸收能力不随迭代阶数变化则停止迭代,获得迭代收敛值,否则则重新执行所述步骤a将修正后的火焰温度与火焰吸收能力作为输入值重新进行迭代;
e.输出迭代收敛值。
本发明提供的计算方法能够通过多光谱成像设备测量所得多光谱辐射强度图像同时重建RBCC发动机燃烧室火焰温度二维空间分布,不依赖迭代初值的精度与任何先验条件,迭代速度快,迭代结果精度高。
附图说明
图1为本发明实施例提供的重建RBCC发动机燃烧室火焰温度二维分布的方法的总体流程图;
图2为本发明实施例提供的RBCC燃烧室火焰多光谱辐射强度二维分布;
图3为本发明实施例提供的不同迭代阶数的重建温度与多光谱辐射强度残差平方和图;
图4为本发明不同阶数由求解结果重新计算出的光谱辐射强度图;
图5为本发明实施例提供的基于牛顿迭代与最小二乘算法重建的RBCC发动机燃烧室火焰温度二维分布图。
具体实施方式
下面结合附图对本发明的具体实施方法进行说明。
本发明所述方法流程为利用多光谱成像设备测量RBCC发动机燃烧室火焰多光谱辐射强度图像,将二维多光谱辐射强度数据作为输入数据并基于牛顿迭代法与最小二乘算法修正火焰温度与火焰吸收能力,不断增加迭代阶数直至迭代收敛,最终收敛结果即为火焰真实温度,总体流程参见图1所示。
实施例
本发明重建RBCC发动机燃烧室火焰温度二维分布的方法,具体实施步骤为:
1.利用多光谱成像设备测量RBCC发动机燃烧室火焰多光谱辐射强度图像,输出多光谱成像设备测量波段范围内16张空间分辨率为512(H)×272(V)的光谱辐射强度图像,如图2所示。
2.火焰光谱辐射强度可表示为:
上式中,I为火焰光谱辐射强度,单位为W/m3/sr;c1=3.742×10-16W·m2,c2=1.4388×10-2m·K;κL为火焰吸收能力,λ为波长,单位为nm;T为温度,单位为K。
3.采用双色法为火焰吸收能力与火焰温度设置任意值为迭代初值,具体公式为:
为了获取最接近测量光谱辐射强度曲线的火焰温度,将测量光谱辐射强度与式(1)相结合构建多光谱辐射强度偏差方程组,通过求解该方程组的解获得火焰温度与火焰吸收能力,多光谱辐射强度偏差方程组如下:
式中,fj为第j波段的偏差方程,k为迭代阶数。
4.采用牛顿迭代法求解多光谱辐射强度偏差方程组,火焰温度与火焰吸收能力的迭代步长与偏差方程有如下关系:
上式中,Δ为迭代步长。
通过在求导点添加微小偏差的方式求解偏差方程对火焰温度与火焰吸收能力的偏导数:
上式中,δ(T)与δ(κL)分别为偏差方程对火焰温度与火焰吸收能力的微小偏差。
根据式(5)求得火焰温度与火焰吸收能力的修正值后,基于下式与迭代步长对火焰温度与火焰吸收能力进行修正:
[κLk+1,Tk+1]=[κLk,Tk]-[ΔκLk,ΔTk](6)
上式中,k为迭代阶数,Δ为迭代步长。
基于迭代后的火焰温度与火焰吸收能力计算RBCC发动机燃烧室火焰多光谱辐射强度曲线并以相对残差的1范数评估该迭代后辐射强度与测量辐射强度的偏差:
上式中,Imeasure为多光谱成像设备测量值,Ical为计算的辐射强度值。
若相对残差大于允许误差则将迭代后的火焰温度与火焰吸收能力作为输入值重新进行迭代,反之,则将迭代后的火焰温度作为迭代收敛值进行输出。
以图2所示RBCC发动机燃烧室内火焰的多光谱辐射强度图像关联数据为计算实例,火焰于RBCC发动机燃烧室内经由等离子体点燃煤油生成,燃料出口流量2.52kg/s,出口总温1642K。
选取燃烧室宽度15mm,高度8mm处的光谱辐射强度谱线进行迭代分析,分析结果分别见图3和图4。图3给出了基于本发明中迭代算法重建的火焰温度与多光谱辐射强度残差平方和随迭代阶数的变化。图4为基于本发明的迭代算法重建的不同阶数下火焰光谱辐射强度。图5给出了基于本发明提出的辐射测温法重建的火焰温度二维空间分布。以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为落入本发明的保护范围。
Claims (3)
1.一种重建RBCC发动机燃烧室火焰温度二维分布的方法,其特征在于,包括如下步骤:
步骤1:利用多光谱成像设备测量RBCC发动机燃烧室火焰二维多光谱辐射强度图像;
步骤2:RBCC发动机燃烧室火焰光谱辐射强度表示为以下公式1:
上式中,I为火焰光谱辐射强度,单位为W/m3/sr;c1=3.742×10-16W·m2,c2=1.4388×10-2m·K;κL为火焰吸收能力,λ为波长,单位为nm;T为温度,单位为K;
步骤3:通过双色法计算火焰温度与火焰吸收能力的迭代初值T0,κL0;
步骤4:构造多光谱辐射强度偏差方程组;
步骤5:采用牛顿迭代法与最小二乘算法迭代求解κLk+1、Tk+1,k为迭代阶数,将上一次的迭代结果作为迭代初值求解火焰温度与火焰吸收能力,直至迭代结果不随迭代阶数的变化而变化,视为迭代收敛,此时的火焰温度视为火焰真实温度。
3.根据权利要求1所述的重建RBCC发动机燃烧室火焰温度二维分布的方法,其特征在于,所述采用牛顿迭代法与最小二乘算法迭代求解κLk+1、Tk+1包括:
a.根据所述公式1和所述公式2构造多光谱辐射强度偏差方程组;
b.构造所述偏差方程对火焰温度与火焰吸收能力的偏导数,求解迭代步长;
c.得到所述迭代步长后,基于以下公式3对火焰吸收能力与火焰温度进行修正:
[κLk+1,Tk+1]=[κLk,Tk]-[ΔkLk,ΔTk]
上式中,k为迭代阶数,Δ为修正值;
d.依据修正后的火焰温度与火焰吸收能力的变化趋势判断迭代过程是否收敛:若修正后的火焰温度与火焰吸收能力不随迭代阶数变化则停止迭代,获得迭代收敛值,否则重新执行所述步骤a将修正后的火焰温度与火焰吸收能力作为输入值重新进行迭代;
e.输出迭代收敛值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210454349.7A CN114894491B (zh) | 2022-04-27 | 2022-04-27 | 一种重建rbcc发动机燃烧室火焰温度二维分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210454349.7A CN114894491B (zh) | 2022-04-27 | 2022-04-27 | 一种重建rbcc发动机燃烧室火焰温度二维分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114894491A CN114894491A (zh) | 2022-08-12 |
CN114894491B true CN114894491B (zh) | 2023-03-10 |
Family
ID=82720539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210454349.7A Active CN114894491B (zh) | 2022-04-27 | 2022-04-27 | 一种重建rbcc发动机燃烧室火焰温度二维分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114894491B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6271522B1 (en) * | 1998-05-16 | 2001-08-07 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Process for the quantitative analysis of gas volumes, specifically exhaust and waste gases from combustion systems or incineration plants, as well as systems for performing these processes |
CN104268938A (zh) * | 2014-09-18 | 2015-01-07 | 内蒙航天动力机械测试所 | 一种三维温度场的重建方法 |
CN104677495A (zh) * | 2015-03-02 | 2015-06-03 | 清华大学 | 一种基于光谱辐射强度测量火焰温度和发射率分布的方法 |
CN108225577B (zh) * | 2017-12-19 | 2020-02-18 | 东南大学 | 基于双光谱辐射信息的火焰三维温度分布重建方法 |
CN111795746A (zh) * | 2020-06-09 | 2020-10-20 | 哈尔滨工业大学 | 基于主被动光学层析融合探测的火焰多参数场协同测量方法 |
-
2022
- 2022-04-27 CN CN202210454349.7A patent/CN114894491B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6271522B1 (en) * | 1998-05-16 | 2001-08-07 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Process for the quantitative analysis of gas volumes, specifically exhaust and waste gases from combustion systems or incineration plants, as well as systems for performing these processes |
CN104268938A (zh) * | 2014-09-18 | 2015-01-07 | 内蒙航天动力机械测试所 | 一种三维温度场的重建方法 |
CN104677495A (zh) * | 2015-03-02 | 2015-06-03 | 清华大学 | 一种基于光谱辐射强度测量火焰温度和发射率分布的方法 |
CN108225577B (zh) * | 2017-12-19 | 2020-02-18 | 东南大学 | 基于双光谱辐射信息的火焰三维温度分布重建方法 |
CN111795746A (zh) * | 2020-06-09 | 2020-10-20 | 哈尔滨工业大学 | 基于主被动光学层析融合探测的火焰多参数场协同测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114894491A (zh) | 2022-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ren et al. | Machine learning applied to retrieval of temperature and concentration distributions from infrared emission measurements | |
Liu et al. | Measurement of distributions of temperature and wavelength-dependent emissivity of a laminar diffusion flame using hyper-spectral imaging technique | |
Fu et al. | Temperature measurements of diesel fuel combustion with multicolor pyrometry | |
Xia et al. | Two-step tomographic reconstructions of temperature and species concentration in a flame based on laser absorption measurements with a rotation platform | |
US11359967B2 (en) | Method for measuring actual temperature of flame by using all information of radiation spectrum and measurement system thereof | |
CN104677495B (zh) | 一种基于光谱辐射强度测量火焰温度和发射率分布的方法 | |
O'Byrne et al. | Dual-pump coherent anti-Stokes Raman scattering measurements in a supersonic combustor | |
Kranendonk et al. | Robust method for calculating temperature, pressure, and absorber mole fraction from broadband spectra | |
Wang et al. | Constraint optimization algorithm for spectral emissivity calculation in multispectral thermometry | |
CN104864977B (zh) | 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 | |
CN111751008B (zh) | 一种基于彩色火焰图像处理的锅炉炉内三维温度场分布检测方法 | |
Wang et al. | Baseline reduction algorithm for direct absorption spectroscopy with interference features | |
CN115014776B (zh) | 一种基于多光谱辐射强度图像计算超燃冲压发动机燃烧室火焰温度与发射率的方法 | |
CN115014775A (zh) | 一种测量超燃冲压发动机燃烧室温度和发射率的方法 | |
KR102637349B1 (ko) | 온도 측정 장치의 교정 방법, 온도 측정 장치의 교정 장치, 물리량 측정 장치의 교정 방법 및, 물리량 측정 장치의 교정 장치 | |
CN114894491B (zh) | 一种重建rbcc发动机燃烧室火焰温度二维分布的方法 | |
CN106896079A (zh) | 耐高温镍基合金材料的光谱发射率建模方法和测量系统 | |
CN117725853B (zh) | 一种液体火箭喷流尾焰辐射热流预测方法及装置 | |
McEvoy et al. | Methodologies and uncertainty estimates for T− T 90 measurements over the temperature range from 430 K to 1358 K under the auspices of the EMPIR InK2 project | |
Hashiguchi et al. | Spectral analysis of H2O near 7180 cm–1 to accurately measure trace moisture in N2 gas: evaluation of line shape profiles using Akaike Information Criterion | |
CN113670191B (zh) | 一种红外高光谱干涉仪关键定标参数在轨优化方法 | |
Yu et al. | Theoretical and experimental verification of wide-spectrum thermometry based on Taylor series de-integration method | |
Rhoby et al. | Hyperspectral Imaging Diagnostics of a Scramjet Cavity-based Flameholder | |
Té et al. | High emissivity blackbody for radiometric calibration near ambient temperature | |
Wood et al. | Optimisation of a tomography sensor for imaging of temperature in a gas turbine engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |