楊慶濤，白菡塵，劉濟(jì)春，王 輝

(1.中國空氣動力研究與發(fā)展中心 超高速所，四川 綿陽621000;2.中國空氣動力研究與發(fā)展中心 超高速所 高超聲速沖壓發(fā)動機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 綿陽621000)

0 引 言

薄壁量熱計(jì)是氣動熱與熱防護(hù)試驗(yàn)中常用的一種瞬態(tài)熱流傳感器，既可用于脈沖式試驗(yàn)，也可用于防熱材料與熱結(jié)構(gòu)試驗(yàn)設(shè)備［1］，其測熱原理是基于后壁絕熱的假定，這在實(shí)際測量中是很難實(shí)現(xiàn)的。后壁與相鄰環(huán)境之間總會存在熱量交換，從而影響薄壁量熱計(jì)的溫度響應(yīng)，造成測量誤差。

熱損失對測量影響的研究主要針對不同邊界條件分別進(jìn)行分析。Mahta R C 等人［2］利用有限元方法求解了一邊為氣動加熱條件，另一邊為輻射邊界條件的二維導(dǎo)熱方程，并給出了考慮后表面輻射影響的熱流換算關(guān)系式。Kidd C T［3］對熱電偶線熱傳導(dǎo)造成的熱損失進(jìn)行了分析，表明薄壁向熱電偶線的導(dǎo)熱會使傳感器響應(yīng)時(shí)間增大。Keltner N R 等人［4］給出了不同薄壁和熱電偶材料組合下，熱電偶導(dǎo)熱造成的溫升率誤差計(jì)算公式。Lewis D R 等人［5］數(shù)值分析了不同直徑的E 型熱電偶導(dǎo)熱對不銹鋼薄壁量熱計(jì)測量的影響，指出熱電偶線的存在起到熱沉的作用，導(dǎo)致測得熱流低于實(shí)際施加熱流，增加量熱計(jì)的響應(yīng)時(shí)間，導(dǎo)熱性更好的熱電偶線誤差會更大。Schultz D L 和Jones T V［6］給出了在薄壁向相鄰物質(zhì)導(dǎo)熱的情況下，薄壁前后表面溫度響應(yīng)的解析表達(dá)式，以及基體導(dǎo)熱熱流的估算表達(dá)式，但文中的表面溫度響應(yīng)表達(dá)式不正確。

本文利用理論分析方法，對比后壁有導(dǎo)熱損失與無導(dǎo)熱損失兩種情況下薄壁的溫度響應(yīng)，研究考慮后壁導(dǎo)熱損失時(shí)的熱流修正方法，并對修正效果進(jìn)行評價(jià)分析。

1 帶基底材料的薄壁量熱計(jì)溫度響應(yīng)

典型的薄壁量熱計(jì)結(jié)構(gòu)如圖1 所示。傳感器后壁采用空氣隔熱時(shí)，如果可以忽略對流換熱的影響，可將空氣作為導(dǎo)熱系數(shù)非常低的隔熱基底材料［7］。如果傳感器所處熱環(huán)境氣體壓力較大，為避免薄壁發(fā)生變形，常采用導(dǎo)熱系數(shù)低的固體基底材料，如水泥、陶瓷、樹脂等。在一定時(shí)間內(nèi)可以將基底材料近似為半無限大體，用圖2 所示的物理模型來表示，其中，薄壁材料厚度為l，密度、比熱和導(dǎo)熱系數(shù)分別記為ρ1，c1，k1;基底材料相應(yīng)的熱物性參數(shù)記為ρ2，c2，

圖1 薄壁量熱計(jì)測量熱流原理示意圖Fig 1 Diagram of heat flux measurement principle of thin skin calorimeter

圖2 兩層介質(zhì)傳熱模型Fig 2 Heat transfer model for two-layer medium

對于一維熱傳導(dǎo)，兩層介質(zhì)導(dǎo)熱微分方程

初始條件為

邊界條件為

式中 T 為材料溫升，α=k/(ρc)為材料的熱擴(kuò)散率，qw為輸入熱流。

利用拉普拉斯變換［6，8］，方程組(1)可以求解。前后表面的溫度變化率利用傅里葉數(shù)Fo 表示:

前表面

后表面

其中，F(xiàn)o=α1t/l2為以薄膜材料的熱物性和厚度表示的傅里葉數(shù)，熱物性比表征薄壁和基底材料熱物性的相對差別。

將式(2)和式(3)對時(shí)間t 進(jìn)行積分，可以得到存在基底材料時(shí)薄壁前后表面溫度隨時(shí)間的變化?？梢娪绊懕”诒砻鏈囟鹊囊蛩赜?熱流qw、薄壁集總熱參數(shù)、基底材料的、薄壁厚度l 和無量綱時(shí)間Fo。

