CN103344339A - 一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 - Google Patents
一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 Download PDFInfo
- Publication number
- CN103344339A CN103344339A CN2013102306248A CN201310230624A CN103344339A CN 103344339 A CN103344339 A CN 103344339A CN 2013102306248 A CN2013102306248 A CN 2013102306248A CN 201310230624 A CN201310230624 A CN 201310230624A CN 103344339 A CN103344339 A CN 103344339A
- Authority
- CN
- China
- Prior art keywords
- temperature
- value
- gray
- formula
- internal combustion
- 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.)
- Granted
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及计算温度场分布的方法,根据内燃机缸内高速摄影得到的照片和压力示功图计算修正缸内温度场分布的方法,主要用于分析内燃机缸内燃烧过程中温度场的分布,以及不同温度区域的面积,为进一步分析内燃机燃烧和排放特性提供定量数据参考。本发明的优点是:不需要复杂的测试硬件,能批量处理各曲轴转角对应的燃烧图像,能很简便快速的得出温度场的分布,并能很好的验证计算的准确性。
Description
技术领域
本发明专利涉及计算温度场分布的方法,特指一种根据内燃机缸内高速摄影得到的照片和压力示功图计算修正缸内温度场分布的方法,主要用于分析内燃机缸内燃烧过程中温度场的分布,以及不同温度区域的面积,为进一步分析内燃机燃烧和排放特性提供定量数据参考。
背景技术
内燃机燃烧过程中温度场的分布情况对内燃机的燃烧效率、零部件热负荷的分布和燃烧产物的生成等有重要影响,因此,温度场的测量和计算已经成为内燃机研究领域的一个热点。随着计算机和光学测量技术的进步,高速摄影技术已经广泛的应用于内燃机的可视化研究中,可以较为方便的获取内燃机缸内燃烧过程的视频。现有的专利大多倾向于搭建温度场测量的平台以及提出与之配套的测量方法,如200910077463.7“一种双CCD温度场测量装置及方法”,200910077464.1“一种三CCD温度场测量装置及方法”、02292184.2“发动机燃烧室温度场旋转测试装置”、94111569.0“电站锅炉炉膛燃烧温度分布检测方法及其装置”、01118260.1“测量锅炉燃烧辐射能及温度场并控制燃烧的方法及其***”等,没有涉及利用高速摄影照片来获取温度场分布的问题。
利用燃烧分析仪可以测得内燃机缸内燃烧压力,并获得示功图,该示功图能够反映燃烧温度与曲轴转角的关系曲线,即内燃机缸内压力随曲轴转角的变化关系,燃烧分析仪的后处理软件,可在缸内压力曲线的基础上,通过设定燃料特性、初始条件、边界条件,得到缸内平均温度随曲轴转角的变化关系。该方法得到的是缸内的平均温度,不能获得缸内温度场的分布情况。
发明内容
本发明的旨在提供一种基于内燃机缸内高速摄影照片计算温度场分布的方法。
本发明是利用内窥镜高速摄影装置,获取内燃机缸内燃烧火焰的照片,提取燃烧过程高速摄影照片上各点的灰度值,根据灰度值与温度的对于关系,计算出各点的温度值,从而得出照片区域温度场的分布。
本发明的基于内燃机缸内燃烧火焰摄影照片计算温度场的方法,包括如下步骤:
步骤1、拍摄内燃机缸内燃烧火焰并获得彩色照片,记录照片对应的曲轴转角度数为通过燃烧分析仪测得内燃机缸内燃烧压力,获得示功图,根据示功图的燃烧温度与曲轴转角的关系曲线,得到曲轴转角所对应的温度,记为Tmc;
步骤2、在彩色照片中,像素点的数量为n,记像素点的灰度值为N,记灰度值最大点的灰度值为NH,则NH=256;
步骤3、将灰度值最大点对应的温度值在1800k~2300k区间内取任意值,记为TH;
步骤4、根据公式(1),
式中:NL为低温点的灰度值,NH为灰度值最大点的灰度值,F(N)为修正系数;
步骤5、根据公式(2)计算低温点的温度值TL:
式中:TL为低温点的温度值,TH为灰度值最大点的温度值;
根据公式(2)计算获得TL,至此获得两组数据(NH,TH)、(NL,TL);
步骤6、根据公式(3)计算:
T=a×Nb …………(3)
式中:N为像素点的灰度值,T为像素点灰度值对应的温度值,a、b分别为常数,
将步骤5获得的两组数据(NH,TH)、(NL,TL)代入公式(3),求解公式(3)获得a和b的值;
步骤7,将步骤2中的n个像素点的灰度值N逐一代入公式(3),计算获得像素点灰度值对应的温度值T,共获得n组灰度值及其对应的温度值数据(N,T);
步骤8、将n个温度值数据分为I个温度区域,记每个温度区域的平均温度为Ti,i=1,2,……I,通过图像处理软件获得每个温度区域内像素点的面积,并根据公式(4)计算获得平均温度Tm:
式中:i为温度区域的序号,i=1,2,……I,I为温度区域的个数,Tm为平均温度,Ti为温度区域的平均温度,Si为温度区域的面积;
步骤9、计算|Tm-Tmc|,若|Tm-Tmc|≥20K,则步骤2中TH在1800k~2300k区间的取值错误,返回步骤2重新取值,若|Tm-Tmc|<20K,则认为步骤2中TH取值正确,进入步骤10;
步骤10、根据公式(3)T=a×Nb,实现彩色照片中任意像素点的灰度值计算获得对应的温度值。
本发明的优点是不需要复杂的测试硬件,能批量处理各曲轴转角对应的燃烧图像,能很简便快速的得出温度场的分布,并能很好的验证计算的准确性。
可将内燃机缸内火焰的辐射当作灰体处理,根据全辐射定律,燃烧火焰的全辐射亮度与其本身的辐射温度存在线性关系,而内燃机燃烧过程高速摄影照片的像素点与亮度灰度值是对应的,因此,可根据高速摄影照片计算出各点的灰度值,进而得到其对应的温度值。这种方法得到的是温度场的分布。
附图说明
图1为本发明引入光源拍摄照片示意图。
图2为本发明拍摄照片时内窥镜在气缸盖上的位置图。
图3为不同曲轴转角拍摄获得的彩色照片。
图4为本发明方法流程图。
图5为采用本发明方法计算后获得的火焰照片和计算的温度场分布拟合效果图。
图中,1、气缸中心,2、燃烧室中心,3、汽缸套,4、引入光源,5、引出光源,6、喷嘴中心,7、内窥镜中心,8、燃烧室中心,9、缸盖底平面。
具体实施方式
本发明是基于内燃机缸内燃烧火焰摄影照片计算温度的方法,以490Q柴油机缸内燃烧火焰的分析来说明本发明的具体实施方式。具体实施包括以下步骤:
步骤1、拍摄内燃机缸内燃烧火焰并获得彩色照片,记录照片对应的曲轴转角度数为通过燃烧分析仪测得内燃机缸内燃烧压力,获得示功图,根据示功图的燃烧温度与曲轴转角的关系曲线,得到曲轴转角所对应的温度,记为Tmc;
如图1所示,本实施例在缸盖上布置了两个内窥镜,向燃烧室内引入光源,能够同时拍摄到喷油油束和燃烧火焰的图像。引入光源为70°视角,引出光源为30°视角。
内窥镜在气缸盖上的位置图如图2所示。摄影装置的主要技术参数如表1所示。试验所用发动机为常州柴油机有限公司生产的CC490Q柴油机,其主要性能参数见表2所示。选择的发动机空载且转速为1500r/min的工况。拍摄得到的部分曲轴转角对应的照片如图3所示。
本实施例中,选取曲轴转角7.9°CA对应的照片为例进行计算。
表1高速摄影装置的主要技术参数
表2柴油机的主要性能参数
步骤2、在彩色照片中,像素点的数量为n,记像素点的灰度值为N,记灰度值最大点的灰度值为NH,则NH=256。
步骤3、将灰度值最大点对应的温度值在1800k~2300k区间内取任意值,记为TH;本实施例中TH=2000k。
步骤4、根据公式(1),
式中:NL为低温点的灰度值,NH为灰度值最大点的灰度值,F(N)为修正系数;
步骤5、根据公式(2)计算低温点的温度值TL:
式中:TL为低温点的温度值,TH为灰度值最大点的温度值;
根据公式(2)计算获得TL,TL=1570k,至此获得两组数据(NH,TH)为(256,2000k)、(NL,TL)为(64,1570k);
步骤6、根据公式(3)计算:
T=a×Nb …………(3)
式中:N为像素点的灰度值,T为像素点灰度值对应的温度值,a、b分别为常数,
将步骤5获得的两组数据(NH,TH)、(NL,TL)代入公式(3),求解公式(3)得到a=759.5408,b=0.1746。
步骤7,将步骤2中的n个像素点的灰度值N逐一代入公式(3),计算获得像素点灰度值对应的温度值T,共获得n组灰度值及其对应的温度值数据(N,T);
步骤8、将n个温度值数据分为I个温度区域,记每个温度区域的平均温度为Ti,i=1,2,……I,通过图像处理软件获得每个温度区域内像素点的面积,并根据公式(4)计算获得平均温度Tm:
式中:i为温度区域的序号,i=1,2,……I,I为温度区域的个数,Tm为平均温度,Ti为温度区域的平均温度,Si为温度区域的面积,单位为mm2。
本实施例中,最低温度从600k开始,最高温度为2300k,以100k为间隔,将温度区域划分为17个,图像处理软件为:MATLAB,美国MathWorks公司出品,2010版,通过软件计算获得温度区域的面积Si,即可通过公式(4)计算得到Tm。
步骤9、计算|Tm-Tmc|,若|Tm-Tmc|≥20K,则认为步骤2中TH在1800k~2300k区间的取值错误,返回步骤2重新取值,若|Tm-Tmc|<20K,则认为步骤2中TH取值正确,进入步骤10;
本实施例中,最后计算得到满足约束条件的灰度值与温度值的关系式如下:
T=754.3278N0.1832 …………(10)整个过程的循环计算流程图如图4所示。
步骤10、根据公式(10)T=754.3278N0.1832,实现彩色照片中任意像素点的灰度值计算获得对应的温度值。火焰照片和计算的温度场分布效果图如图5所示。
Claims (1)
1.一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法,其特征在于包括以下步骤:
步骤1、拍摄内燃机缸内燃烧火焰并获得彩色照片,记录照片对应的曲轴转角度数为通过燃烧分析仪测得内燃机缸内燃烧压力,获得示功图,根据示功图的燃烧温度与曲轴转角的关系曲线,得到曲轴转角所对应的温度,记为Tmc;
步骤2、在彩色照片中,像素点的数量为n,记像素点的灰度值为N,记灰度值最大点的灰度值为NH,则NH=256;
步骤3、将灰度值最大点对应的温度值在1800k~2300k区间内取任意值,记为TH;
步骤4、根据公式(1),
式中:NL为低温点的灰度值,NH为灰度值最大点的灰度值,F(N)为修正系数;
步骤5、根据公式(2)计算低温点的温度值TL:
式中:TL为低温点的温度值,TH为灰度值最大点的温度值;
根据公式(2)计算获得TL,至此获得两组数据(NH,TH)、(NL,TL);
步骤6、根据公式(3)计算:
T=a×Nb …………(3)
式中:N为像素点的灰度值,T为像素点灰度值对应的温度值,a、b分别为常数,
将步骤5获得的两组数据(NH,TH)、(NL,TL)代入公式(3),求解公式(3)获得a和b的值;
步骤7,将步骤2中的n个像素点的灰度值N逐一代入公式(3),计算获得像素点灰度值对应的温度值T,共获得n组灰度值及其对应的温度值数据(N,T);
步骤8、将n个温度值数据分为I个温度区域,记每个温度区域的平均温度为Ti,i=1,2,……I,通过图像处理软件获得每个温度区域内像素点的面积,并根据公式(4)计算获得平均温度Tm:
式中:i为温度区域的序号,i=1,2,……I,I为温度区域的个数,Tm为平均温度,Ti为温度区域的平均温度,Si为温度区域的面积;
步骤9、计算|Tm-Tmc|,若|Tm-Tmc|≥20K,则步骤2中TH在1800k~2300k区间的取值错误,返回步骤2重新取值,若|Tm-Tmc|<20K,则认为步骤2中TH取值正确,进入步骤10;
步骤10、根据公式(3)T=a×Nb,实现彩色照片中任意像素点的灰度值计算获得对应的温度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310230624.8A CN103344339B (zh) | 2013-06-09 | 2013-06-09 | 一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310230624.8A CN103344339B (zh) | 2013-06-09 | 2013-06-09 | 一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103344339A true CN103344339A (zh) | 2013-10-09 |
CN103344339B CN103344339B (zh) | 2015-09-30 |
Family
ID=49279155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310230624.8A Active CN103344339B (zh) | 2013-06-09 | 2013-06-09 | 一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103344339B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846366A (zh) * | 2018-06-21 | 2018-11-20 | 江苏大学 | 计算内燃机缸内浓度场的方法及装置 |
CN109764820A (zh) * | 2018-12-24 | 2019-05-17 | 西华大学 | 一种定容燃烧火焰传播半径的测量角度步长确定方法 |
CN114152347A (zh) * | 2021-09-30 | 2022-03-08 | 国网黑龙江省电力有限公司电力科学研究院 | 一种变电站电力设备故障定位及火险研判综合检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308224A (zh) * | 2001-03-27 | 2001-08-15 | 华中科技大学 | 炉膛燃烧温度图像检测方法 |
US20030123518A1 (en) * | 2002-01-03 | 2003-07-03 | Abbasi Hamid A. | Dual wavelength thermal imaging system for surface temperature monitoring and process control |
CN102798475A (zh) * | 2012-08-02 | 2012-11-28 | 电子科技大学 | 一种根据标准红外视频图像确定目标温度的方法 |
-
2013
- 2013-06-09 CN CN201310230624.8A patent/CN103344339B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308224A (zh) * | 2001-03-27 | 2001-08-15 | 华中科技大学 | 炉膛燃烧温度图像检测方法 |
US20030123518A1 (en) * | 2002-01-03 | 2003-07-03 | Abbasi Hamid A. | Dual wavelength thermal imaging system for surface temperature monitoring and process control |
CN102798475A (zh) * | 2012-08-02 | 2012-11-28 | 电子科技大学 | 一种根据标准红外视频图像确定目标温度的方法 |
Non-Patent Citations (3)
Title |
---|
梁其诚等: "基于彩色火焰图像测量柴油机燃烧火焰温度", 《工业安全与环保》 * |
毛功平等: "生物柴油缸内燃烧火焰温度对NO_x生成的影响", 《内燃机工程》 * |
马凡华等: "一种求取柴油机燃烧火焰温度的新方法", 《内燃机工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846366A (zh) * | 2018-06-21 | 2018-11-20 | 江苏大学 | 计算内燃机缸内浓度场的方法及装置 |
CN109764820A (zh) * | 2018-12-24 | 2019-05-17 | 西华大学 | 一种定容燃烧火焰传播半径的测量角度步长确定方法 |
CN109764820B (zh) * | 2018-12-24 | 2020-08-11 | 西华大学 | 一种定容燃烧火焰传播半径的测量角度步长确定方法 |
CN114152347A (zh) * | 2021-09-30 | 2022-03-08 | 国网黑龙江省电力有限公司电力科学研究院 | 一种变电站电力设备故障定位及火险研判综合检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103344339B (zh) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111429540B (zh) | 温度场和变形场同步测量装置和方法 | |
CN102706529B (zh) | 超声速流场密度场的校准以及测量超声速密度场的方法 | |
Fuhrmann et al. | Two-dimensional cycle-resolved exhaust valve temperature measurements in an optically accessible internal combustion engine using thermographic phosphors | |
CN105606222A (zh) | 一种火焰三维温度场测量的成像装置、测量装置及测量方法 | |
CN107401977B (zh) | 高温双目立体视觉测量中考虑折射偏差的成像补偿方法 | |
CN203375949U (zh) | 高温适性红外观瞄***光轴平行性及mrtd野外检测仪 | |
CN103792090A (zh) | 一种柴油机缸内喷雾与燃烧可视化测试装置及测试方法 | |
CN203745204U (zh) | 一种柴油机缸内喷雾与燃烧动态可视化测试装置 | |
CN106404410A (zh) | 一种同时测量柴油喷雾结构及燃烧特性的装置及方法 | |
Shawal et al. | High-speed imaging of early flame growth in spark-ignited engines using different imaging systems via endoscopic and full optical access | |
CN103344339B (zh) | 一种基于内燃机缸内燃烧火焰摄影照片计算温度场的方法 | |
CN105973852B (zh) | 燃料射流浓度场分布测试装置及其实施方法 | |
CN101718620A (zh) | 多光谱动态调制传递函数测量方法与装置 | |
Goschütz et al. | Endoscopic imaging of early flame propagation in a near-production engine | |
Reggeti et al. | Two-color pyrometry system to eliminate optical errors for spatially resolved measurements in flames | |
Pastor et al. | Schlieren methodology for the analysis of transient diesel flame evolution | |
Lillo et al. | Plenoptic single-shot 3D imaging of in-cylinder fuel spray geometry | |
CN110376129B (zh) | 一种燃烧温度场及碳烟浓度场同步测量的方法及装置 | |
Payri et al. | Diesel spray analysis by means of planar laser-induced exciplex fluorescence | |
van Overbrüggen et al. | Holographic particle image velocimetry measurements in a four-valve combustion engine | |
Friedrich et al. | Mixture formation in a CNG-DI engine in stratified operation | |
Chen et al. | Compression and end-gas temperatures from iodine absorption spectra | |
Zhou et al. | Research on a quantitative method for three-dimensional computed tomography of chemiluminescence | |
Boxx et al. | Laser-induced incandescence measurements in a fired diesel engine at 3 kHz | |
Singh et al. | Tomographic particle image velocimetry for flow analysis in a single cylinder optical engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201123 Address after: No. 159, Chengjiang Middle Road, Jiangyin City, Wuxi City, Jiangsu Province Patentee after: Jiangyin Intellectual Property Operation Co., Ltd Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Patentee before: JIANGSU University |