CN107843347B - 一种多孔介质三维温度分布测量方法 - Google Patents
一种多孔介质三维温度分布测量方法 Download PDFInfo
- Publication number
- CN107843347B CN107843347B CN201711100574.6A CN201711100574A CN107843347B CN 107843347 B CN107843347 B CN 107843347B CN 201711100574 A CN201711100574 A CN 201711100574A CN 107843347 B CN107843347 B CN 107843347B
- Authority
- CN
- China
- Prior art keywords
- porous media
- temperature
- heat
- thermal conductivity
- fluid
- 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
- 238000000691 measurement method Methods 0.000 title claims abstract description 6
- 238000009825 accumulation Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000012546 transfer Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 239000000523 sample Substances 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 12
- 239000011347 resin Substances 0.000 claims description 10
- 229920005989 resin Polymers 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 9
- 230000001052 transient effect Effects 0.000 abstract description 6
- 238000009529 body temperature measurement Methods 0.000 abstract description 4
- 238000001931 thermography Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 3
- 239000011521 glass Substances 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 102000016938 Catalase Human genes 0.000 description 1
- 108010053835 Catalase Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000003729 cation exchange resin Substances 0.000 description 1
- 229940023913 cation exchange resins Drugs 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000011257 shell material Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012956 testing procedure Methods 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提出一种多孔介质三维温度分布测量方法,利用瞬态导热系数测量装置在线测量堆积多孔介质的有效导热系数等热物理参数和热中心温度,同时利用非接触测温装置和技术测量堆积多孔介质的外表面温度分布;最后根据建立的导热微分方程组和单值性条件求解得到三维温度场。本发明采用瞬态法在线测量有效导热系数,不受化学反应干扰,可在线完成;采用非接触测温,可降低和消除非接触测温干扰;应用传热学知识得到的导热微分方程组适用在不同领域和物理背景,计算误差较小。本发明可以广泛应用在涉及微小尺度多孔介质的温度测量以及传热传质研究。
Description
技术领域
本发明涉及一种多孔介质三维温度分布测量方法,特别是一种用于微小尺度堆积多孔介质的三维温度分布测量方法。
背景技术
一直以来,多孔介质的传热传质研究是传热传质领域关注的焦点和难点,目前,尤其较少文献涉及多孔介质的温度测量,并且微小尺度极大地提升了难度,相对于一般尺度而言,微小尺度极大地提升了难度,接触式测温会产生相对较大的干扰和误差,并且容易移动位置,同时非接触测温不透光,核磁、层析和磁扫描测温技术远远还不成熟,并且成本很高,故对于微小尺度多孔介质三维温度场的测量困难重重。
发明内容
本发明的目的是解决现有技术的缺陷,提供一种适用于微小尺度堆积多孔介质三维温度分布测量方法。
本发明的技术方案是要测量微小尺度堆积多孔介质的三维温度分布,包括如下步骤,首先在堆积多孔介质的热中心(几何中心)布置瞬态导热系数测量装置,在堆积多孔介质的四周布置非接触测温装置;然后利用瞬态导热系数测量装置在线测量堆积多孔介质的有效导热系数等热物理参数和堆积多孔介质的热中心温度,同时利用非接触测温装置和技术测量堆积多孔介质的外表面温度分布;最后,基于上述步骤得到的有效导热系数、热中心温度等热物理参数和外表面温度分布,根据建立的导热微分方程组和单值性条件求解得到堆积多孔介质的三维温度场。
所述瞬态导热系数测量装置可以是探针法和热线法等非稳态测量导热系数的装置。
所述非接触测温装置可以是红外热像仪等非接触方法测量温度分布的装置。
所述导热微分方程组和单值性条件是依据传热学理论,根据堆积多孔介质的几何模型和物理模型建立的,并可通过数学解析方法和数值方法求解。
本发明的有益效果是,本发明采用瞬态法在线测量多孔介质的有效导热系数,不受化学反应干扰,在线完成测试,不受多孔介质堆积状态影响;采用非接触测温,提高测温精度,无需接触式布点测温便可得到温度分布,可降低和消除非接触测温干扰,不受微小尺度限制,测温精度高;应用传热学知识得到的导热微分方程组可适用在不同领域和物理背景的多孔介质,数学求解可以到解析解,即便复杂问题也可得到数值解,方法成熟可靠,计算误差较小。本发明可以广泛应用在涉及微小尺度多孔介质的温度测量以及传热传质研究。
附图说明
附图1是本发明实施例的原理示意图。
附图2 是本发明实施例得到的不同时刻的热像图。
附图3是本发明实施例得到的典型时刻的表面温度分布。
附图4是本发明实施例得到的典型时刻的三维温度分布。
其中,1、探针,2、红外热像仪探头,3、微小反应器,4、恒温器,5、进样阀,6、恒流泵,7、数据采集处理模块,8、计算机,9、稳压电源。
具体实施方式
以下结合附图对本发明的技术方案进行详细的描述,以使本领域技术人员能够更加清楚的理解本发明,但并不因此限制本发明的保护范围。
应用本发明技术方案,采用红外热成像测温手段,与探针测量热物理参数相结合,实施了多孔介质堆积固定床微小反应器三维温度分布的测量和重构。微小反应器3为圆柱型,采用同轴套管的形式,套管材料为透红外玻璃管(5~13μm波长范围红外光透过率为0.93),有效长度为100mm,内径分别为10mm和25mm,壁厚为1mm,且内外管间布置间隔5mm环状空气层起热绝缘作用,同时对红外热成像法中辐射的影响可以忽略不计,内外管两端由硅胶垫圈与带有螺纹的SLA堵头密封并固定,堵头内固定有微孔滤布以防止泄露。微小反应器3固定在可调节支架上,周围布置绝热板,减小环境温度辐射对其的影响。***原理如图1所示,在微小反应器3内部堆积填充多孔介质(吸附过氧化氢酶的732型阳离子交换树脂颗粒,0.8~0.9mm)为固定床形式,连接***流路,调节恒温器4控制微小反应器3入口处的温度,进样阀5用于进样,实现微小反应器3内的可控的化学反应,恒流泵6以恒定流量将溶液注入微小反应器3。具体测试步骤如下:首先在微小反应器3中心轴线放置探针1用以在线测量有效导热系数,探针1的加热丝所需恒定电压由稳压电源9提供,在微小反应器3的四周均匀布置三组红外热像仪探头2以全面测量堆积多孔介质的外表面温度;然后利用探针1在线测量微小反应器3内堆积多孔介质的有效导热系数等热物理参数,并得到中心轴线位置的温度,同时调整红外热像仪探头2的角度和焦距,测量微小反应器3外表面的温度分布和变化,由于采用的是透红外玻璃,故红外热像仪探头2能够感测到微小反应器3内堆积多孔介质的外表面温度分布和变化,上述信号经数据采集模块8引入计算机9处理,得到如图2和图3所示微小反应器内堆积多孔介质的典型温度热像图和表面温度分布;最后,基于上述步骤得到的中心温度、有效导热系数等热物理参数和外表面温度分布,依据传热学理论,微小反应器为圆柱型固定床,并假定多孔介质材料均匀且各向同性,流速较小流动为稳态层流,并且为不可压缩流体,流体和多孔介质固相骨架处于局部热平衡状态,得到如下导热微分方程组和单值性条件,数学解析求解后得到微小反应器内堆积多孔介质的三维温度场,如图4所示。
其中:
T为流体的温度;
t为多孔介质(树脂)的温度;
Cg为流体的比热;
Cs为多孔介质(树脂)的比热;
Kc为有效导热系数;
ρg为流体的密度;
ρs为多孔介质(树脂)的比热;
G流体的质量流速;
h为流体与多孔介质(树脂)的换热系数;
ε为多孔介质(树脂)的空隙率;
γ化学反应速率;
ΔH化学反应焓变;
t0初始温度;
ti中心轴线温度;
x为轴向坐标;
r为径向坐标。
τ为时间;
R为微小反应器的半径。
Claims (1)
1.一种多孔介质三维温度分布测量方法,其特征在于,应用于多孔介质堆积固定床微小反应器的三维温度分布测量,包括如下步骤:
1)在圆柱型微小反应器(3)中心轴线放置探针(1)用以在线测量有效导热系数,在所述微小反应器(3)的四周均匀布置三组红外热像仪探头(2)以全面测量外表面温度;
2)利用所述探针(1)在线测量所述微小反应器(3)内多孔介质的有效导热系数,并得到中心轴线位置的温度,同时调整所述红外热像仪探头(2)的角度和焦距,测量所述微小反应器(3)内堆积多孔介质的外表面温度分布;
3)基于上述步骤得到的所述中心轴线的温度、所述有效导热系数和所述外表面温度分布,依据传热学理论,得到导热微分方程组和单值性条件,所述导热微分方程组和单值性条件为:
其中:
T为流体的温度;
t为多孔介质(树脂)的温度;
Cg为流体的比热;
Cs为多孔介质(树脂)的比热;
Kc为有效导热系数;
ρg为流体的密度;
ρs为多孔介质(树脂)的比热;
G流体的质量流速;
h为流体与多孔介质(树脂)的换热系数;
ε为多孔介质(树脂)的空隙率;
γ化学反应速率;
ΔH化学反应焓变;
t0初始温度;
ti中心轴线温度;
x为轴向坐标;
r为径向坐标;
τ为时间;
R为微小反应器的半径;
4)数学解析求解后得到所述微小反应器(3)内堆积多孔介质的三维温度场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711100574.6A CN107843347B (zh) | 2017-11-09 | 2017-11-09 | 一种多孔介质三维温度分布测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711100574.6A CN107843347B (zh) | 2017-11-09 | 2017-11-09 | 一种多孔介质三维温度分布测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107843347A CN107843347A (zh) | 2018-03-27 |
CN107843347B true CN107843347B (zh) | 2019-07-30 |
Family
ID=61681534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711100574.6A Expired - Fee Related CN107843347B (zh) | 2017-11-09 | 2017-11-09 | 一种多孔介质三维温度分布测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107843347B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110987503A (zh) * | 2019-12-19 | 2020-04-10 | 吉林建筑大学 | 一种基于多孔导热材料的热管散热*** |
CN112763014B (zh) * | 2021-02-05 | 2021-09-21 | 江苏华尔威科技集团有限公司 | 基于物联网技术的油井测量*** |
CN113125495B (zh) * | 2021-03-17 | 2022-11-11 | 北京理工大学 | 一种堆积含能材料热导率测量及修正方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710004A (zh) * | 2009-12-11 | 2010-05-19 | 重庆理工大学 | 基于光纤Bragg光栅的生物膜式反应器内温度场测量*** |
CN201732071U (zh) * | 2010-06-10 | 2011-02-02 | 上海理工大学 | 一种测定不良导热材料热导率的装置 |
CN103324781A (zh) * | 2013-05-27 | 2013-09-25 | 华中科技大学 | 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法 |
CN103820631A (zh) * | 2014-02-21 | 2014-05-28 | 中南大学 | 立式淬火炉构件温度场分布检测*** |
CN106508040B (zh) * | 2012-12-31 | 2014-08-20 | 中国人民解放军国防科学技术大学 | 一种多孔热解材料传热过程计算方法 |
CN106248725A (zh) * | 2016-09-16 | 2016-12-21 | 北京工业大学 | 一种多孔介质等效导热系数测量方法 |
CN106320951A (zh) * | 2016-08-15 | 2017-01-11 | 西北工业大学 | 一种内嵌金属泡沫和相变微胶囊材料的强化蓄热百叶窗 |
CN106897537A (zh) * | 2017-03-14 | 2017-06-27 | 清华大学 | 含三维或曲面外形结构体的温度场与热流同时重构方法 |
CN106913316A (zh) * | 2015-12-28 | 2017-07-04 | 精工爱普生株式会社 | 内部温度测定装置、手腕安装式装置及内部温度测定方法 |
CN106960089A (zh) * | 2017-03-14 | 2017-07-18 | 清华大学 | 含内部复杂边界结构体的温度场和热流同时重构方法 |
-
2017
- 2017-11-09 CN CN201711100574.6A patent/CN107843347B/zh not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710004A (zh) * | 2009-12-11 | 2010-05-19 | 重庆理工大学 | 基于光纤Bragg光栅的生物膜式反应器内温度场测量*** |
CN201732071U (zh) * | 2010-06-10 | 2011-02-02 | 上海理工大学 | 一种测定不良导热材料热导率的装置 |
CN106508040B (zh) * | 2012-12-31 | 2014-08-20 | 中国人民解放军国防科学技术大学 | 一种多孔热解材料传热过程计算方法 |
CN103324781A (zh) * | 2013-05-27 | 2013-09-25 | 华中科技大学 | 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法 |
CN103820631A (zh) * | 2014-02-21 | 2014-05-28 | 中南大学 | 立式淬火炉构件温度场分布检测*** |
CN106913316A (zh) * | 2015-12-28 | 2017-07-04 | 精工爱普生株式会社 | 内部温度测定装置、手腕安装式装置及内部温度测定方法 |
CN106320951A (zh) * | 2016-08-15 | 2017-01-11 | 西北工业大学 | 一种内嵌金属泡沫和相变微胶囊材料的强化蓄热百叶窗 |
CN106248725A (zh) * | 2016-09-16 | 2016-12-21 | 北京工业大学 | 一种多孔介质等效导热系数测量方法 |
CN106897537A (zh) * | 2017-03-14 | 2017-06-27 | 清华大学 | 含三维或曲面外形结构体的温度场与热流同时重构方法 |
CN106960089A (zh) * | 2017-03-14 | 2017-07-18 | 清华大学 | 含内部复杂边界结构体的温度场和热流同时重构方法 |
Non-Patent Citations (1)
Title |
---|
探针法测量微细颗粒固定床有效热导率;张心怡 等;《化工进展》;20170605;第36卷(第6期);第2010-2016页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107843347A (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107843347B (zh) | 一种多孔介质三维温度分布测量方法 | |
CN104569045B (zh) | 圆柱套筒壁间结合面接触热阻测试方法及装置 | |
CN101113963A (zh) | 一种测量液体导热系数的方法及其装置 | |
CN102539019A (zh) | 空间真空环境下的温度测量与校准平台 | |
CN103728337B (zh) | 测量物体内部热流密度的热流密度探针以及测量方法 | |
Estrada et al. | Heat transfer analysis in a calorimeter for concentrated solar radiation measurements | |
CN115452180B (zh) | 一种高焓气流恢复温度测量方法及测量装置 | |
CN106153672A (zh) | 基于一维导热原理的多孔粉末材料热导率测量装置及方法 | |
CN103983365B (zh) | 多测头瞬态辐射热流计及热辐射热流密度的测定方法 | |
Chouikh et al. | Experimental study of the natural convection flow around an array of heated horizontal cylinders | |
CN102778475B (zh) | 一种上下恒温参数辨识法测固-固接触热阻 | |
CN206756728U (zh) | 一种良导体导热系数实验仪 | |
Reynard-Carette et al. | Review of nuclear heating measurement by calorimetry in France and USA | |
CN108445042A (zh) | 一种测量建筑外表面对流换热系数的方法 | |
CN207689406U (zh) | 固体材料热扩散系数随温度变化的测量装置 | |
Jaremkiewicz et al. | Measuring transient temperature of the medium in power engineering machines and installations | |
Ritterath et al. | Thermo-resistive mesh sensors (TMS) for temperature field measurements | |
CN107064548A (zh) | 一种传感器装置及测量方法 | |
Sima et al. | Accurate numerical model for characteristic temperature acquisition of miniature fixed-point blackbodies | |
Hohmann et al. | Calibration of heat flux sensors with small heat fluxes | |
CN108152326A (zh) | 固体材料热扩散系数随温度变化的测量装置及方法 | |
CN202947983U (zh) | 表面传热系数的测量装置 | |
CN102062651B (zh) | 热量表温度动态检定方法 | |
CN110376244A (zh) | 一种导热系数测量装置 | |
Dietrich et al. | Thermal conductivity of ceramic sponges at temperatures up to 1000 C |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190730 Termination date: 20201109 |
|
CF01 | Termination of patent right due to non-payment of annual fee |