高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于数字全息的甲烷-氧气预混火焰温度场研究

刘柳 姚燕 蔡晋辉 梁晓瑜

引用本文:
Citation:

基于数字全息的甲烷-氧气预混火焰温度场研究

    作者简介: 刘柳(1997-), 女, 硕士研究生, 现主要从事光学信息处理和燃烧应用的研究.
    通讯作者: 姚燕, yaoyan@cjlu.edu.cn
  • 基金项目:

    国家重点研发计划资助项目 2018YFF01014002

  • 中图分类号: O438.1

Study on methane-oxygen premixed flame temperature field based on digital holography

    Corresponding author: YAO Yan, yaoyan@cjlu.edu.cn ;
  • CLC number: O438.1

  • 摘要: 为了获取相对稳定的温度场, 通过实验得到了不同比例下甲烷-氧气预混火焰温度场的变化情况, 并分析了甲烷质量分数对该预混火焰温度场的影响。借助Mckenna燃烧器, 采用数字全息技术, 获取不同甲烷-氧气比下预混火焰温度场的干涉条纹图像, 利用巴赫沃斯低通滤波降低经过预处理的图像的散斑噪声; 通过改进的四向最小二乘解包裹法提取温度场相位分布信息, 根据相位与温度的理论关系, 取得了相应的温度信息数据, 并利用B型热电偶进行了实验验证。结果表明, 利用该方法测得的温度与相同条件下热电偶测得的温度变化基本一致, 证明了数字全息技术测量温度场的可行性; 甲烷与氧气质量分数之比为0.9时, 其温度变化约为10K, 相对其它工况, 其温度场分布最为稳定。该研究为甲烷-氧气预混燃气的相关研究和Mckenna燃烧器的应用提供了理论基础。
  • Figure 1.  Schematic diagram of digital holography optical path

    a—schematic diagram of digital holography b—digital holographic physical map(1:He-Ne laser; 2:polarizer; 3, 7, 9:attenuator 1, 2, 3;4:beam expander; 5:diaphragm; 6, 12:spectroscope 1, 2;8, 10:retroreflector 1, 2;11:Mckenna burner; 13:CCD camera)

    Figure 2.  Interferogram of the temperature field at 15mm in the axial direction and 2mm from the center in the radial direction under working condition 1

    Figure 3.  Temperature field distribution change diagram of the measured area under the five working conditions (not unwrapped)

    a—condition 1 b—condition 2 c— condition 3 d—condition 4 e—condition 5

    Figure 4.  Temperature field distribution change diagram of the measured area under five working conditions (unwrapping 2-D)

    a—condition 1 b—condition 2 c—condition 3 d—condition 4 e—condition 5

    Figure 5.  Temperature field distribution change diagram of the measured area under five working conditions (unwrapping 3-D)

    a—condition 1 b—condition 2 c—condition 3 d—condition 4 e—condition 5

    Figure 6.  Line graph of average temperature of type B thermocouple under five working conditions

    Figure 7.  Box diagram of average temperature of type-B thermocouple under five working conditions

    Table 1.  Volume flow of each gas inlet under different working conditions

    CH4/
    (L·min-1)
    premixed air/
    (L·min-1)
    mass score ratio
    condition 1 1 11.90 0.8
    condition 2 1 10.58 0.9
    condition 3 1 9.52 1.0
    condition 4 1 8.66 1.1
    condition 5 1 7.94 1.2
    下载: 导出CSV
  • [1]

    RASTOGI A. Optical methods in heat transfer and fluid flow[J]. Optics and Lasers in Engineering, 2006, 44(3): 155-158.
    [2]

    HERRÁEZ J V, BELADA R. A study of free convection in air around horizontal cylinders of different diameters based on holographic interferometry. Temperature field equations and heat transfer coefficients[J]. International Journal of Thermal Sciences, 2002, 41(3): 261-267. doi: 10.1016/S1290-0729(01)01314-X
    [3]

    WANG H Y. Imaging theory and experimental demonstration of digital holographic microscopy[D]. Beijing: Beijing Industry University, 2008: 11-28(in Chinese).
    [4]

    TAKEDA M, INA H, KOBAYASHI S. Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry[J]. Review of Scientific Instruments, 1982, 72(1): 156-160. 
    [5]

    KERR D, F MENDOZASANTOYO, TYRER J R. Manipulation of the Fourier components of speckle fringe patterns as part of an interfe-rometric analysis process[J]. Optica Acta International Journal of Optics, 2017, 36(2): 195-203. 
    [6]

    YAMAGUCHI I. Fringe formation in speckle photography[J]. Journal of the Optical Society of America, 1984, A1(1): 81-86.
    [7]

    CHEN H, DENG Zh H, YU H M. Principle and application of thermocouple temperature measurement system[J]. Manufacturing Automation, 2004, 26(9): 68-71(in Chinese).
    [8]

    ZHANG L L. Study of temperature measurement applications based on the digital holographic interferometry and fluorescent material[D]. Ji'nan: Shandong Normal University, 2016: 19-28(in Chinese).
    [9]

    ZHU J R, DAI J M, YAO Y Ch, et al. Reconstruction of temperature field around horizontal cylinder based on double-exposure holography[J]. Laser Technology, 2015, 39(4): 541-544(in Chinese). 
    [10]

    QI J A, WONG W O, LEUNG C W, et al. Temperature field mea-surement of a premixed butane/air slot laminar flame jet with Mach-Zehnder Interferometry[J]. Applied Thermal Engineering, 2008, 28(14/15): 1806-1812.
    [11]

    DONG L L, CHEUNG C S, LEUNG C W. Heat transfer from an impinging premixed butane/air slot flame jet[J]. International Journal of Heat and Mass Transfer, 2002, 45(5): 979-992. doi: 10.1016/S0017-9310(01)00215-0
    [12]

    QI J A, LEUNG C W, WONG W O, et al. Temperature-field mea-surements of a premixed butane/air circular impinging-flame using reference-beam interferometry[J]. Applied Energy, 2006, 83(12): 1307-1316. doi: 10.1016/j.apenergy.2006.04.001
    [13]

    HAWORTH D C, BLINT R J, CUENOT B, et al. Numerical simulation of turbulent propane-air combustion with nonhomogeneous reactants[J]. Combustion & Flame, 2000, 121(3): 395-417.
    [14]

    WANG Y T. Numerical simulation on premixed combustion of methane-air under variable turbulence[D]. Beijing: Beijing Jiaotong University, 2020: 1-9 (in Chinese).
    [15]

    HERNANDEZ-MONTES M D S, MENDOZA-SANTOYO F, MORENO M F, et al. Macro to nano specimen measurements using photons and electrons with digital holographic interferometry: A review[J]. Journal of the European Optical Society Rapid Publications, 2020, 16(1): 16-47. doi: 10.1186/s41476-020-00133-8
    [16]

    LIU H, XIAO Y L, TIAN J L, et al. Nonlinear reconstruction for off-axis Fresnel digital holography with deep learning[J]. Acta Photonica Sinica, 2020, 49(7): 709001(in Chinese). doi: 10.3788/gzxb20204907.0709001
    [17]

    KATZ J, SHENG J. Applications of holography in fluid mechanics and particle dynamics[J]. Annual Review of Fluid Mechanics, 2010, 42(1): 531-555. doi: 10.1146/annurev-fluid-121108-145508
    [18]

    JIA J Ch, MO W D, LI Zh X, et al. Measurement of temperature field based on digital holography technology[J]. Progress in Laser and Optoelectronics, 2012, 49(3): 030903(in Chinese). doi: 10.3788/LOP49.030903
    [19]

    VANDERWEGE B A, O'BRIEN C J, HOCHGREB S. Quantitative shearography in axisymmetric gas temperature measurements[J]. Optics & Lasers in Engineering, 1999, 31(1): 21-39.
    [20]

    REUSS D L. Temperature measurements in a radially symmetric flame using holographic interferometry[J]. Combustion and Flame, 1983, 49(1/3): 207-219.
    [21]

    STELLA A, GU J G, GIAMMARTINI S. Measurement of axisymmetric temperature fields using reference beam and shearing interferometry for application to flames[J]. Experiments in Fluids, 2000, 29(1): 1-12. doi: 10.1007/s003480050420
    [22]

    YANG X, CHEN B, LI X Y, et al. Design and experimental analysis of the long distance digital off-axis holography[J]. Laser Journal, 2016, 37(10): 74-77(in Chinese).
    [23]

    DING W, KONG Y, DU T Y, et al. Research on digital holographic imaging denoising based on time-space domain depolarization[J]. Laser Technology, 2021, 45(1): 115-120(in Chinese).
    [24]

    KANG Q, WANG L, XU Y, et al. A new approach for image de-noising algorithm based on bayesian estimation[C]// Proceedings of 2016 3rd International Conference on Information and Communication Technology for Education. Cook, USA: American Applied Sciences Research Institute, 2016: 4.
    [25]

    ARAI Y. Improvement of measuring accuracy of spatial fringe analysis method using only two speckle patterns in electronic speckle pattern interferometry[J]. Optical Engineering, 2014, 53(3): 034107. doi: 10.1117/1.OE.53.3.034107
    [26]

    GUO Y, WU Q, M Q, et al. Four directions least-square phase unwrapping algorithm based on lateral shearing interferometry[J]. Journal of Optoelectronics·Laser, 2015, 26(10): 1953-1959(in Chinese). 
    [27]

    ZHU T, WANG Zh Y, XIE M, et al. Research of four directions weighted least-square algorithm for phase unwrapping[J]. Journal of Sichuan University(Natural Science Edition), 2009, 46(2): 372-376(in Chinese). 
    [28]

    PENG J, ZHANG R Zh. Comparison on the noise processing ability of common unwrapping algorithms in interference detection[J]. Optics and Optoelectronic Technology, 2019, 17(4): 77-83(in Chinese).
    [29]

    ZHANG X, ZHANG X, XU M, et al. Phase unwrapping in digital holography based on non-subsampled contourlet transform[J]. Optics Communications, 2018, 407: 367-374. doi: 10.1016/j.optcom.2017.09.057
    [30]

    CAI Z, LIU X, PEDRINI G, et al. Structured-light-field 3-D imaging without phase unwrapping-ScienceDirect[J]. Optics and Lasers in Engineering, 2020, 129(3): 106047.
  • [1] 陈希慧焦春妍李俊昌 . 空间载波相移法用于全息CT测量气体温度场. 激光技术, 2006, 30(4): 412-414.
    [2] 朱进容代金梅姚育成成纯富 . 双曝光全息术用于水平圆管温度场的重建. 激光技术, 2015, 39(4): 541-544. doi: 10.7510/jgjs.issn.1001-3806.2015.04.025
    [3] 杨富黄伟张彬楚晓亮刘志国蔡邦维 . CO2激光辐照下光学薄膜的温度场与热畸变. 激光技术, 2004, 28(3): 255-258.
    [4] 王占亮张永安钱晓凡 . 数字全息和拉东变换重构温度场的研究. 激光技术, 2011, 35(5): 589-592. doi: 10.3969/j.issn.1001-3806.2011.05.004
    [5] 魏亚运周昕王梦婷呼有军 . 基于光学扫描全息测量相位物体的相位分布. 激光技术, 2016, 40(6): 820-824. doi: 10.7510/jgjs.issn.1001-3806.2016.06.010
    [6] 盛兆玄孙新利 . 计算全息技术在光学相关识别中的应用. 激光技术, 2008, 32(1): 61-63.
    [7] 徐青曹娜曹亮韩长材雷岚王清华 . 静电雾化场的数字全息实验研究. 激光技术, 2013, 37(2): 143-146. doi: 10.7510/jgjs.issn.1001-3806.2013.02.002
    [8] 姜雨甫周昕王运郭雪帆 . 采用二值化相位光瞳的光学扫描全息方法. 激光技术, 2013, 37(6): 723-727. doi: 10.7510/jgjs.issn.1001-3806.2013.06.004
    [9] 乐建明周昕孙阿茜张鹏 . 改善光学扫描全息聚焦层图像效果的实验研究. 激光技术, 2017, 41(3): 332-336. doi: 10.7510/jgjs.issn.1001-3806.2017.03.006
    [10] 李雪梅朱春英马友光 . 激光全息干涉法测量维生素B6的扩散系数. 激光技术, 2008, 32(3): 281-283.
    [11] 任真樊则宾李俊昌邴小林李兴华 . 精密测量微小转角的数字全息方法及应用. 激光技术, 2012, 36(6): 798-801. doi: 10.3969/j.issn.1001-3806.2012.06.021
    [12] 易亨瑜彭勇田小强胡晓阳吕百达 . 光腔衰荡法的单波长反射率测量实验研究. 激光技术, 2005, 29(4): 337-339.
    [13] 谭中奇龙兴武 . 用于连续波腔衰荡法测量的LD驱动电路设计. 激光技术, 2008, 32(1): 27-29,56.
    [14] 刘仲寿万雄王庆危荃 . 基于DSP的温度场重建. 激光技术, 2009, 33(1): 53-56.
    [15] 胡增荣周建忠郭华锋杜建钧 . 应用ABAQUS模拟激光焊接温度场. 激光技术, 2007, 31(3): 326-329.
    [16] 曾大文谢长生 . 复合涂层激光熔池温度场及流场的数值模拟. 激光技术, 2000, 24(6): 370-374.
    [17] 张君魏晓峰张雄军吴登生田晓琳 . 重频应用下电光晶体温度场应变场影响因素. 激光技术, 2009, 33(1): 27-31.
    [18] 李明海柳爱国宋耀祖 . 激光放大介质温度场和热应力场的数值模拟. 激光技术, 2002, 26(2): 86-89.
    [19] 王秀凤吕晓东陈光南胡世光 . 激光强化温度场的数值模拟与校验. 激光技术, 2004, 28(2): 162-165.
    [20] 张建宇高立新崔玲丽吴迪平杨久霞王会刚 . 激光强化温度场的理论解析与实验论证. 激光技术, 2006, 30(1): 56-59.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  3051
  • HTML全文浏览量:  2182
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-22
  • 录用日期:  2021-05-04
  • 刊出日期:  2022-05-25

基于数字全息的甲烷-氧气预混火焰温度场研究

    通讯作者: 姚燕, yaoyan@cjlu.edu.cn
    作者简介: 刘柳(1997-), 女, 硕士研究生, 现主要从事光学信息处理和燃烧应用的研究
  • 中国计量大学 计量测试工程学院, 杭州 310018
基金项目:  国家重点研发计划资助项目 2018YFF01014002

摘要: 为了获取相对稳定的温度场, 通过实验得到了不同比例下甲烷-氧气预混火焰温度场的变化情况, 并分析了甲烷质量分数对该预混火焰温度场的影响。借助Mckenna燃烧器, 采用数字全息技术, 获取不同甲烷-氧气比下预混火焰温度场的干涉条纹图像, 利用巴赫沃斯低通滤波降低经过预处理的图像的散斑噪声; 通过改进的四向最小二乘解包裹法提取温度场相位分布信息, 根据相位与温度的理论关系, 取得了相应的温度信息数据, 并利用B型热电偶进行了实验验证。结果表明, 利用该方法测得的温度与相同条件下热电偶测得的温度变化基本一致, 证明了数字全息技术测量温度场的可行性; 甲烷与氧气质量分数之比为0.9时, 其温度变化约为10K, 相对其它工况, 其温度场分布最为稳定。该研究为甲烷-氧气预混燃气的相关研究和Mckenna燃烧器的应用提供了理论基础。

English Abstract

    • 洞悉温度场的时空分布特性是燃烧学领域中的重要研究内容之一。数字全息技术作为近年来迅速成长起来的一个新颖的科学计量手段,可实现可视化整个区域的实时流体运动和温度分布[1-2],有着实时监测、非接触测量等[3]优点,它还能存储数据,通过测量对象的幅度和相位,得到折射率、表面变形、表面形状、振动条件等[4-6]信息,并且能全程迅速记录瞬时温度场。传统的火焰温度传感器,如热电偶[7]在测温时由于其侵入式测量的特点,会对火焰产生干扰,且传统手段无法提供火焰温度的立体分布。对此,国内外学者开展了开拓性的研究。ZHANG[8]和ZHU等人[9]利用数字全息干涉技术成功实现对轴对称型温度场的测量。HERRÁEZ等人[2]基于全息干涉法探讨了表面温度和不同直径对水平加热筒周围空气对流的影响,并通过实验定义相应的温度场方程和传热系数。全息干涉术可通过可视化流体中的密度场,进而深入了解其温度分布。QI等人[10]利用马赫-曾德尔干涉法测定不同等效比率和不同雷诺数的预混丁烷/空气槽层火焰射流的温度场。根据现有的研究结果不难得出: 数字全息技术利用数字再现,通过编程手段消除在记录过程中引入的噪声等不利因素的影响,能完整快速地再现温度场,大大提高了再现质量; 且其测量过程无滞后、无干扰,节省了大量处理数据的时间和空间,显著提高了测量精度。

      由于利用预混燃气(如丁烷-空气、丙烷-空气)的燃烧器,能在低雷诺数和低压下完全燃烧和快速传热,因此,利用预混燃气的燃烧器被广泛用于轻型工业领域[11-13]。Mckenna燃烧器作为能模拟实际电厂燃烧环境的新型装置,是实验室燃烧实验的首选燃烧器。以甲烷-氧气为主的预混气是目前应用较为广泛的气体燃料[14]。在国内外研究中,利用数字全息技术探索不同甲烷-氧气比对甲烷-氧气预混火焰温度稳定性的研究尚且不足。因此,本文中借助Mckenna燃烧器,利用数字全息法得到不同甲烷-氧气比下的预混湍流火焰的干涉条纹图,并借助巴赫沃斯低通滤波、改进的四向最小二乘法等手段获取原始干涉条纹对应的2维相位分布图。基于温度与相位的关系,通过实验得到的2维和3维的相位分布,研究不同甲烷-氧气比下,其预混火焰在高度15mm处的温度分布,并用B型热电偶进行验证,成功找到使得预混火焰温度分布稳定的甲烷-氧气比。本实验对研究甲烷-氧气预混湍流火焰的燃烧温度场的测量和Mckenna燃烧器的应用具有实际指导借鉴意义。

    • 数字全息技术以传统光学全息为理论基础,其记录光路与传统光学全息基本相同。不同的是,传统光学全息以全息干板为介质记录干涉图,利用光学手段(如曝光、显影、定影等)对干涉图进行处理与再现,其过程繁琐、对环境要求较高,难以做到实时记录与再现。而数字全息以CCD等光学探测器件作为记录干涉图的介质,将干涉图直接输入计算机进行数字处理与再现。数字全息较传统光学全息最为优越的一点是,它能连续记录变化物场的多幅全息图,利用定量得到的被记录物体再现像的振幅和相位信息,实现对变化物场的实时测量。此处物场为甲烷-氧气预混火焰的燃烧场,通过处理燃烧场的相位即可得到其温度场变化情况。

      光是一种电磁波,在介质中传播时,其相位周期变化。当光束向前一个波长时,相位变化为2π,由此可得,对某束光而言,其光程与相位的对应关系为:

      $ \varphi=\frac{2 {\rm{ \mathsf{ π} }}}{\lambda} \int_{L} n \mathrm{d} l $

      (1)

      式中,l为理论光程距离,L为实际光程距离,λ为光波长,n为光经过的场的气体折射率。当波长相同的两束光发生干涉时,由(1)式可得到相位差:

      $ \Delta \varphi=\varphi_{1}-\varphi_{2}=\frac{2 {\rm{ \mathsf{ π} }}}{\lambda}\left(\int_{L_{1}} n_{1} \mathrm{d} l-\int_{L_{2}} n_{2} \mathrm{d} l\right) $

      (2)

      式中,φ1φ2分别为物光和参考光的相位,n1n2分别为物光和参考光所经过的透明场的气体折射率,L1L2分别为物光和参考光的传播距离,如图 1所示,L1=L2=L1′+L2′。该相位差被记录在通过CCD得到的干涉图中。

      Figure 1.  Schematic diagram of digital holography optical path

      物光和参考光在CCD感光面上叠加,得到干涉条纹的光强[15-17]为:

      $ \begin{array}{l} \ \ \ \ \ \ \ \ E\left(x_{H}, y_{H}\right)=\left|E_{\mathrm{r}}\left(x_{H}, y_{H}\right)\right|^{2}+\left|E_{\mathrm{o}}\left(x_{H}, y_{H}\right)\right|^{2}+ \\ E_{\mathrm{r}}\left(x_{H}, y_{H}\right) E_{\mathrm{o}}{ }^{*}\left(x_{H}, y_{H}\right)+E_{\mathrm{r}}{ }^{*}\left(x_{H}, y_{H}\right) E_{\mathrm{o}}\left(x_{H}, y_{H}\right) \end{array} $

      (3)

      式中,Er(xH, yH)为参考光复振幅,Eo(xH, yH)为物光复振幅,*表示复共轭振幅,xHyH是图像全息平面H上的xy坐标,(3)式的最后两项包含与使用快速傅里叶变换方法检索物光的振幅和相位所相应的信息。用振幅和相位表示物光和参考光:

      $ \left\{\begin{array}{l} E_{\mathrm{o}}\left(x_{H}, y_{H}\right)=e_{\mathrm{o}}\left(x_{H}, y_{H}\right) \mathrm{e}^{\mathrm{i} \varphi\left(x_{H}, y_{H}\right)} \\ E_{\mathrm{r}}\left(x_{H}, y_{H}\right)=e_{\mathrm{r}}\left(x_{H}, y_{H}\right) \mathrm{e}^{-2 {\rm{ \mathsf{ π} }} \mathrm{i}\left(f_{x} x_{H}+f_{y} y_{H}\right)} \end{array}\right. $

      (4)

      式中,er(xH, yH)和eo(xH, yH)分别为参考光和物光记录在CCD上的振幅,φ(xH, yH)表示物光记录在CCD上的相位分布,fxfyx, y方向的光波频率,fxxHfyyH分别表示沿xHyH方向的空间载波频率。

      将(4)式代入(3)式得:

      $ \begin{aligned} E\left(x_{H}, y_{H}\right)=& \alpha\left(x_{H}, y_{H}\right)+c\left(x_{H}, y_{H}\right) \mathrm{e}^{2 {\rm{ \mathsf{ π} }} \mathrm{i}\left(f_{x} x_{H}+f{_{y}{y}_{H}}\right)}+\\ & c^{*}\left(x_{H}, y_{H}\right) \mathrm{e}^{-2 {\rm{ \mathsf{ π} }} \mathrm{i}\left(f_{x} x_{H}+f_{y} y_{H}\right)} \end{aligned} $

      (5)

      其中,

      $ \left\{\begin{array}{l} \alpha\left(x_{H}, y_{H}\right)=e_{\mathrm{o}}{}^{2}\left(x_{H}, y_{H}\right)+e_{\mathrm{r}}{}^{2}\left(x_{H}, y_{H}\right) \\ c\left(x_{H}, y_{H}\right)=e_{\mathrm{o}}\left(x_{H}, y_{H}\right) e_{\mathrm{r}}\left(x_{H}, y_{H}\right) \mathrm{e}^{\mathrm{i} \varphi\left(x_{H}, y_{H}\right)} \end{array}\right. $

      (6)

      式中,α(xH, yH)表示物光和参考光的强度分布,c(xH, yH)表示干涉干扰项,c*c的复共轭振幅。(5)式的后面两项表示由于干扰信号而引入项φ(xH, yH)对空间载波进行相位调制。为了获取光学相位,对(5)式进行2维傅里叶变换得:

      $ \begin{gathered} \operatorname{FT}\left\{E\left(x_{H}, y_{H}\right)\right\}=I\left(f_{x_{H}}, f_{y_{H}}\right)+C\left(f_{x_{H}}-f_{x} x_{H}, f_{y_{H}}-\right. \\ \left.f_{y} y_{H}\right)+C^{*}\left(f_{x_{H}}+f_{x} x_{H}, f_{y_{H}}+f_{y} y_{H}\right) \end{gathered} $

      (7)

      式中,fxH, fyH为光波在图像全息平面H上的频率,I(fxH, fyH)是频率坐标为(fxH, fyH)的低频背景照明;CC*分别表示cc*的傅里叶变换,且每个都包含相同的相位信息φ(xH, yH),由于它们在傅里叶频谱中空间分离,因此可使用滤波器消除IC*项。C项以频率为中心,且对每个滤波后的全息图进行傅里叶逆变换操作。接着可以用:

      $ \varphi\left(x_{H}, y_{H}\right)=\arctan \left[\frac{\operatorname{Im} c\left(x_{H}, y_{H}\right)}{\operatorname{Re} c\left(x_{H}, y_{H}\right)}\right] $

      (8)

      式中,Im c(xH, yH)和Re c(xH, yH)是逆傅里叶全息图C项的虚部和实部,变化前后两个相位相减得到包裹相位图,其相对相位差为:

      $ \Delta \varphi\left(x_{H}, y_{H}\right)=\varphi\left(x_{H}, y_{H}\right)-\varphi^{\prime}\left(x_{H}, y_{H}\right) $

      (9)

      对得到的包裹相位图进行相位展开,即相位解包裹,就能从图中提取相位变化,原理如(9)式所示。由于相位差与折射率满足[7, 18]

      $ \Delta \varphi=2 {\rm{ \mathsf{ π} }} d(n-1) / \lambda $

      (10)

      式中,d为有效温度场的长度。通过(10)式可得温度场内折射率n的变化。

      通过VANDERWEGE,REUSS和STELLA等人[19-21]的研究发现,将甲烷-氧气预混火焰的温度场视为轴对称且介质透明均匀,对结果不会产生重大误差。因此,此处可根据Lorenz-Lorentz关系式和Gladstone-Dale公式得:

      $ 1 / T_{\mathrm{r}}-1 / T_{0}=R\left(n_{\mathrm{r}}-1\right) /(p K M) $

      (11)

      式中,Trnr为某点的温度和气体折射率,T0为初始温度;p=1.01×105Pa为气体压强;M=29g/mol为空气摩尔质量;K=2.43×10-4m3/kg,是Gladstone-Dale常数;R=8.31J/mol/K,是气体普适常量。通过(11)式进一步得到温度Tr

      结合(10)式和(11)式不难看出,相位变化与温度变化存在一致性,因此通过观察相位变化即可得到温度变化。

    • 本实验中采用的马赫-曾德尔干涉光路,是一种常用的数字全息光路,如图 1所示。以He-Ne激光器(波长λ=632.8nm)为光源,光束通过偏振器、衰减器和扩束准直器后,形成一束平行光,调节光阑,得到合适直径大小(约3mm)的光束,使其通过可调分光镜1后,分成物光和参考光。其中,在CCD相机前的分光镜2处,参考光和载有燃烧温度场信息的物光汇合。经过分光镜的两束光以一定夹角在CCD上形成干涉条纹图,并存储于PC端。实验对象为带中心管的Mckenna燃烧器,甲烷和空气按一定流量配比混合燃烧以获得稳定的平面火焰,为后续的生物质样本燃烧提供高温环境。

      所用CCD感光面尺寸为512mm×512mm,分辨率为2016pixel×2016pixel,像素尺寸为3.1μm。拍摄全息时,物参光夹角的选取需要考虑两个条件[22]:(1)全息图的±1级和0级频谱可以在空间分离;(2)记录在CCD上的微干涉条纹的空间频率fmax必须大于CCD截止频率1/(2Δx),Δx为CCD的像素尺寸,即每一个条纹周期内最少需要两个采样点。满足上述条件后,记录在CCD上的物参光夹角满足θλ/(2Δx)。实验中采用He-Ne激光器,波长为λ=632.8nm,所用CCD的像素尺寸为3.1μm。由于所用光束直径较小,仅为3mm,而燃烧器中心送样空气出口直径为6mm,因此本实验中得到的全息干涉图是燃烧温度场的一部分。为方便衍射计算,在实际计算中,使用MATLAB软件定位截取原始全息图(2016pixel×2016pixel)的中间部分(256pixel×256pixel)进行数据处理。

    • 实验室燃烧环境一般分为3类,分别是富氧环境(甲烷-氧气质量分数之比小于1.0)、完全反应环境(甲烷-氧气质量分数比为1.0)和富燃料环境(甲烷-氧气质量分数比大于1.0)。本次实验中设计5种不同的工况,5种工况的各气体进口体积流量见表 1

      Table 1.  Volume flow of each gas inlet under different working conditions

      CH4/
      (L·min-1)
      premixed air/
      (L·min-1)
      mass score ratio
      condition 1 1 11.90 0.8
      condition 2 1 10.58 0.9
      condition 3 1 9.52 1.0
      condition 4 1 8.66 1.1
      condition 5 1 7.94 1.2

      实验室测得燃烧环境的冷态温度为293K,以轴向15mm、径向距中心2mm处的燃烧温度场为例,得到未通入甲烷燃烧时(即初始状态下的全息图)和燃烧器在工况1、工况2、工况3、工况4和工况5下的6幅全息干涉图。由于6幅干涉图显示相似,此处以工况1为例说明,将MATLAB软件定位截取256pixel×256pixel的全息图放大观察,如图 2所示,实际数据处理时不进行放大。

      Figure 2.  Interferogram of the temperature field at 15mm in the axial direction and 2mm from the center in the radial direction under working condition 1

      图 2干涉图明暗条纹相间,反映了物体周围温度随燃烧工况的变化,温度相同处光程相同,折射率也相同,即每一干涉条纹就是一条等温线。由于温度变化前后物光波的相位分布直接决定干涉条纹的分布,利用(11)式就可以精确计算出物场中每一像素所对应得空间位置处得温度变化。

    • 在记录数字全息干涉图时,会引入多种噪声,如散斑噪声、环境噪声和拍摄噪声等。其中,尽可能减小散斑噪声的影响成为全息技术的一个研究重点[23-25]。高度相干光源产生的光容易和其它光发生干涉,导致形成散斑噪声,可表示为:

      $ g(a, b)=f(a, b) \cdot u(a, b)+v(a, b) $

      (12)

      式中,f(a, b)为原始全息图在坐标(a, b)处的数值,g(a, b)为原始全息图像受到散斑噪声污染后的图像,u(a, b),v(a, b)分别是和原始图像分布相互独立的乘性噪声分量和加性噪声分量。本文中利用傅里叶变换、巴赫沃斯低通滤波和傅里叶逆变换等手段降低散斑噪声,即采用各种图像处理的算法对全息图进行处理,并通过软件处理图像达到降噪的目的。将未燃烧状态下的干涉图分别和工况1~工况5的干涉图做数字处理得到相位被压缩在(-π, +π)区间内的包裹相位图,如图 3a~图 3e所示,接着分别将其解包裹相位[26-30]分别得到图 4图 5,图中色度条表示相位变化,单位为rad。

      Figure 3.  Temperature field distribution change diagram of the measured area under the five working conditions (not unwrapped)

      Figure 4.  Temperature field distribution change diagram of the measured area under five working conditions (unwrapping 2-D)

      Figure 5.  Temperature field distribution change diagram of the measured area under five working conditions (unwrapping 3-D)

      图 3图 4中得到的被测空气区域的相位变化Δφ表示未通入甲烷燃烧时与燃烧器分别在工况1~工况5情况下该区域温度变化的情况,其中,横、纵坐标分别表示数字矩阵的行列数。由于燃烧器的燃气是从面板上的多个小孔输出,即会产生无数个小火焰从而形成燃烧器整个平面的大火焰,因此图 3中的每一个封闭等温环均表示为一个小火焰。图 3中条纹密集处表示温度梯度大,温度变化剧烈,反之亦然。火焰干涉条纹的明亮区域被认为基本不受火焰的影响,即温度不变。由图 4可以看出,相比较而言,工况1、工况2和工况4下温度干涉图中的明亮区域较多,且工况1、工况2和工况4下的等温环分布较为均匀,因此较工况3和工况5而言,工况1、工况2和工况4的燃烧温度较为稳定。

      图 5中可以看到,工况1(即甲烷-氧气质量分数之比为0.8)时,被测区域温度场的相位差约为1.6rad;工况2(即甲烷-氧气质量分数之比为0.9)时,被测区域温度场的相位差约为1.2rad;工况3(即甲烷-氧气质量分数之比为1.0)时,被测区域温度场的相位差约为3.0rad;工况4(即甲烷-氧气质量分数之比为1.1)时,被测区域温度场的相位差约为3.0rad;工况3(即甲烷-氧气质量分数之比为1.2)时,被测区域温度场的相位差约为2.5rad。由于相位变化幅度体现出温度变化趋势,两者变化一致,因此可以看出,随着甲烷-氧气比的增加,即甲烷含量变多,甲烷-氧气预混火焰在15mm高度处的温度场变化情况是先变小后变大。不难发现,工况2温度场的最大相位与最小相位之差最小,即甲烷-氧气比为0.9时,该区域温度变化最平缓。

      将3支B型热电偶以120°夹角分布在燃烧器周围,采用温度记录仪记录不同工况下燃烧器出口15mm高度处的温度值,每隔0.005s记录一次,在45s内记录9000次,并对3支热电偶测得数据平均处理,得到温度变化曲线图,由于曲线数值点过多,此处采用拟合平滑处理,如图 6所示,更有利于查看温度曲线变化趋势。

      Figure 6.  Line graph of average temperature of type B thermocouple under five working conditions

      由于热电偶从冷态突然进入燃烧环境,热电偶在0s~10s左右时间内测得的温度起伏大,如图 6所示,无实际参考意义,因此对10s后的温度曲线进行观察。由于温度点过于密集,10s后的温度变化采用箱线图表示,如图 7所示。从图 7中不难发现,相比其它工况而言,工况2在10s后的温度变化范围最小,约为10K,验证了从图 5得出的工况2的温度变化最为平缓的结论。从实验结果中发现,随着甲烷含量变多,甲烷-氧气预混火焰燃烧情况并不是一直趋向平稳,甲烷含量过少或过多,都会导致甲烷-氧气预混火焰的温度变化激烈。

      Figure 7.  Box diagram of average temperature of type-B thermocouple under five working conditions

    • 采用数字全息技术,得到不同甲烷-氧气比的预混火焰的干涉条纹图,并对所得条纹图进行傅里叶变换,获取对应的频域图像。接着利用巴赫沃斯低通滤波对频域图像进行滤波和傅里叶逆变换处理,获得相应包裹相位。最后借助改进的四向最小二乘法等手段获取原始干涉条纹对应的2维和3维的相位图。

      本实验中借助Mckenna燃烧器,通过2维和3维的相位图研究不同比例下甲烷-氧气预混火焰的空间温度分布,得到了在工况2(即甲烷-氧气质量分数之比为0.9)条件下,15mm高度处的甲烷-氧气预混火焰温度场最为稳定,约为10K的结论,并用B型热电偶验证该结论的准确性。实验成功验证了数字全息技术对于检测温度场等流场分布测量在技术上的优越性和灵敏度较高的特点。

参考文献 (30)

目录

    /

    返回文章
    返回