CN108256166A - 一种用于热图技术的数据处理方法 - Google Patents
一种用于热图技术的数据处理方法 Download PDFInfo
- Publication number
- CN108256166A CN108256166A CN201711420400.8A CN201711420400A CN108256166A CN 108256166 A CN108256166 A CN 108256166A CN 201711420400 A CN201711420400 A CN 201711420400A CN 108256166 A CN108256166 A CN 108256166A
- Authority
- CN
- China
- Prior art keywords
- test
- model
- thermo
- model surface
- data processing
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/065—Measuring arrangements specially adapted for aerodynamic testing dealing with flow
- G01M9/067—Measuring arrangements specially adapted for aerodynamic testing dealing with flow visualisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明涉及一种用于热图技术的数据处理方法,通过来流建立前多幅图像算术平均,降低初始图像噪声,通过给定函数形式的多参数最小二乘优化算法降低试验过程中噪声,最终在降低CMOS相机试验结果噪声的同时,保证试验数据精度和热流变化细节。本发明采用CMOS相机,相对于以往的单点式传感器,能够获得相机视场范围内的全流场热流数据,数据更详细,不仅能够获得能够获得热流的定量数据,而且能够获得热流的分布趋势。
Description
技术领域:
本发明涉及一种用于热图技术的数据处理方法,尤其指一种应用于脉冲风洞中模型表面热流密度大面积非接触测量的数据处理方法,属于图像处理技术领域。
背景技术:
针对高超声速飞行器的热环境预测,传统的试验测量手段采用点式或整体式热流传感器,这种方式只能获得离散点的热流分布情况,而且容易破坏模型表面结构,对流场产生干扰;热图技术采用非接触测量方式,可以在不干扰流场的前提下获得大面积的热流数据,成为了热流测量技术发展的重要方向。
用于热图技术的相机主要分为两大类,CCD相机和CMOS相机;其中CCD相机量子效率高,噪声少,得到了广泛的运用。但是由于结构的原因CCD相机的采集速度很难做到高速,而CMOS图像传感器电荷读出速度更快,使得高速相机绝大多数采用CMOS感光芯片。随着高超声速飞行器的快速发展,对试验结果的精细化程度要求越来越高,而对于热图技术常用的炮/激波风洞,风洞运行时间只有20ms左右。而现有的CCD相机最快的采集帧速只有200fps,在风洞有效时间内,只能获取1-2张有效数据,已经逐渐不能满足试验的需求。CMOS图像传感器的主要问题就是噪声超过CCD芯片,这使得采用CMOS相机进行热图试验时,试验结果中混杂有大量的噪声信号,不能得到有效的试验结果。
发明内容:
为了克服CMOS相机在热图试验中噪声大的不足,本发明提供了一种用于热图技术的数据后处理方法,在降低结果噪声的同时,保证试验数据精度。
本发明的技术方案如下:
提供一种用于热图技术的数据处理方法,包括如下步骤:
(1)试验模型的表面喷涂温敏涂层;通过CMOS相机连续获得试验模型的荧光图像;
(2)获取试验前多幅荧光图像的光强分布,通过算术平均获得来流前模型表面的初始平均光强分布;
(3)开始试验,获取试验中多幅荧光图像的光强分布;
(4)将每幅试验中荧光图像的光强分布与初始平均光强分布进行对比,获得模型表面光强分布变化率;
(5)将每幅模型表面光强分布变化率与标定曲线对比,获得试验过程中的模型表面温度;
(6)利用试验过程中的模型表面温度获取试验过程中模型表面热流。
优选的,步骤(3)中还包括对CMOS相机获取的荧光图像进行中值滤波,获取的强度分布为滤波后的荧光图像光强分布。
优选的,步骤(6)中获取试验过程中模型表面热流的具体方法为:
qw(i)=A(B-θm(i))(i=1、2……m)
假设从试验开始到试验结束依次采集m张图像,设为m个采样点,步骤(5)中获得的第i个采样点处的表面温度数据为θm(i)(i=1、2……m);A为对流换热系数,B为绝热壁温,qw(i)为模型表面热流。
优选的,A、B的获取方式如下:
构造优化函数J为:
采用最小二乘优化的算法求解使得J最小的A、B值;其中ρ为模型材料密度,c为模型比热、k为模型导热系数,t为时间。
本发明相对于现有技术的优点在于:
(1)本发明采用CMOS相机,相对于以往的单点式传感器,能够获得相机视场范围内的全流场热流数据,数据更详细,不仅能够获得能够获得热流的定量数据,而且能够获得热流的分布趋势。
(2)本发明通过多幅荧光强度算术平均获得来流前模型表面的初始平均光强分布,降低初始帧的噪声水平,试验过程中对每幅图像进行中值滤波,降低试验过程中图像的空间噪声,对各个图像在时间序列上进行多变量最小二乘优化,降低图像的时间噪声。通过以上三个流程的去噪,最终降低了热流分布云图噪声,提高了热流计算精度。
(3)本发明通过采用优化算法,在去除试验结果噪声的同时保证了热流数据变化的真实细节。
附图说明:
图1是本发明数据处理方法流程。
具体实施方式:
本发明用于热图技术的数据处理方法,包括如下步骤:
(1)试验模型的表面喷射温敏涂层,在紫外光源的激发下辐射荧光,通过高速CMOS相机连续获得试验模型的荧光图像;
(2)获取试验前多幅荧光图像的光强分布,通过算术平均获得来流前模型表面的初始平均光强分布,通过算术平均,降低了试验前图像的随机噪声,保证了初始帧的精度;
(3)开始试验,试验过程中试验模型的温度升高,表面荧光强度下降,对高速CMOS相机获取的荧光图像进行中值滤波,降低试验过程中图像的空间噪声;获取试验中多幅滤波后的荧光图像的光强分布;
(4)将每幅试验中荧光图像的光强分布与初始平均光强分布进行对比,获得模型表面光强分布变化率;
(5)将每幅模型表面光强分布变化率与标定曲线对比,获得模型表面温度;通过试验中多幅荧光图像的强度分布变化率,获取试验过程中模型表面温度随时间的变化数据;所述标定曲线为试验模型的表面喷涂温敏涂层的材料光强变化率随温度的曲线,在试验前通过对温敏涂层的材料进行标定获得;
(6)获取试验过程中模型表面热流值随时间的变化数据。
热图技术热流计算基于一维半无限大假设。由于在脉冲设备中试验时间在20ms左右,试验时间非常短,可以认为在整个试验内热流为一恒定值,控制方程为:
飞行器周围空气和飞行器壁面间进行对流换热,由牛顿冷却公式可得:
qw=h(Taw-Tw) (2)
其中,qw模型表面热流,h为风洞来流和模型壁面间的对流换热系数,它与流体的物性,流动状态和模型表面情况等有关,T(x,t)为模型温度,x为该位置距离模型表面的距离,x=0表明为模型表面,t为试验时间,Taw为模型表面的绝热壁温,Tw为模型壁面温度等于T(0,t)。忽略试验过程中模型材料密度ρ(t)、比热c(t)、导热系数k(t)的变化,认为为常数,表示为ρ、比热c、导热系数k。将式(2)和一维导热微分方程组联立,可求出该边界条件下导热微分方程的解析解:
其中erfc为误差函数,其变量为t为试验时间,T0为模型初始温度,即试验开始时的温度。设A=h,B=Taw-T0,θ(x,t)=T(x,t)-T0,当x=0时,θw=θ(0,t)=T(0,t)-T0。则(3)式可化为:
该式可以看作关于A、B的多变量优化问题。假设从试验开始到试验结束依次采集m张图像,设为m个采样点,步骤(5)中获得的每个采样点处的表面温度数据为θm(i)(i=1、2……m),i=1为试验开始后采集的第一张图像,m为试验中采集的最后一张图像。计算时首先给A、B赋初值,A的初值为0至100之间的常数,B的初值为试验时流场的总温。
然后由式(4)可以计算出不同时刻的表面温度为θw(i)(i=1、2……m),此时构造优化函数J为:
这里采用最小二乘优化的算法,寻找最优的A、B值,使得目标函数J达到极小值,此时式(5)就是试验观测数据的最佳逼近函数。求出A、B以后由式(6)即可求出不同时刻该点的热流值为:
qw(i)=A(B-θm(i))(i=1、2……m) (6)
利用公式(6)获得模型表面的热流分布。
采用本发明的数据处理方法,热流分布噪声水平由12.5‰降低到4.3‰,大幅提高了测量精度。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (4)
1.一种用于热图技术的数据处理方法,其特征在于,包括如下步骤:
(1)试验模型的表面喷涂温敏涂层;通过CMOS相机连续获得试验模型的荧光图像;
(2)获取试验前多幅荧光图像的光强分布,通过算术平均获得来流前模型表面的初始平均光强分布;
(3)开始试验,获取试验中多幅荧光图像的光强分布;
(4)将每幅试验中荧光图像的光强分布与初始平均光强分布进行对比,获得模型表面光强分布变化率;
(5)将每幅模型表面光强分布变化率与标定曲线对比,获得试验过程中的模型表面温度;
(6)利用试验过程中的模型表面温度获取试验过程中模型表面热流。
2.如权利要求1所述的用于热图技术的数据处理方法,其特征在于,步骤(3)中还包括对CMOS相机获取的荧光图像进行中值滤波,获取的强度分布为滤波后的荧光图像光强分布。
3.如权利要求1或2所述的用于热图技术的数据处理方法,其特征在于,步骤(6)中获取试验过程中模型表面热流的具体方法为:
qw(i)=A(B-θm(i))(i=1、2……m)
假设从试验开始到试验结束依次采集m张图像,设为m个采样点,步骤(5)中获得的第i个采样点处的表面温度数据为θm(i)(i=1、2……m);A为对流换热系数,B为绝热壁温,qw(i)为模型表面热流。
4.如权利要求3所述的用于热图技术的数据处理方法,其特征在于,A、B的获取方式如下:
构造优化函数J为:
采用最小二乘优化的算法求解使得J最小的A、B值;其中ρ为模型材料密度,c为模型比热、k为模型导热系数,t为试验时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711420400.8A CN108256166A (zh) | 2017-12-25 | 2017-12-25 | 一种用于热图技术的数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711420400.8A CN108256166A (zh) | 2017-12-25 | 2017-12-25 | 一种用于热图技术的数据处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108256166A true CN108256166A (zh) | 2018-07-06 |
Family
ID=62723805
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711420400.8A Pending CN108256166A (zh) | 2017-12-25 | 2017-12-25 | 一种用于热图技术的数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108256166A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110057399A (zh) * | 2019-03-28 | 2019-07-26 | 东南大学 | 一种基于3d-dic的温度场与位移场同步测量系统及测量方法 |
CN112504611A (zh) * | 2020-12-29 | 2021-03-16 | 中国航天空气动力技术研究院 | 用于常规风洞的变姿态热图测量方法、系统及装置 |
CN112699583A (zh) * | 2020-12-28 | 2021-04-23 | 中国航天空气动力技术研究院 | 一种用于热图技术的数据空间分块处理方法 |
CN112834159A (zh) * | 2020-12-30 | 2021-05-25 | 中国航天空气动力技术研究院 | 一种翼舵缝隙内部及舵轴热流测量方法 |
CN112834158A (zh) * | 2020-12-30 | 2021-05-25 | 中国航天空气动力技术研究院 | 一种内通道热流面测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101594457A (zh) * | 2008-05-30 | 2009-12-02 | 深圳艾科创新微电子有限公司 | 一种分级降噪系统及方法 |
CN103616156A (zh) * | 2013-12-11 | 2014-03-05 | 中国航天空气动力技术研究院 | 脉冲风洞热流测量系统、方法及其标定系统、方法 |
CN106610316A (zh) * | 2016-12-29 | 2017-05-03 | 重庆工商大学 | 基于热波动耦合红外成像的薄壁局部换热系数测量方法 |
CN103776611B (zh) * | 2013-12-19 | 2017-07-28 | 中国航天空气动力技术研究院 | 脉冲风洞热流测量系统、方法及其标定系统、方法 |
-
2017
- 2017-12-25 CN CN201711420400.8A patent/CN108256166A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101594457A (zh) * | 2008-05-30 | 2009-12-02 | 深圳艾科创新微电子有限公司 | 一种分级降噪系统及方法 |
CN103616156A (zh) * | 2013-12-11 | 2014-03-05 | 中国航天空气动力技术研究院 | 脉冲风洞热流测量系统、方法及其标定系统、方法 |
CN103776611B (zh) * | 2013-12-19 | 2017-07-28 | 中国航天空气动力技术研究院 | 脉冲风洞热流测量系统、方法及其标定系统、方法 |
CN106610316A (zh) * | 2016-12-29 | 2017-05-03 | 重庆工商大学 | 基于热波动耦合红外成像的薄壁局部换热系数测量方法 |
Non-Patent Citations (3)
Title |
---|
毕志献 等: "磷光热图测热技术研究", 《实验流体力学》 * |
郭爱芳: "《传感器原理及应用》", 31 May 2007 * |
韩曙光 等: "磷光热图技术在常规高超声速风洞热环境实验中的应用", 《气体物理》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110057399A (zh) * | 2019-03-28 | 2019-07-26 | 东南大学 | 一种基于3d-dic的温度场与位移场同步测量系统及测量方法 |
CN112699583A (zh) * | 2020-12-28 | 2021-04-23 | 中国航天空气动力技术研究院 | 一种用于热图技术的数据空间分块处理方法 |
CN112699583B (zh) * | 2020-12-28 | 2023-06-20 | 中国航天空气动力技术研究院 | 一种用于热图技术的数据空间分块处理方法 |
CN112504611A (zh) * | 2020-12-29 | 2021-03-16 | 中国航天空气动力技术研究院 | 用于常规风洞的变姿态热图测量方法、系统及装置 |
CN112834159A (zh) * | 2020-12-30 | 2021-05-25 | 中国航天空气动力技术研究院 | 一种翼舵缝隙内部及舵轴热流测量方法 |
CN112834158A (zh) * | 2020-12-30 | 2021-05-25 | 中国航天空气动力技术研究院 | 一种内通道热流面测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108256166A (zh) | 一种用于热图技术的数据处理方法 | |
CN106501227B (zh) | 基于压力敏感涂料探针分子荧光寿命的测量方法 | |
CN102297660B (zh) | 一种盾构隧道衬砌管片接缝张开宽度的测量方法及装置 | |
Ni et al. | Lagrangian acceleration measurements in convective thermal turbulence | |
CN104807611B (zh) | 基于视频的烟气流速场、湍流度场实验测量装置及方法 | |
CN104976960B (zh) | 一种雨滴物理特性观测方法 | |
CN108225718A (zh) | 一种基于温度效应的荧光油膜厚度标定方法 | |
CN113588204B (zh) | 一种进气道激波边界层干扰特性测量方法 | |
CN102795627B (zh) | 多晶硅还原炉多参量在线监测与优化控制装置及方法 | |
CN103196564A (zh) | 一种通过图像分割修正表面发射率的红外热成像测温方法 | |
CN103698256A (zh) | 一种全场彩虹在线测量液体喷雾的方法和装置 | |
CN104613888B (zh) | 一种火焰烟雾环境下透过烟雾物体变形测量方法 | |
Tokgoz et al. | Temperature and velocity measurements in a fluid layer using background-oriented schlieren and PIV methods | |
Zhang et al. | Development of a CCD-based pyrometer for surface temperature measurement of casting billets | |
CN108050955B (zh) | 基于结构光投影与数字图像相关的高温空气扰动滤除方法 | |
CN105965017A (zh) | 一种slm成形过程的温度场监控装置及方法 | |
CN111189541B (zh) | 一种基于红外热像扫描的在线热损失测量方法 | |
CN109100022B (zh) | 测温方法和系统 | |
CN104748793B (zh) | 铝电解槽熔体温度与流速实时组合测量装置及其测量方法 | |
CN106885824A (zh) | 一种基于红外图像处理的建筑墙体传热系数检测方法 | |
He et al. | Velocity measurement of blast furnace molten iron based on local multi-feature correction using multi-stage filtered high-speed camera | |
CN207610670U (zh) | 一种液膜厚度测量装置 | |
Stelter et al. | Thermographic 3D particle tracking velocimetry for turbulent gas flows | |
CN209198266U (zh) | 颗粒聚光辐射吸收特性测量装置 | |
Petracci et al. | Combined PIV/LIF measurements in a Rayleigh-Bénard convection cell |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180706 |
|
RJ01 | Rejection of invention patent application after publication |