常見薄壁量熱計(jì)的熱流數(shù)據(jù)處理假定除感應(yīng)面以外的其他面與外界絕熱，氣流傳遞給量熱計(jì)的熱量僅造成量熱計(jì)材料的內(nèi)能增加。忽略基底材料導(dǎo)熱的影響，式(2)、式(3)可以簡化為

其中，Tf為前表面溫度;Tb為后表面溫度。非穩(wěn)態(tài)導(dǎo)熱進(jìn)入正規(guī)狀況階段后，材料內(nèi)部的溫度分布只取決于邊界條件和材料物性，不再受初始溫度分布影響，材料內(nèi)部任意位置的溫度隨時(shí)間的變化處處相等［9］。此時(shí)式(4)和式(5)所得前后表面變化率相等，qw可表示為

式(6)是工程上最常用的薄壁熱流測量公式［10］，也是式(2)和式(3)在理想條件下的一個(gè)特例。

利用軸對稱有限元分析模擬了帶基底材料的薄壁量熱計(jì)的溫度響應(yīng)，所選用薄壁材料為不銹鋼，入射熱流為0.5 MW/m2，厚度l=1 mm。數(shù)值計(jì)算結(jié)果與上述公式的計(jì)算結(jié)果符合良好。

2 后壁導(dǎo)熱損失的影響與修正方法

2.1 后壁損失的影響估算

利用薄壁量熱計(jì)測量熱流，通常利用式(6)處理溫度數(shù)據(jù)，在相同熱流輸入條件下，后壁損失造成的測量誤差可以用后壁絕熱條件下的平均溫度變化率和真實(shí)情況下的溫度變化率測量值來表示

其中，(dTm/dt)|no－sub為后壁絕熱條件下薄壁平均溫度變化率，當(dāng)qw為階躍熱流時(shí)為定值;(dT/dt)|sub為存在導(dǎo)熱損失時(shí)的溫度變化率，可取前、后表面溫度或平均溫度。

進(jìn)行熱流測量時(shí)，輸入熱流qw為待測量，(dTm/dt)|no－sub亦為未知數(shù)，因此，式(7)并不能直接用于誤差估算。薄壁厚度比較小，如果假定后壁溫度與薄壁平均溫度Tm相等，則E 等于后壁導(dǎo)熱熱流qsub與入射熱流qw的比值［6］

可見E 取決于傅利葉數(shù)Fo 和集總熱參數(shù)比a。

實(shí)際測量中一般利用熱電偶測量后壁溫度Tb，用于計(jì)算qw。利用式(3)求得dTb/dt，令dT/dt|sub=dTb/dt，代入式(7)可求出E 值隨Fo 的變化。圖3 為薄壁材料為不銹鋼、l=1 mm 時(shí)，E 隨Fo 的變化?？梢奱 值對測量誤差有明顯影響，E 隨著a 的增大而增大。a=0.1 時(shí)，F(xiàn)o=0.64 時(shí)E 最小，約為10.1%;a=0.25 時(shí)，F(xiàn)o=0.57 時(shí)E 最小，約為21.5%。達(dá)到最小值后，測量誤差隨著Fo 增大而增大。圖4 為測量誤差最小值及對應(yīng)的傅利葉數(shù)隨a 的變化。

圖3 測量誤差E 隨a 值的變化Fig 3 Measurement error E value change with a value

圖4 最小測量誤差E 及對應(yīng)傅里葉數(shù)隨a 值的變化Fig 4 The minimum measurement error E and Fo change with a value

2.2 后壁損失的修正

在實(shí)際的測量過程中，薄壁量熱計(jì)一般測量的是金屬薄壁后壁的溫度。近似認(rèn)為后壁溫度變化率等于薄壁平均溫度變化率，根據(jù)式(8)，可得到修正的薄壁溫度變化率

式中 (dT/dt)meas為實(shí)際測得的薄壁后壁溫度變化率;(dT/dt)|corrected為修正后的溫度變化率。

利用式(9)得到的(dT/dt)|corrected代入式(6)，可得到考慮基底材料導(dǎo)熱的修正熱流。

對于薄壁材料為不銹鋼，l=1 mm，qw=5×105W/m3，a=0.1和a=0.3 兩種條件下，利用式(3)計(jì)算出dTb/dt后，利用式(9)計(jì)算出修正后的溫度變化率，代入式(6)計(jì)算出修正的熱流值。修正前后的溫度變化率和歸一化熱流值(計(jì)算熱流與入射熱流之比:q/qw)見圖5?？梢娎檬?9)修正出的熱流值大大改善了熱流測量準(zhǔn)確度。由于向基底傳熱造成的測量誤差，修正后a=0.1 時(shí)小于1%(修正前＞10%);a=0.3 時(shí)小于5%(修正前＞25%)。因?yàn)楹雎粤撕蟊跍囟扰c薄壁平均溫度的差異，所以，修正的熱流不能完全消除基底導(dǎo)熱測量誤差。

