CN104864977B - 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 - Google Patents
一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 Download PDFInfo
- Publication number
- CN104864977B CN104864977B CN201410820718.5A CN201410820718A CN104864977B CN 104864977 B CN104864977 B CN 104864977B CN 201410820718 A CN201410820718 A CN 201410820718A CN 104864977 B CN104864977 B CN 104864977B
- Authority
- CN
- China
- Prior art keywords
- wavelength
- temperature
- spectral
- spectral radiance
- rocket engine
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法,本发明基于普朗克定律,通过采集燃气辐射光谱波长与对应的光谱辐射强度数据,将燃气光谱波长和光谱强度的实测关系曲线与不同温度下的普朗克曲线进行匹配,根据匹配程度,确定火箭发动机燃气实际温度;本发明所能测试温度最高可达5000K,远远超过热电偶的测温范围,更适应于火箭发动机燃气高温测量;本发明能够采集多达数百个波长及所对应的光谱辐射强度数据,提高了波长和光谱辐射强度数据拟合的精度,具有较强的抗干扰能力。相比于双波长测温法和8波长光谱仪测温,能更准确计算燃气温度,适用性更强。
Description
技术领域
本发明涉及火箭发动机燃气温度的测试方法,属于温度测量技术领域,具体涉及一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法。
背景技术
火箭发动机燃气温度在推进剂燃烧机理研究、发动机结构设计和发动机防烧蚀设计方面具有重要的参考价值。燃气温度,一方面是火箭发动机推进剂的重要能量参数,为推进剂配方调节提供设计指导;另一方面,是表征流场特征及火箭发动机内弹道性能的主要参考依据;同时,燃气温度也是影响火箭发动机内绝热层烧蚀性能的重要参数,为发动机热防护提供设计依据。但是,由于火箭发动机工作环境恶劣,且燃气温度非常高(通常大于3000K),燃气温度的测试是一项要求高、难度大的复杂技术,多年来一直没有真正解决。目前,燃气温度测试方法基本上可以分为接触式测温法和非接触式测温法两大类。
接触式测温法如热电偶法,其局限性在于测量温度范围十分有限,一般用于2500K以下的温度测量,应用到火箭发动机燃气温度测量时,热电偶无法测量高温燃气温度;其次,热电偶测温响应慢,具有一定的滞后性,不适合火箭发动机燃气温度的动态快速响应。
非接触测温方法中,适合火箭发动机燃气温度测试的方法并不多:
西北工业大学的徐朝启等人在《固体推进剂燃烧温度的双波长测试方法》提出了一种双波长法测试固体推进剂燃烧温度。选取两个波长区间,得到两个区间光谱辐射量的比值,拟合光谱辐射量比值和温度的关系式,在实验测得光谱辐射量的前提下,通过拟合关系式计算推进剂燃烧温度。这种双波长测温方法仍有一定的缺陷:(1)双波长测温法只能选取两个特定的波长区间,不能表征其他波段的光谱辐射量与温度的关系;(2)该方法主要测量推进剂凝相产物(主要是Al2O3)的光谱辐射量,忽略大量气相产物光谱辐射量的影响,只能适用于特定种类的推进剂;(3)该方法要求使用者预先测定其他推进剂的光谱辐射量和燃气温度,从而拟合出光谱辐射量与温度的关系式,应用于未知推进剂燃烧温度的测量,过程比较繁琐,无法进一步推广应用。
哈尔滨工业大学的孙晓刚等人在《基于多光谱法的固体火箭发动机羽焰温度测量》中,采用8波长高温光谱仪对中国某型号固体火箭发动机羽焰温度进行测量。这种方法的不足是:(1)光谱仪只能采集8个波长的光谱辐射强度数据,数据较少,容易受到非测量对象辐射光谱的干扰;(2)该方法需要预先假设被测对象的发射率模型,当被测目标实际的发射率模型与假定情形不相符时,误差很大。
发明内容
为解决上述问题,本发明提供了一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法,基于普朗克定律,以波长、温度和光谱辐射强度的关系作为数据处理依据,采集200nm-2000nm范围内光谱波长和对应的光谱辐射强度,建立波长和光谱辐射强度的数学模型,根据建立的数学模型计算火箭发动机燃气温度。
考虑多波长光谱辐射的火箭发动机燃气温度测试方法基本原理如下:
对于黑体,普朗克定律表达式如下:
式中:
Eb(λ,T)—-黑体的光谱辐射强度,[W/m2·μm];
T——绝对温度,[K];
λ——光谱辐射波长,[μm];
c1——第一辐射常数,
c2——第二辐射常数,c2=hc0/k=1.4388×104[μm·K];
h——普朗克常数,h=6.626176×10-34[J·s];
k——玻尔兹曼常数,k=1.380662×10-23[J/K];
c0——真空中的光速,c0=2.99792458×1014[μm·s];
根据普朗克定律,预先计算获得每一个温度点(间隔1K取一个温度测量点)下波长与光谱辐射强度的普朗克曲线;然后利用多波长光谱仪实测燃气的波长与对应光谱辐射强度数据,拟合得到实测关系曲线;将实测关系曲线与不同温度下的普朗克曲线进行匹配,当两者相一致时,即认为实测关系曲线所代表的燃气温度为普朗克曲线所表征的温度。
本发明采取的技术方案为:
一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法,包括如下步骤:
S1、构建普朗克曲线数据库:采用编程的方式,根据普朗克定律,预先计算获得1000K-5000K范围内每一个温度点(间隔1K取一个温度测量点)下波长与对应光谱辐射强度的普朗克曲线,波长范围是200nm-2000nm;
S2、多波长光谱仪的标定:考虑到光谱经过发动机玻璃窗口、光学镜头和光纤传递过程中的能量损失,以及光谱仪采集处理信号的失真因素,实际数据处理过程中,1000K-3000K温度范围内,采用黑体炉对多波长光谱仪进行标定,获得光谱辐射强度的损失函数A(λ),3000K以上,损失函数A(λ)采用外推法处理,从而得到普朗克定律修正关系式,如下:
式(2)中,A(λ)是只与波长有关而与温度无关的损失函数,它与该波长下光谱仪的响应率、光学元件透过率、几何尺寸及第一辐射常数有关。
S3、利用多波长光谱仪,采集火箭发动机燃气辐射光谱中多个波长和相应的光谱辐射强度数据,λi和Ea(λi,Ti),波长范围限定于200nm-2000nm;
S4、根据标定的损失函数A(λ)和公式(2),计算获取光谱辐射强度的补偿值Ec(λi,Ti),计算式为Ec(λ,T)=Ea(λ,T)/A(λ);
S5、结合实测波长λi和光谱辐射强度补偿值Ec(λi,Ti),拟合光谱辐射强度公式Ec(λ,T)=g(λ,T),从而得到燃气辐射光谱波长与对应光谱辐射强度的实测关系曲线;
S6、基于最小二乘原理,将光谱波长与对应光谱辐射强度的实测关系曲线与不同温度下的普朗克曲线进行匹配,当曲线的吻合程度最高时,即认为实测温度为普朗克曲线所表征的温度。
本发明具有以下有益效果:
本发明在国内首次提出,基于普朗克定律,采集燃气辐射光谱波长与对应的光谱辐射强度数据,拟合得到实测关系曲线,将实测关系曲线与不同温度下的普朗克曲线进行匹配,根据匹配程度,确定火箭发动机燃气实际温度。
本发明考虑了多波长光谱仪的温度范围是1000K-5000K,这个温度范围基本涵盖了火箭发动机燃气温度。相比于热电偶测温,本发明所能测试温度最高可达5000K,远远超过热电偶的测温范围,更适应于火箭发动机燃气高温测量。
本发明考虑了燃气辐射光谱波长范围是200nm-2000nm。针对燃烧室内气相产物和凝相产物共存的复杂的燃烧环境,本发明能够采集多达数百个波长及所对应的光谱辐射强度数据,提高了波长和光谱辐射强度数据拟合的精度,具有较强的抗干扰能力。相比于双波长测温法和8波长光谱仪测温,能更准确计算燃气温度,适用性更强。
本发明采用黑体炉在3000K以下对光谱仪进行标定。本发明综合考虑了数据采集过程中光谱能量损失和各元件带来的误差,对所采集的光谱辐射强度进行补偿,减小了因设备原因导致的光谱辐射强度损失。
附图说明
图1是本发明的实验装置示意图;
图中:1-火箭发动机;2-透明窗;3-光学镜头;4-多波长光谱仪;
5-计算机数据处理系统。
图2是本发明的数据采集处理流程图。
具体实施方式
为了使本发明的目的及优点更加清楚明白,以下结合实施例对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例
如图1所示,本实施例中的一种测试火箭发动机燃气温度的实验装置,主要包括:火箭发动机1,透明窗2,光学镜头3,多波长光谱仪4,计算机数据处理系统5。
如图2所示,具体的实施步骤如下:
(1)多波长光谱仪的标定。光谱辐射源是最高温度3000K的黑体炉,采用图1所示的光谱采集处理系统,在1000K-3000K的温度范围内,间隔1K取一个温度测量点,采集每一个温度点下波长λi与对应光谱辐射强度Ea(λi,Ti),将实测光谱辐射强度Ea(λi,Ti)与相应波长下普朗克定律计算得到的光谱辐射强度Eb(λi,Ti)进行对比分析,得到光谱辐射强度的损失函数A(λ),在3000K以上,损失函数A(λ)采用外推法处理;
(2)如图1所示,搭建实验装置系统。将光学镜头和光谱仪之间通过光纤连接,将光谱仪和计算机数据处理系统之间通过数据线连接。光学镜头应对准透明窗,保持同轴,以便光学镜头可以捕捉到燃气辐射的光谱。
(3)发动机点火,同时光谱仪开始采集燃气波长和光谱辐射强度数据。燃气辐射的光谱通过发动机透明窗、光学镜头以及光纤传递进入多波长光谱仪中,进行分光处理,光谱仪采集各个波长数据λi及其对应的光谱辐射强度数据Ea(λi,Ti)。
(4)数据处理,计算出发动机燃气温度。根据实际采集到的光谱辐射强度Ea(λi,Ti),结合光谱辐射强度损失函数A(λ),计算获取光谱辐射强度的补偿值Ec(λi,Ti);根据所采集光谱波长λi和光谱辐射强度的补偿值Ec(λi,Ti)的一一对应关系,拟合出波长和光谱辐射强度的关系式Ec(λ,T)=g(λ,T);将关系式g(λ,T)所对应的实测关系曲线与不同温度下的普朗克曲线进行匹配,当两者吻合程度最高时,普朗克曲线所代表的温度即为火箭发动机燃气温度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法,其特征在于,包括如下步骤:
S1、构建普朗克曲线数据库:根据普朗克定律,采用编程的方式,在200nm-2000nm的波长范围内,预先计算获得1000K-5000K范围内每一个温度点下波长与对应光谱辐射强度的普朗克曲线;
S2、多波长光谱仪的标定:考虑到光谱强度在数据采集、处理过程中的损失,在1000K-3000K温度范围内,采用黑体炉对多波长光谱仪进行标定,获得光谱辐射强度的损失函数,3000K以上,损失函数采用外推法处理,进而得到普朗克定律修正关系式;
S3、在200nm-2000nm的波长范围内,采集火箭发动机燃气辐射光谱中多个波长和相应的光谱辐射强度数据;
S4、根据标定的损失函数,计算获取光谱辐射强度的补偿值;
S5、根据燃气辐射光谱波长和对应的辐射强度补偿值,拟合实测光谱辐射强度公式,得到燃气辐射光谱波长与对应光谱辐射强度的实测关系曲线;
S6、基于最小二乘原理,将光谱波长与对应光谱辐射强度的实测关系曲线与不同温度下的普朗克曲线进行匹配,当曲线的吻合程度最高时,即认为实测温度为普朗克曲线所表征的温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410820718.5A CN104864977B (zh) | 2014-12-17 | 2014-12-17 | 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410820718.5A CN104864977B (zh) | 2014-12-17 | 2014-12-17 | 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104864977A CN104864977A (zh) | 2015-08-26 |
CN104864977B true CN104864977B (zh) | 2018-02-06 |
Family
ID=53910965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410820718.5A Expired - Fee Related CN104864977B (zh) | 2014-12-17 | 2014-12-17 | 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104864977B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107843405B (zh) * | 2016-09-21 | 2019-06-11 | 北京空天技术研究所 | 试验件和发动机燃气对飞行器底部辐射热流的获取方法 |
CN109540297B (zh) * | 2018-10-23 | 2021-01-08 | 昆山优尼电能运动科技有限公司 | 基于fpa温度的红外热像仪标定方法 |
CN110095248B (zh) * | 2019-04-28 | 2020-08-28 | 中国科学院力学研究所 | 一种高频感应风洞流场非平衡性诊断系统和方法 |
CN111551255B (zh) * | 2020-04-01 | 2021-05-18 | 华北电力大学 | 一种基于多光谱测量生物质火焰发射率的方法 |
CN112964662B (zh) * | 2021-02-07 | 2022-07-15 | 中国科学院长春光学精密机械与物理研究所 | 一种航空发动机高温燃气浓度及温度测量方法 |
CN115979425B (zh) * | 2022-12-21 | 2023-08-15 | 中国矿业大学 | 一种多波长移动窄带窗口寻优光谱测温方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86105818A (zh) * | 1985-07-15 | 1987-06-10 | 凯诺股份公司 | 物质、材料和物体的特性鉴定和检验的方法和装置 |
CN1712915A (zh) * | 2004-06-21 | 2005-12-28 | 何民才 | 动态辐射测温技术 |
CN1995971A (zh) * | 2006-12-19 | 2007-07-11 | 中国科学院安徽光学精密机械研究所 | 高精度多功能宽光谱辐射比较系统 |
CN103900723A (zh) * | 2014-04-25 | 2014-07-02 | 刘友祥 | 利用辐射谱的全部信息测火焰实际温度的方法及测量系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060190211A1 (en) * | 2001-07-23 | 2006-08-24 | Schietinger Charles W | In-situ wafer parameter measurement method employing a hot susceptor as radiation source for reflectance measurement |
-
2014
- 2014-12-17 CN CN201410820718.5A patent/CN104864977B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86105818A (zh) * | 1985-07-15 | 1987-06-10 | 凯诺股份公司 | 物质、材料和物体的特性鉴定和检验的方法和装置 |
CN1712915A (zh) * | 2004-06-21 | 2005-12-28 | 何民才 | 动态辐射测温技术 |
CN1995971A (zh) * | 2006-12-19 | 2007-07-11 | 中国科学院安徽光学精密机械研究所 | 高精度多功能宽光谱辐射比较系统 |
CN103900723A (zh) * | 2014-04-25 | 2014-07-02 | 刘友祥 | 利用辐射谱的全部信息测火焰实际温度的方法及测量系统 |
Non-Patent Citations (3)
Title |
---|
"多光谱辐射测温技术测量火工烟火药剂燃烧温度";李占英等;《光谱学与光谱分析》;20101231;第30卷(第8期);第2062-2064页 * |
"太阳光谱测温处理方法与校正技术的研究";王丽玲;《中国优秀硕士论文全文数据库》;20140815(第8期);正文第13、24-28页 * |
基于辐射强度多波长分析的燃烧火焰温度测量方法的实验研究;周洁;《动力工程》;19991231;第19卷(第6期);第460-465页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104864977A (zh) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104864977B (zh) | 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法 | |
Mekhrengin et al. | Multispectral pyrometer for high temperature measurements inside combustion chamber of gas turbine engines | |
CN107152972B (zh) | 一种航空发动机涡轮叶片温度监测装置 | |
CN102830064B (zh) | 一种中高温红外发射率测试装置 | |
Fu et al. | Temperature measurements of diesel fuel combustion with multicolor pyrometry | |
Liang et al. | Development of a new fiber-optic multi-target multispectral pyrometer for achievable true temperature measurement of the solid rocket motor plume | |
CN109655161B (zh) | 基于红外热像仪的目标红外积分辐射强度测试方法及设备 | |
Lee et al. | Calibration of a mid-IR optical emission spectrometer with a 256-array PbSe detector and an absolute spectral analysis of IR signatures | |
CN101793563A (zh) | 多波段红外辐射自动测量系统 | |
Lü et al. | Ambient temperature-independent dual-band mid-infrared radiation thermometry | |
CN105784333A (zh) | 红外窗口的高温光谱发射率测量方法及装置 | |
Fu et al. | Temperature measurements of high-temperature semi-transparent infrared material using multi-wavelength pyrometry | |
CN107328478B (zh) | 一种基于三波段辐射光测量涡轮叶片温度及发射率的方法 | |
CN102620841A (zh) | 基于光谱技术的瞬态温度测量与反演系统 | |
CN106482790B (zh) | 基于火焰辐射的固体火箭推进剂燃烧测量装置和测量方法 | |
Sankaranarayanan et al. | Investigation of sooting flames by color-ratio pyrometry with a consumer-grade DSLR camera | |
Wang et al. | The development of a multispectral pyrometer for achievable true temperature field measurements of the explosion flame | |
CN105890776A (zh) | 一种分光耦合装置 | |
CN103644972B (zh) | 一种多光谱辐射测温装置及方法 | |
Deep et al. | Temperature characterization of a radiating gas layer using digital-single-lens-reflex-camera-based two-color ratio pyrometry | |
CN201892573U (zh) | 一种近红外辐射温度计 | |
CN112964367B (zh) | 一种航空发动机涡轮叶片辐射测温方法 | |
CN106768366A (zh) | 一种利用微型光谱仪测量高温的方法及装置 | |
Sun et al. | Research on error correction method of filter assembly based on abnormal trend recognition of sequence data | |
Liang et al. | Direct calibration method of effective wavelength for fiber-optic multi-wavelength pyrometer based on optimization function |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180206 Termination date: 20181217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |