一種金屬表面發射率測量的新方法

張航 李圣慧 嚴小亮 盧欣欣

摘 要：【目的】針對測量管殼式熱交換器的管束表面發射率精度低、成本高的問題，提出了一種廉價且精準的表面發射率測量的新方法。【方法】該方法以馬弗爐內部空間為獨立的封閉系統，根據能量守恒定律推導出系統內被測對象的表面發射率的函數關系式。對溫度采樣的過程中的低頻噪聲污染，采用兩種平滑算法進行降噪。【結果】溫度設置為200 ～ 500 ℃，對三種不同材質、不同表面處理方式的鋼材進行了溫度測量和發射率計算與誤差分析。對比結果后發現：隨溫度和尺寸的升高，奧氏體不銹鋼發射率逐漸減小，而碳素鋼發射率逐漸增大。【結論】新方法可實現表面發射率的在線測量，且誤差低。通過溫度數據的降噪優化算法可以提高精度。奧氏體不銹鋼的發射率與樣品的尺寸有關，尺寸越大，發射率越大。

關鍵詞：管殼式換熱器；鋼材料；表面發射率；輻射傳熱

中圖分類號：TG806? ? ?文獻標志碼：A? ? ?文章編號：1003-5168（2024）07-0081-07

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.07.016

A New Method for Metal Surface Emissivity Measurement

ZHANG Hang1，2 LI Shenghui1 YAN Xiaoliang1 LU Xinxin1

（1.School of Aviation and Transportation， Jiangsu College of Engineering and Technology， Nantong 226006， China; 2.Liupanshui Shuangyuan Aluminum Industry Limited， Liupanshui 553022， China）

Abstract： [Purposes] Aiming at the problem of low precision and high cost of measuring the surface emissivity of tube bundles of shell-and-tube heat exchangers， this paper proposes a new method of cheap and accurate surface emissivity. [Methods] The method takes the internal space of the muffle furnace as an independent closed system， and derives the functional equation of the surface emissivity of the measured object in the system according to the law of energy conservation. Two smoothing algorithms are used for noise reduction of low-frequency noise pollution in the process of temperature sampling. [Findings] Temperature measurement and emissivity calculation and error analysis were carried out for three different materials and surface treatments of steel in the temperature range of 200 ～ 500 ℃. Comparison of the results shows that the emissivity of austenitic stainless steel gradually decreases with increasing temperature and size， while the emissivity of carbon steel gradually increases. [Conclusions] The new method can realize online measurement of surface emissivity with low error. The accuracy can be improved by the noise reduction optimization algorithm of temperature data. The emissivity of austenitic stainless steel is related to the size of the sample， and the larger the size， the larger the emissivity.

Keywords： shell-and-tube heat exchanger；steel material； surface emissivity； radiative heat transfer

0 引言

管殼式熱交換器是煙氣余熱回收利用的主要技術和設備之一，其換熱效率直接取決于管束材質的傳熱性能，即鋼材的升溫速率［1］。熱輻射是其升溫的主要方式之一，因此，鋼的輻射性能直接關系到煙氣熱交換器的工作效率［2］。

發射率是描述物體熱輻射特性的重要參數，發射率越大，管束表面的熱輻射越強烈，熱量傳遞的速度也越快［3-4］。選用具有高發射率的材料可以增強熱交換器的熱輻射能力，提高設備的能源利用效率。發射率的大小與材料組分、溫度、表面粗糙度和氧化程度等因素密切相關，實際情況非常復雜［5-8］。目前，獨立黑體法是發射率測量技術中應用最廣泛的方法之一，但它仍不能滿足工業上的精度要求［9-12］。例如，Cai等［13］基于能量比較法搭建了一套發射率測量系統，其測量不確定度為4.2%。Adibekyan等［14］提出用液氮冷卻關鍵光路組件以減少背景輻射，從而提高發射率精確度的方法。這種方法對于碳化硅的方向光譜發射率的測量標準不確定度可以達到≤0.01的水平。但該方法成本昂貴，需要光譜儀才能進行。不僅如此，上述研究及其他研究者［15-18］采用的是基于能量比較法的間歇式測量，測量裝置技術參數見表1。操作人員在高溫下處理材料不方便，也不安全，且這類方法對波長范圍和溫度范圍的應用非常有限。

以上方法在過去粗放式的生產中可以解決一定的問題。隨著節能環保要求的不斷提高，急需開發出一種更精確的發射率分析方法，以便充分利用鋼材料在室溫至工作溫度范圍內的輻射特性。本研究開發一種簡單廉價的方法來連續測定鋼材料的發射率，并找到它們與溫度的函數關系。該方法將室溫中樣品迅速放入高溫的爐中，記錄樣品的升溫曲線，并使用推導出的材料發射率公式計算得出發射率。

1 物理與數學模型

將室溫下的球形金屬物體，在初始時刻置入一個加熱的球形空腔中，試驗原理示意如圖1所示。金屬球的初始溫度設為T1，空腔初始溫度設為T2。T1和T2都是時間t的函數。

1.1 熱輻射

現對兩個同心球面之間的輻射傳熱進行計算。內球面的平均溫度設為T1，外球面的平均溫度設為T2。內球面的表面積和發射率分別設為A1和ε1，外球面的表面積和發射率分別設為A2和ε2。假設球面是漫反射灰體。使用凈輻射方法求解輻射傳熱量［19］，內球面的凈輻射熱流見式（1）。

[Qrad=AσT42-T411ε2T2+A1A21ε2T2-1] （1）

式中：σ為Stefan-Boltzmann常數，5.67×10-8 w/（m2K4）；Qrad正值表示內球面是接收熱量（T2>T1），而Qrad負值表示內球面是輸出熱量（T2

[Qrad=ε1T1AσT42-T41] （2）

即凈輻射熱流與外球面的發射率無關。

1.2 熱傳導

現計算兩個同心球面之間的熱傳導。球面之間的流體熱導率為λA。內球面和外球面的直徑分別為D1和D2。內球面的熱傳導表達式見式（3）。

[Qcond=2πλAD1T2-T11D1-1D2=2πλAD1T2-T11-D1D2] （3）

當外球面的直徑遠大于內球面時，式（3）可以簡化為式（4）。

[Qcond=2πλAD1T2-T1] （4）

1.3 熱對流

假設該球體放置在大空間的流體中。在遠離物體的地方，假設流體靜止，且溫度為T2。流體和物體之間的對流傳熱表達式見式（5）。

[Qconv=hconvA1T2-T1] （5）

式中：hconv為流體和球體之間的對流傳熱系數。當球體吸熱時，符號為正。流體與流體中球體之間的自然對流傳熱方程給出對流傳熱系數的求解方法［20］，見式（6）。

[Nu=1.414+0.387Ra1/61+0.492Pr9/168/272] （6）

式中：Nu為Nusselt數；Ra為Rayleigh數，Pr為Prandtl數。它們的定義見式（7）至式（10）。

[Nu=hConvLλA] （7）

式中：L為球體的特征尺寸，即球體的直徑。λA為流體的熱導率。

[Ra=Gr?Pr] （8）

式中：Gr為Grashof數。如果流體可以被假定為理想氣體，那么Gr的定義見式（9）。

[Gr=gL3v2AT2-T1T2] （9）

式中：g為萬有引力常數（9.81 m/s2），νA為流體的運動黏度。Pr的定義見式（10）。

[Pr=μACpAλA] （10）

式中：μA為動力黏度；CpA為等壓比熱容。將式（7）至式（10）代入式（6）即可求解對流傳熱系數hconv。流體的所有熱學性質都是在流體和球體之間的平均溫度下取得的。下面比較對流傳熱和傳導傳熱的相對大小。將式（7）代入式（5），可得到式（11）.

[QConv=hConvA1T2-T1=NuπλAD1T2-T1]

（11）

熱對流與熱傳導的比值見式（12）。

[QConvQcond=Nu2] （12）

當T1趨近于T2，溫差引起的流動降低，QCond接近QConv。當Gr趨近于0時，Nu趨近于2，則 QCond ≈ QConv。當ΔT=T1-T2變大，Gr和Nu值增大，這反映了對流換熱增強。

1.4 能量守恒方程

在任一時刻，球體都處于能量平衡狀態，能量平衡方程見式（13）。

[m1cp1dT1dt=iQi] （13）

式中：m1為物體的質量，cp1為物質的比熱，t為時間，Qi為物體與環境之間的換熱量。Qi包括輻射熱 Qrad和對流熱Qconv。因此，式（13）可以變換為式（14）。

[m1cp1dT1dt=QRad+QConv=ε1T1A1σT42-T41+hconvA1T2-T1] （14）

由此可得發射率ε1（T1）的表達式見式（15）。

[ε1T1=m1cp1dT1dt-hconvA1T2-T1A1σT42-T41] （15）

