CN109974893B - 一种梯度折射率火焰三维温度场测量方法 - Google Patents
一种梯度折射率火焰三维温度场测量方法 Download PDFInfo
- Publication number
- CN109974893B CN109974893B CN201910241921.XA CN201910241921A CN109974893B CN 109974893 B CN109974893 B CN 109974893B CN 201910241921 A CN201910241921 A CN 201910241921A CN 109974893 B CN109974893 B CN 109974893B
- Authority
- CN
- China
- Prior art keywords
- flame
- refractive index
- distribution
- radiation
- radiation intensity
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/006—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of the effect of a material on microwaves or longer electromagnetic waves, e.g. measuring temperature via microwaves emitted by the object
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明介绍一种梯度折射率火焰高精度三维温度场测量方法,包括以下步骤,利用“Z”字型纹影技术对轴对称火焰成像并定量化得到光线在折射率场中的偏转,由起偏角计算出火焰内部的折射率场;遵循菲涅尔折射定律逆向追踪光线,然后根据梯度折射率介质下的辐射传输理论,结合源项有限体积法建立火焰广义辐射源项和火焰出射辐射强度之间的关系;最后采用LSQR反演算法迭代计算,得到收敛解,重建出火焰的三维温度场分布。本方法基于“Z”字型纹影技术得到的火焰内部非均匀分布折射率场,考虑了梯度折射率介质下火焰对光的折射作用,从而可以实现精度更高的火焰三维温度场测量。
Description
技术领域
本发明涉及一种梯度折射率火焰高精度三维温度场测量方法,属于测量技术领域。
背景技术
在高温工业设备中,如锅炉、燃气轮机、发动机燃烧室等,为了提高燃烧效率和加强燃烧安全防护,通常需要实时监测各项燃烧参数。温度是高温热处理的重要热工参数,实时而准确的温度场检测是预防工业事故、高温测量、低碳环保、节约能源的重要手段。
图像法作为一种被动式非接触测量手段,在高温火焰的温度监测中得到了越来越多的关注,它需要建立火焰中的辐射传输模型。目前图像测温法通常假定介质的折射率均匀分布,即假定了光线的传播轨迹可以按直线处理,然而在实际生产中,火焰内部的温度分布不均会导致折射率分布不均。尤其是在大型火焰中,火焰内部存在多个高温区,如果不考虑由于折射率分布不均带来的光线偏折的影响,由此得到的温度测量结果会造成较大的误差。
折射率的测量方法通常有纹影法、光学干涉法、阴影法、全息干涉法等,而其中纹影法具有测量精度高、测量装置简单的优点,它是一种将相分布转换为可见图像的光学方法,通常可以用于测量一个区域内的介质折射率,然后根据折射率与温度的G-D关系式来测量火焰的温度场,但不同燃料和燃烧工况下的温度和折射率之间的这种状态关系是半经验的理论,只能针对少量工况有较好的测量精度。
发明内容
本发明所要解决的技术问题是针对上述现有技术存在的问题,而提出一种基于纹影法与图像法的融合测试技术,旨在测量梯度折射率火焰的三维温度场,以解决传统火焰温度场计算模型因忽略火焰对光的折射作用而无法保证足够精度的问题。
本发明的技术方案如下:
一种梯度折射率火焰高精度三维温度场测量方法,包括以下几个步骤:
步骤一:利用同步控制器连接两台相机,对轴对称燃烧火焰实现同时采集其“Z”字型纹影成像和火焰自发光灰度图像。
步骤二:根据纹影成像中光线经过火焰由于折射率引起的起偏角,获得火焰的折射率分布,偏转角θ和折射率n之间的关系:
其中x方向为光线入射方向,y方向为入射方向的垂直方向。
步骤三:根据火焰的折射率分布,利用菲涅尔折射定律,从相机的探测器逆向追踪光线,记录光线在火焰中的穿行方向和控制体编号。
步骤四:利用黑体炉在相同工作距离和曝光条件下标定图像灰度与温度之间的关系,获得火焰的出射辐射强度分布。
步骤五:根据光线在火焰中的传播路径,利用源项有限体积法建立火焰广义辐射源项和相机探测到的火焰出射辐射强度之间的关系:
其中,I(rw,Ω)表示探测到的火焰出射辐射强度,S(r,Ω)表示火焰广义辐射源项,β表示火焰的衰减系数,M表示光线在火焰中穿过的控制体总个数,Δs表示火焰在每个控制体中穿行的距离,Ω表示某根光线代表的探测方向。
步骤六:根据广义辐射源项和火焰出射辐射强度的关系,采用LSQR反演算法进行迭代计算,重建出火焰的三维辐射源项分布。
步骤七:将广义源项分布代入梯度折射率介质的辐射传递方程,利用直角坐标系下的有限体积法计算出梯度折射率火焰的辐射强度分布,直角坐标系下梯度折射率辐射传递方程形式为:
步骤八:根据广义源项的定义和辐射强度分布,计算出火焰的自发辐射强度,进而根据黑体发光或化学发光中温度与自发辐射强度之间的关系获得火焰的温度分布。广义辐射源项定义为:
其中,Ib(r)表示火焰的自发辐射强度分布,I(r,Ω′)表示火焰在位置r沿方向Ω′的辐射强度大小Φ(Ω,Ω′)表示由Ω方向散射出去并且由Ω′方向入射的散射相函。dΩ′表示立体角大小。
有益效果
与现有技术相比,本发明有以下优点:
1、考虑了火焰的折射率分布,建立了梯度折射率火焰的成像模型,比当前未考虑折射率影响的图像法温度测量方法,从模型上来说更加精确。
2、利用纹影法测量得到了火焰的折射率分布,没有进一步采用传统纹影法测温的折射率和温度场的关系或者温度场,而是与图像法测温相结合,避免了传统的纹影法测温时由于使用半经验公式而只能在少数工况下适用的范围限制。
3、由于采用纹影法与图像法相融合的测温方法,突破了原来纹影法测温中只考虑气体而忽略颗粒散射的特性,测温对象可以从原来的化学发光火焰拓展到黑体发光与化学发光共存的大部分火焰。
附图说明
图1为实验系统示意图;
图2为采集的蜡烛火焰图像;
图3为采用纹影法得到的火焰纹影图;
图4为在某一控制体界面处的光线偏折示意图;
图5为光线在梯度折射率介质中的轨迹追踪示意图;
图6为黑体炉标定的灰度值与强度值的关系曲线图;
图7为考虑梯度折射率后的火焰温度场三维重建结果。
具体实施方式
下面结合附图,对本发明的具体实施方式作进一步详细描述。
具体包括以下几个步骤:
步骤一:搭建如图1所示的实验系统,首先利用同步控制器连接纹影相机以及图像法测温相机,然后对燃烧的蜡烛火焰实现同时采集其“Z”字型纹影成像和火焰自发光灰度图像(如图2,图3所示)。
步骤二:根据纹影成像中光线经过火焰由于折射率引起的起偏角,根据偏转角θ和折射率n之间的关系,获得火焰的折射率分布,
其中,x方向为光线入射方向,y方向为入射方向的垂直方向。
步骤三:根据火焰的折射率分布,利用菲涅尔折射定律,从相机的探测器逆向追踪光线,记录光线在火焰中的穿行方向和控制体编号,菲涅尔折射定律为:
式中,j为控制体编号;
光线在其中某一个控制体界面的偏折如图4所示;而当用I(rw,Ω)表示探测到的火焰出射辐射强度,Δs表示火焰在每个控制体中穿行的距离时,光线整个的偏折轨迹如图5所示。
步骤四:打开黑体炉和相机,将相机与黑体炉的距离、曝光条件等参数调整到与采集火焰图像时相同,黑体炉温度从600℃到850℃,每隔50℃拍摄十张黑体炉图像,根据普朗克定律,通过拟合得到强度值与灰度值的关系,完成标定(如图6),获得火焰的出射辐射强度分布。
步骤五:根据光线在火焰中的传播路径,利用源项有限体积法建立火焰广义辐射源项和相机探测到的火焰出射辐射强度之间的关系:
其中,I(rw,Ω)表示探测到的火焰出射辐射强度,S(r,Ω)表示火焰广义辐射源项,β表示火焰的衰减系数,M表示光线在火焰中穿过的控制体总个数,Δs表示火焰在每个控制体中穿行的距离,Ω表示某根光线代表的探测方向;i,j为控制体编号;
步骤六:在吸收系数已知的条件下,步骤五中的关系式可简化为如下形式:
AX=B
其中,A,B均为常数矩阵,X为所求辐射源项分布矩阵。
接下来采用LSQR反演算法进行计算,通过初始化、变量迭代循环计算以及参数修正等步骤对该方程进行求解,最终求得火焰的三维辐射源项分布结果。
步骤七:将广义源项分布代入梯度折射率介质的辐射传递方程,利用直角坐标系下的有限体积法计算出梯度折射率火焰的辐射强度分布,直角坐标系下梯度折射率辐射传递方程形式为:
步骤八:根据广义源项的定义和辐射强度分布,计算出火焰的自发辐射强度,进而根据黑体发光或化学发光中温度与自发辐射强度之间的关系获得火焰的温度分布(如图5所示)。广义辐射源项定义为:
其中,Ib(r)表示火焰的自发辐射强度分布,I(r,Ω′)表示火焰在位置r沿方向Ω′的辐射强度大小Φ(Ω,Ω′)表示由Ω方向散射出去并且由Ω′方向入射的散射相函;dΩ′表示立体角大小。
研究对象为黑体发光火焰时,辐射强度与温度的关系为:
式中,Ibλ是波长λ的光谱辐射强度值,W/(sr·m3);ε是黑体的发射率;T是黑体的温度值,K;c1为第一辐射常数,c2为第二辐射常数,分别为3.7418×10-16W·m2和1.4388×10- 2m·K。
研究对象为化学发光火焰时,根据化学发光测温原理,分子中不同能级间的跃迁产生相应的谱线强度与温度的关系为:
式中,h为普朗克常数,c为光速,v为振动能级,J′为高能级转动量子数;N为高能级分子数,A为自发发射爱因斯坦跃迁概率,v为跃迁波数;N0为分子总数,k为玻尔兹曼常数,ge为电子态简并度,Te、Tv、Tr分别为电子、振动和转动温度,Ee、Ev、Er分别为电子态、振动态和转动态能量,Qe、Qv、Qr分别为电子态、振动态和转动态配分函数。
Claims (4)
1.一种梯度折射率火焰三维温度场测量方法,其特征在于,包括以下步骤:
步骤一:利用同步控制器连接两台相机,对轴对称燃烧火焰实现同时采集燃烧火焰的“Z”字型纹影成像和火焰自发光灰度图像;
步骤二:根据纹影成像得到光线在火焰中的传播路径;根据火焰自发光灰度图像获得火焰的出射辐射强度分布;
步骤三:根据光线在火焰中的传播路径,利用源项有限体积法建立火焰广义辐射源项和相机探测到的火焰出射辐射强度之间的关系;
步骤四:根据广义辐射源项和火焰出射辐射强度的关系,采用LSQR反演算法进行迭代计算,重建出火焰的三维辐射源项分布;
步骤五:根据重建得到的火焰的三维辐射源项分布获得火焰的温度分布;
步骤三中,建立的火焰广义辐射源项和相机探测到的火焰出射辐射强度之间的关系为:
其中,I(rw,Ω)表示探测到的火焰出射辐射强度,S(r,Ω)表示火焰广义辐射源项,β表示火焰的衰减系数,M表示光线在火焰中穿过的控制体总个数,Δs表示火焰在每个控制体中穿行的距离,Ω表示某根光线代表的探测方向;
所述火焰广义辐射源项定义为:
其中,Ib(r)表示火焰的自发辐射强度分布;I(r,Ω′)表示火焰在位置r沿方向Ω′的辐射强度大小Φ(Ω,Ω′)表示由Ω方向散射出去并且由Ω′方向入射的散射相函;dΩ′表示立体角大小;
ka为吸收系数;ks表示散射系数;n表示折射率;
所述步骤五,根据重建得到的火焰的三维辐射源项分布获得火焰的温度分布的方法是:
将广义源项分布代入梯度折射率介质的辐射传递方程,利用直角坐标系下的有限体积法计算出梯度折射率火焰的辐射强度分布;
根据广义源项的定义和辐射强度分布,计算出火焰的自发辐射强度;根据温度与自发辐射强度之间的关系获得火焰的温度分布;
所述直角坐标系下梯度折射率辐射传递方程形式为:
其中,I(r,Ω)表示火焰在位置r沿方向Ω的辐射强度大小;
2.根据权利要求1所述的火焰三维温度场测量方法,其特征在于:所述步骤二中,根据纹影成像得到光线在火焰中的传播路径的方法是:
根据纹影成像中光线经过火焰由于折射率引起的起偏角,获得火焰的折射率分布;
根据火焰的折射率分布,利用菲涅尔折射定律,从相机的探测器逆向追踪光线,记录光线在火焰中的穿行方向和控制体编号:
控制体沿着某个方向串起来一条曲线得到传播路径。
4.根据权利要求1所述的火焰三维温度场测量方法,其特征在于:所述步骤二中,根据火焰自发光灰度图像获得火焰的出射辐射强度分布的方法是:
利用黑体炉在相同工作距离和曝光条件下标定图像灰度与温度之间的关系,获得火焰的出射辐射强度分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910241921.XA CN109974893B (zh) | 2019-03-27 | 2019-03-27 | 一种梯度折射率火焰三维温度场测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910241921.XA CN109974893B (zh) | 2019-03-27 | 2019-03-27 | 一种梯度折射率火焰三维温度场测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109974893A CN109974893A (zh) | 2019-07-05 |
CN109974893B true CN109974893B (zh) | 2020-11-03 |
Family
ID=67081052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910241921.XA Active CN109974893B (zh) | 2019-03-27 | 2019-03-27 | 一种梯度折射率火焰三维温度场测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109974893B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110595632A (zh) * | 2019-08-26 | 2019-12-20 | 南京理工大学 | 一种集成阵列式热电偶束的三维立体测温装置及方法 |
CN116576983B (zh) * | 2023-05-12 | 2023-10-27 | 南京航空航天大学 | 一种电弧放电的二维温度测量方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10185694A (ja) * | 1996-12-20 | 1998-07-14 | Tokyo Gas Co Ltd | 火炎温度測定方法およびその装置 |
JP2012013580A (ja) * | 2010-07-01 | 2012-01-19 | Central Res Inst Of Electric Power Ind | 粒子及び液滴の形状・径と温度との同時計測システム並びにそのプログラム |
CN102353478A (zh) * | 2011-10-10 | 2012-02-15 | 哈尔滨工业大学 | 半透明介质环境下非接触测温的校正方法 |
CN102393027A (zh) * | 2011-10-18 | 2012-03-28 | 沈阳理工大学 | 一种炉膛三维温度场重建方法 |
CN105547469A (zh) * | 2015-12-04 | 2016-05-04 | 哈尔滨工业大学 | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 |
CN105606222A (zh) * | 2015-09-06 | 2016-05-25 | 东南大学 | 一种火焰三维温度场测量的成像装置、测量装置及测量方法 |
CN105866970A (zh) * | 2016-05-27 | 2016-08-17 | 武汉理工大学 | 可矫正色差的反射式定容燃烧弹高速摄像纹影系统 |
CN108731837A (zh) * | 2018-05-21 | 2018-11-02 | 中国计量大学 | 一种双光路光学结构的火焰测温系统的测量方法 |
CN108917943A (zh) * | 2018-08-14 | 2018-11-30 | 岭南师范学院 | 一种基于投影式背景纹影技术的火焰温度场测量装置及方法 |
-
2019
- 2019-03-27 CN CN201910241921.XA patent/CN109974893B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10185694A (ja) * | 1996-12-20 | 1998-07-14 | Tokyo Gas Co Ltd | 火炎温度測定方法およびその装置 |
JP2012013580A (ja) * | 2010-07-01 | 2012-01-19 | Central Res Inst Of Electric Power Ind | 粒子及び液滴の形状・径と温度との同時計測システム並びにそのプログラム |
CN102353478A (zh) * | 2011-10-10 | 2012-02-15 | 哈尔滨工业大学 | 半透明介质环境下非接触测温的校正方法 |
CN102393027A (zh) * | 2011-10-18 | 2012-03-28 | 沈阳理工大学 | 一种炉膛三维温度场重建方法 |
CN105606222A (zh) * | 2015-09-06 | 2016-05-25 | 东南大学 | 一种火焰三维温度场测量的成像装置、测量装置及测量方法 |
CN105547469A (zh) * | 2015-12-04 | 2016-05-04 | 哈尔滨工业大学 | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 |
CN105866970A (zh) * | 2016-05-27 | 2016-08-17 | 武汉理工大学 | 可矫正色差的反射式定容燃烧弹高速摄像纹影系统 |
CN108731837A (zh) * | 2018-05-21 | 2018-11-02 | 中国计量大学 | 一种双光路光学结构的火焰测温系统的测量方法 |
CN108917943A (zh) * | 2018-08-14 | 2018-11-30 | 岭南师范学院 | 一种基于投影式背景纹影技术的火焰温度场测量装置及方法 |
Non-Patent Citations (2)
Title |
---|
基于光场成像的火焰三维温度场测量方法研究;孙俊;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20190115(第12期);第33、43-44页 * |
梯度折射率介质内辐射传递方程数值模拟的有限元法;张琳;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20111115(第11期);第11-12页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109974893A (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhao et al. | Optical sectioning tomographic reconstruction of three-dimensional flame temperature distribution using single light field camera | |
Hossain et al. | Three-dimensional reconstruction of flame temperature and emissivity distribution using optical tomographic and two-colour pyrometric techniques | |
Sun et al. | Three-dimensional temperature field measurement of flame using a single light field camera | |
Sun et al. | A novel calibration method of focused light field camera for 3-D reconstruction of flame temperature | |
CN102706459B (zh) | 一种单ccd成像系统的炉膛内三维温度场检测装置及方法 | |
Zhou et al. | 3-D reconstruction of flame temperature distribution using tomographic and two-color pyrometric techniques | |
Li et al. | Estimation of radiative properties and temperature distributions in coal-fired boiler furnaces by a portable image processing system | |
Li et al. | Simultaneous measurement of flame temperature and absorption coefficient through LMBC-NNLS and plenoptic imaging techniques | |
CN109974893B (zh) | 一种梯度折射率火焰三维温度场测量方法 | |
Zhang et al. | An inverse method for flue gas shielded metal surface temperature measurement based on infrared radiation | |
Zhang et al. | An improved colorimetric method for visualization of 2-D, inhomogeneous temperature distribution in a gas fired industrial furnace by radiation image processing | |
Li et al. | Simulation of calibration process in flame measurement by plenoptic camera | |
Niu et al. | Temperature field reconstruction of 3D luminous flames based on light field tomography theory | |
Qi et al. | Flame temperature reconstruction through a multi-plenoptic camera technique | |
Liu et al. | Volumetric imaging of flame refractive index, density, and temperature using background-oriented Schlieren tomography | |
CN105973473A (zh) | 基于目标表面温度分布的空间温度场重构方法及装置 | |
Qi et al. | Soot temperature measurement within 3D flame by light-field imaging based on wave optics theory | |
Shi et al. | Three-dimensional temperature reconstruction of diffusion flame from the light-field convolution imaging by the focused plenoptic camera | |
Li et al. | Effect of nonuniform radiation properties on flame temperature reconstruction based on light field imaging | |
Zhang et al. | Development of a CCD-based pyrometer for surface temperature measurement of casting billets | |
Li et al. | Meshed axisymmetric flame simulation and temperature reconstruction using light field camera | |
Hossain et al. | Three-dimensional reconstruction of flame temperature and emissivity through tomographic imaging and pyrometric measurement | |
Dai et al. | Experimental and simulation investigation of 3-D soot temperature and volume fraction fields of afterburner flame | |
Wang et al. | Reconstruction of temperature field in coal-fired boiler based on limited flame image information | |
Shi et al. | Reconstruction of 3D flame temperature and absorption coefficient field by the hybrid light-field imaging and laser extinction technique |
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 |