圖5 修正前后的溫度變化率和歸一化熱流Fig 5 Temperature changing rate and normalized heat flux before and after correction

3 試驗(yàn)與結(jié)果分析

在中國空氣動力研究與發(fā)展中心的600 mm 脈沖燃燒風(fēng)洞上，利用薄壁量熱計(jì)和同軸熱電偶測量了軸對稱發(fā)動機(jī)模型(圖6)壁面上的熱流。試驗(yàn)中空氣與氫氣在燃燒器中燃燒升溫，通過噴管膨脹加速后進(jìn)入燃燒室。在不往燃燒室噴注燃料時(shí)，燃燒室壁面壓力保持基本穩(wěn)定，在來流參數(shù)穩(wěn)定的情況下，其壁面熱流接近階躍輸入熱流。

試驗(yàn)中所用的薄壁量熱計(jì)材料為純銅，厚度為0.5 mm，基底隔熱材料為耐溫粘接劑。試驗(yàn)中相同軸向位置上的同軸熱電偶和薄壁量熱計(jì)測得的溫升曲線見圖7。

圖6 軸對稱發(fā)動機(jī)模型Fig 6 Axisymmetric engine model

圖7 溫升曲線Fig 7 Temperature rising curve

相同軸向位置上薄壁量熱計(jì)所得熱流與同軸熱電偶所得熱流對比見圖8，為便于分析對比，圖中是利用快速傅里葉變換(FFT)對曲線進(jìn)行平滑處理后的結(jié)果。其中同軸熱電偶利用Cook－Felderman 公式［10］計(jì)算熱流，薄壁量熱計(jì)利用前述方法分別得到修正前后熱流。試驗(yàn)中同軸熱電偶的響應(yīng)比薄壁量熱計(jì)要快;在線性響應(yīng)段，同軸熱電偶測得平均熱流為1.51 MW/m2，薄壁量熱計(jì)測得平均熱流修正前為1.34 MW/m2，修正后為1.56 MW/m2。

圖8 熱流測量結(jié)果Fig 8 Heat flux measurement results

4 結(jié) 論

1)如果忽略向基底材料的導(dǎo)熱，實(shí)際的薄壁量熱計(jì)測量的熱流會低于真實(shí)入射熱流。

2)后壁導(dǎo)熱造成的誤差取決于傅里葉Fo 和薄壁與基底材料的熱物性比

3)一定時(shí)間以后，導(dǎo)熱相對熱損失達(dá)到最小值，之后隨著Fo 增大而增大。導(dǎo)熱相對損失隨a 的增大而增大。

4)利用修正公式對測量溫度數(shù)據(jù)進(jìn)行處理，可以減小薄壁量熱計(jì)的測量誤差。

［1］ ASTM standard E459－97.Standard test method for measuring heat transfer rate using a thin－skin calorimeter［S］.

［2］ Mehta R C，Jayachandran T，Sastri V M K，et al.Finite element analysis of conductive and radiative heating of a thin skin calorimeter［J］.W?rme－und Stofffübertragung，1988(22):227－230.

［3］ Kidd C T.Thin－skin technique heat－transfer measurement errors due to heat conduction into thermocouple wires［J］.ISA Transactions，1985，24(2):97－99.

［4］ Keltner N R，Beck J V.Surface temperature measurement errors［J］.Journal of Heat Transfer，1983，105(5):312－318.

［5］ Lewis D R，Norris J D.Techniques for thin－skin heat－transfer measurements for surface roughness at AEDC tunnel 9［C］∥26th AIAA Aerodynamic Measurement Technnology and Ground Testing Conference，Seattle，Washington，2008.

［6］ Shultz D L，Jones T V.Heat－transfer measurements in short－duration hypersonic facilities［R］.AGARD－AG－165，1973.

［7］ 胡 芃，陳則韶.量熱技術(shù)和熱物性測定［M］.2 版.合肥:中國科學(xué)技術(shù)大學(xué)出版社，2009.

［8］ Carslaw H S，Jaeger J C.Conduction of heat in solids［M］.2nd ed.Oxford:Clarendon Press，1986.

［9］ Incropera F P，DeWitt D P，Bergman T L，et al.New York:Fundamentals of heat and mass transfer［M］.6th ed.New York:John Wiley＆Sons Inc，2006.

［10］劉初平.氣動熱與熱防護(hù)試驗(yàn)熱流測量［M］.北京:國防工業(yè)出版社，2013.