式中：m1，σ，A1和cp1已知。在試驗過程中記錄下T1（t）和T2（t）后，hconv通過式（6）計算得到，就可以確定材料的發射率ε1（T1）。

2 試驗過程設計

2.1 試驗方案

馬弗爐的爐膛尺寸：長度400 mm，寬度200 mm，高度180 mm，表面積為0.376 m2。馬弗爐采用雙層陶瓷，四面加熱，配備數字式PID控制器，可調節升溫曲線，溫度上限為1 000 ℃，如圖2所示。

2.2 樣品處理

試驗樣品將使用三種不同的材料，分別為ST-1，ST-2和CT-3。測試材料ST-1是從厚度為5.0 mm的奧氏體不銹鋼上切割下來的。測試材料ST-2是從厚度為2.0 mm的薄板上切下的冷軋不銹鋼。測試材料CT-3是從厚度為2.0 mm的薄板上切下的冷軋碳鋼。表面處理工藝見表2。

樣品制備：在材料厚度方向上切一個3 mm深的矩形槽，用來放置K型熱電偶，如圖3所示。

采用K型熱電偶測量爐內空氣和樣品的溫度，用于測量爐溫的三個熱電偶被固定在鋼絲網。

2.3 數據處理

在爐溫到達設定溫度15 min后打開爐門，將樣品及其支撐架放入，關閉爐門。當樣品的溫度開始接近爐溫時，再次打開爐門，將樣品和支撐架從爐中取出。最后，終止記錄數據。當爐門打開時存在低頻噪音的污染，因此采用了兩種不同的方法（常規方法和優化方法）進行額外的平滑處理。

常規方法：①每個樣品稱重并記錄。②游標卡尺測量樣品尺寸，計算表面積。③假設測試材料ST-1和ST-2的比熱是溫度的線性函數［21］。而CT-3的比熱是溫度的三次多項式擬合［22］。④樣品溫度由貼在樣品上的熱電偶測得。⑤爐溫是位于爐膛三個不同高度的熱電偶測得的溫度平均值。⑥樣品溫度dT1/dt的時間導數按以下方式計算：首先，每個時間點及其前后各四個單位時間的平均溫度作為當前時間的樣品溫度；其次，在ti時刻的樣品溫度導數由公式dT1/dt=（T1，i+1-T1，i-1）/（2Δt）算出。⑦利用式（2）、式（5）和式（6）至式（10）計算出樣品與爐膛內空氣之間的對流傳熱系數。爐內空氣的熱力學參數是在爐膛和樣品的平均溫度下取得［23］。

優化方法：步驟①～④、⑦與上述過程一致，主要區別在于步驟⑤、⑥。樣品溫度是樣品表面溫度與爐溫的平均值。爐溫的采樣點是三個不同高度，三個不同深度的空間點。樣品溫度dT1/dt的時間導數按以下方式計算：首先，每個時間點的溫度是當前溫度與前后各4個點的溫度變化率的平均值的乘積，再加上當前溫度作為當前時間的樣品溫度。然后，在ti時刻的樣品溫度導數由公式dT1/dt=（T1，i+1-T1，i-1）/（2Δt）算出。

ST-1的升溫曲線及發射率如圖4所示。常規方法獲得的溫度，如圖4（a）所示。在爐門開啟時，樣品溫度受熱電偶采樣方式的影響出現劇烈下降，這顯然與正常現象相違背。在爐門關閉后，樣品溫度測量受熱電偶性能的影響出現輕微浮動。使用常規算法對噪聲污染的數據進行分析得到發射率的溫度在達到400 ℃后，出現振蕩，說明常規方法對噪聲的處理不能滿足要求。

從優化方法的處理結果可以看出，樣品進出爐時的溫度變化過程更加緩慢。優化算法的結果比常規算法得到的結果更平滑，如圖4（b）所示。

2.4 誤差分析

誤差分析將計算發射率值的不確定度u（ε），作為式（15）右側測量值不確定性的函數。測量不確定度是一個參數，它表征了測量值的合理范圍。如果樣品溫度值的不確定度小于4 ℃，發射率的不確定度小于±2.0%被認為是合理的［24］。對于如下關系：

[y=x1+x2] （16）

與y相關的相對不確定度見式（17）。

[uy=ux12+ux22] （17）

對于如下關系：

[y=xm1xn2] （18）

與y相關的相對不確定度見式（19）。

[uyy=m2u（x1）x12+n2u（x2）x22] （19）

將式（17）和式（19）代入式（15），可得到發射率值的不確定度。輸入量的不確定性估計值見表3。不確定度的主要來源是爐溫的測量、樣品溫度的時間導數的估算和對流傳熱系數的估算。

在被測溫度范圍內，樣品溫度升高，不確定度開始增加，在480 ℃時達到3.0%，如圖5所示。原因是隨著樣品溫度開始接近爐溫，樣品溫度的時間導數變得很小，因此樣品溫度時間導數的相對不確定度成為主導因素。

3 結果與分析

3.1 爐溫的影響

在試驗中，爐溫在200～500 ℃之間增加時，樣品ST-1和ST-2發射率是逐漸降低的，而樣品CT-3的卻逐漸升高。熱輻射換熱計算的方程中包含溫度的四次方，這意味著在溫度的微小變化對熱輻射換熱的影響巨大，如圖6所示。其次，優化算法對于樣品ST-1和ST-2在爐溫變化時的發射率仍然波動，但相比常規算法有明顯的改進。但是對于樣品CT-3，兩種算法的發射率變化的差異并不明顯。

3.2 樣品大小的影響

樣品ST-1的長度分為25 mm和50 mm時的發射率如圖7所示。長度上的差異會導致對流傳熱系數發生相當大的變化。對于ST-1，長度為50 mm和25 mm時，優化算法獲得的發射率ε的范圍分別為0.265和0.29。

3.3 不同樣品材料的對比

不同測試材料的發射率如圖8所示。由圖8可知，樣品ST-1和ST-2的發射率是溫度的弱遞減函數。另一組材料CT-3在低溫下發射率相當低，并且與溫度無關。當溫度超過350 ℃時，CT-3的發射率隨之增加，直到溫度超過500 ℃時，發射率與溫度的遞增關系消失。

圖8（a）的結果表明：優化算法得到ST-1的發射率ε≈0.305，常規算法得到εn≈0.22～0.28，常規算法獲得的結果不確定度約為12.5%～25.0%，而優化算法的不確定度約為1.5%～2.5%。圖8（b）表示的是樣品ST-2的發射率。結果表明，優化算法得到ε≈0.255，常規算法得到εn≈0.20～0.23。同樣地，優化算法的不確定度約為4.5%～7.2%。圖8（c）表示的是CT-3樣品的發射率。結果表明，發射率與溫度范圍有明顯關聯。在溫度高于330 ℃時，發射率開始增加。這可能是碳素鋼的表面處理工藝導致樣品發射率在該階段溫度的異常表現。

4 結論

輻射是管殼式換熱器的一種主要換熱形式，發射率是輻射換熱的一個重要參數，因此發射率的測量對于提高管殼式換熱的工作效率具有重要意義。本研究提出了一種測量與計算鋼材料發射率的新方法，并通過試驗得出如下結論。

①對于奧氏體不銹鋼，測定的溫度范圍從200 ℃至500 ℃，發射率逐漸降低，此外，樣品的長度和厚度尺寸發生變化時，樣品的發射率也會有所變化。長度越長，發射率越大；厚度越厚，發射率越大。通過優化算法計算的發射率不確定度約為4.5%～7.2%，相同情況下常規算法的發射率不確定度約為12.5%～25.0%。優化算法的發射率不確定度顯著低于常規算法。

②對于碳素鋼，當溫度區間為200 ℃至500 ℃時，通過優化算法和常規算法計算得到的發射率不確定度接近。發射率在從250 ℃至450 ℃區間，隨著溫度的升高而增加，在其他溫度下則非常小。

③通過優化方法計算得到的發射率普遍高于常規方法。

參考文獻：

［1］邱爽.管殼式換熱器失效原因及其預防措施［J］.化學工程與裝備，2022（4）：190-191，212.

［2］張禹，褚浩然，徐衛，等.煙氣水冷換熱器考慮熱輻射影響時的數值模擬［J］.化學工程，2020，48（3）：42-47，78.

［3］董珂馨，李寶寬，徐健祥，等.雙排料步進梁式加熱爐壁面發射率對鋼坯溫度均勻性影響［J］.中國冶金，2023，33（9）：118-127.

［4］莊亞平，董海濤，潘姚凡.鹽霧對鋼制安全殼涂層傳熱性能影響研究［J］.核科學與工程，2022，42（1）：185-191.

［5］馮馳，劉思源.渦輪葉片發射率測量［J］.哈爾濱商業大學學報（自然科學版），2020，36（5）：576-581.

［6］熊磊，曾鳴，周波，等.涂層半球發射率影響因素的研究［J］.節能技術，2021，39（5）：426-431.

［7］王鎮，周海洋，羅蛟.反射隔熱涂層半球發射率的影響因素［J］.上海涂料，2021，59（2）：14-16.

［8］張嵐，蔡靜，路林鋒.影響輻射測溫用黑體輻射源應用的因素分析［J］.計測技術，2020，40（4）：1-4.

［9］林芳，胡進坤，徐海濤，等.不同表面粗糙度金屬的輻射發射率實驗設計［J］.實驗技術與管理，2018，35（12）：43-47，61.

［10］龍超，陳軍燕，楊雨川，等.溫度變化下研磨金屬表面反射率和發射率的測量［J］.發光學報，2016，37（12）：1566-1570.

［11］魏薇，趙卉.金屬材料表面發射率測量方法研究［J］.工業計量，2015，25（6）：55-56，59.

［12］呂方明，王坤，黃樹紅，等.超臨界水蒸氣氧化的12Cr鋼發射率試驗研究［J］.機械工程學報，2012，48（18）：124-129.

［13］CAI J， YANG Y， LIAO L， et al. Material Spectral Emissivity Measurement Based on Two Reference Blackbodies［J］.InternationalJournalofThermophysics，2015，36：3288-3296．

［14］ADIBEKYAN A， MONTE C， KEHRT M， et al. Emissivity Measurement Under Vacuum from 4μm to 100μm? and from ?40℃ to 450℃ at PTB［J］.InternationalJournalofThemophysics， 2015， 36：283-289.

［15］HONNEROV? P， MARTAN J， KU?ERA M， et al. Newexperimentaldeviceforhigh-temperaturenormalspectral emissivity measurements of coatings［J］. Measurement ScienceandTechnology， 2014， 25（9）：30-37.

［16］WEN C D. Study of Steel Emissivity Characteristics and Application of Multispectral Radiation Thermometry （MRT）［J］. Journal of Material Engineering and Perform， 2011， 20： 289-297.

［17］H?SER D， WALLIMANN R， VONROHR P R. Uncertaintyanalysis for emissivity measurement at elevated temperatureswithan infraredcamera［J］. International Journal ofThermophysics，2016，37：14.

［18］WOODS S， JUNG T， LY G， KAPLAN S， et al. Broadband emissivity calibration of highly reflective samples at cryogenic temperatures［J］. Metrologia，2012， 49： 737-744.

［19］王怡，王宇琨，孟曉靜.高溫熱源工業建筑室內輻射熱在各壁面分布研究［J］.西安建筑科技大學學報（自然科學版），2020，52（2）：302-308.

［20］PATEL P A， BHATNAGAR D， KUMAR S R， et al. Numerical study on turbulent natural convection and radiation heat transfer of nanofluids in a differentially heated square enclosure［J］. Journal of Thermal Analysis and Calorimetry，2020，147（1）： 675-684.

［21］陳鵬，王如梅，王征，等.金屬比熱容測量和溫度傳感器測試熱學實驗的對比教學探討［J］.大學物理實驗，2021，34（1）：27-30，59.

［22］JIAO Y， LIU C F， CUI X P， et al. A new approach for measurement of the low-temperature specific heat capacity［J］. Measurement，2022，203（15）： 111892.

［23］DEHAGHANI A H S， SOLEIMANI R， MOHAMMADI A H. Screening of Important Parameters in Optimal Design of Compressed Air Energy Storage System Using an Ensemble Learning Method［J］.Journal of Energy Storage，2022，48：104023.

［24］袁良，袁林光，董再天，等.高溫狀態下的材料法向光譜發射率測量［J］.應用光學，2023，44（3）：580-585.

收稿日期：2023-11-21

基金項目：江蘇省普通高校自然科學研究計劃重大項目（23KJA460004）；江蘇省普通高校自然科學研究計劃面上項目（22KJD480001）；江蘇省高職院校青年教師企業實踐培訓項目（2023QYSJ021）；南通市基礎科學研究計劃項目青年科技人才創新專項（JC12022019）；江蘇工程職業技術學院自然科學研究基金項目（GYKY/2022/6、GYKY/2022/3）

作者簡介：張航（1987—），男，博士，講師，研究方向：新能源裝備與控制工程。

通信作者：李圣慧（1989—），女，講師，研究方向：新能源政策分析。