熱電偶絲溫度響應(yīng)過程數(shù)值計算分析

摘要:高焓高速氣流總溫主要應(yīng)用總溫探針測量脈沖類高焓風洞對溫度響應(yīng)時間要求較高而探針溫度響應(yīng)主要是熱電偶絲的溫度響應(yīng)過程.為了分析影響熱電偶絲溫度響應(yīng)的過程針對電偶絲進行了理論分析獲得了熱電偶絲溫度響應(yīng)的理論公式并分析了熱電偶絲測溫存在的誤差.隨后利用有限差分法對熱電偶絲的溫度響應(yīng)過程進行數(shù)值計算獲得了熱電偶絲溫度響應(yīng)曲線并且結(jié)果表明:熱電偶絲直徑越小響應(yīng)時間越短;熱電偶絲結(jié)點越小其響應(yīng)時間越短;相同直徑的熱電偶絲越長其響應(yīng)越快.
0引言
　　在空氣動力學(xué)領(lǐng)域高焓高超聲速氣流的溫度測量-般指氣流總溫的測量.總溫指氣流在絕熱滯止狀態(tài)下所能達到的溫度而在實際的總溫測量過程中完全絕熱滯止是無法實現(xiàn)的,因而即使不存在使用誤差的情況下探針的指示溫度也總是小于總溫而高于氣流靜溫這就給總溫測量帶來了困難.高焓高速氣流的總溫測量必須采用侵入式方法才可以測量而采用非接觸方法只能得到氣流的靜溫因此以上方法中只有熱電偶法適用而且需要把熱電偶做成總溫探針的樣式進行高焓高超聲速氣流總溫的測量.
　　在國外總溫探針廣泛應(yīng)用于高焓高超聲速地面設(shè)備的流場校測.在AEDC/VKF研究中心C風洞的M4氣動熱風洞中采用單屏蔽熱電偶探針進行總溫測量["]該總溫探針在風洞穩(wěn)定段進行總溫測量時誤差不超過3%.AEDC高超聲速風洞(噴管馬赫數(shù)6.8.10)在1991年校測1]時校測排架上安裝有43只皮托壓力探針.25只單屏蔽總溫探針和3只流向角探針在每種馬赫數(shù)下總溫測量的誤差小于1%而沒有修正的總溫測量結(jié)果誤差為+2%.國內(nèi)對總溫測量進行了許多研究但用于高焓高超聲速風洞的測量卻不多.701所饒文成國曾在FD-20風洞上利用研制出的-種耐沖擊的快速響應(yīng)熱電偶對前室總溫進行直接測量熱電偶材料為鎳鉻-鎳硅從試驗結(jié)果來看測量誤差比較大達到10%~16%.中科院力學(xué)所王世芬[59在JF8高焓激波風洞上分別采用激波馬赫數(shù)、駐點熱流率和皮托壓力.帶屏蔽熱電偶(熱電偶材料為鎳鉻-鎳鋁)三種測試技術(shù)來測量氣流在馬赫數(shù)6.5狀態(tài)下總溫目的是比較三種測試技術(shù)及其精度同時考核自制的總溫探針特性.三種方法的試驗結(jié)果基本相符進一步證實了總溫探針在脈沖型風洞中應(yīng)用的可行性和可靠性.
　　脈沖類風洞的試驗時間比較短這要求總溫探針的溫度響應(yīng)過程要短.總溫探針的溫度響應(yīng)過程包括探針內(nèi)部流場的建立過程.滯止室溫度達到穩(wěn)定的過程以及電偶絲自身的溫度響應(yīng)過程.總溫探針放置于高速氣流中其內(nèi)部流場建立的過程耗時不到Ims可以不予考慮;滯止室的溫度達到穩(wěn)態(tài)的過程所涉及的因素較多通過理論分析較為困難,只能根據(jù)實際經(jīng)驗通過試驗來進行改進.電偶絲的溫度響應(yīng)過程對整個探針溫度響應(yīng)過程影響較大,下面通過理論分析和數(shù)值計算兩方面來研究。
1理論分析
　　假設(shè)熱電偶初始溫度和支撐物溫度相等分析其對氣流溫度(T,)響應(yīng)就是對-階躍溫度的響應(yīng)“7因。此過程的控制方程以及初始條件和邊界條件如下:
 
 
　　其響應(yīng)曲線如圖1所示.圖中溫度響應(yīng)曲線的漸近線是(Tƒ-Th)ψ.,且響應(yīng)曲線的初始斜率是(T,-T)/T其值與導(dǎo)熱誤差因子ψ。無關(guān).溫度響應(yīng)的有效時間常數(shù)為r(I-ψ).因此可以說當存在不可忽略的導(dǎo)熱傳熱時響應(yīng)的最終溫度相對于“真實”的氣流溫度T。減小了(T。-T)ψ。相應(yīng)的響應(yīng)時間常數(shù)相對于不存在導(dǎo)熱傳熱也減小了相應(yīng)的比例.在不考慮熱傳導(dǎo)誤差和輻射誤差時熱電偶的溫度響應(yīng)時間常數(shù)為T,它表示在溫度階躍響應(yīng)中電偶絲結(jié)點達到周圍氣流溫度的63.2%所需的時間.
 
2.根據(jù)以上分析的熱電偶傳熱模型如果忽略輻射傳熱可以得到電偶絲的一維非穩(wěn)態(tài)傳熱方程:
 
　　式中△x為離散網(wǎng)格長度.
　　電偶絲的離散模型如圖2所示.在計算電偶絲結(jié)點單元(M<i<M)時由于單元左右兩端的導(dǎo)熱面積不等要分別計算..數(shù)值計算中時間步長不能取太大否則會導(dǎo)致物理上的不真實解.對于本例導(dǎo)熱系數(shù)為常數(shù)方程不存在熱源項而且離散采用均勻網(wǎng)格此時時間步長要滿足不等式:
 
　　利用建立的電偶絲溫度響應(yīng)計算程序?qū)Σ煌睆健⒔Y(jié)點大小和長徑比的電偶絲進行數(shù)值計算此較它們的響應(yīng)時間.電偶絲材料選擇鉑-銠30初始溫度500K階躍至溫度1800K.響應(yīng)曲線如圖3~5.
 
3結(jié)束語
　　根據(jù)以上計算結(jié)果可知:
(1)在其他條件相同時電偶絲直徑越小響應(yīng)時間越短;
(2)不同大小的電偶絲焊接結(jié)點會導(dǎo)致溫度響應(yīng)時間的差別:結(jié)點越小其響應(yīng)時間越短這是因為結(jié)點小熱容量就小，溫升就越快;
(3)電偶絲的長度也影響溫度響應(yīng)時間:偶絲越長其響應(yīng)越快;數(shù)值計算方法為研究電偶絲的響應(yīng)過程提供了一種手段盡管數(shù)值計算中存在誤差但可以定性地分析不同的熱電偶結(jié)構(gòu)尺寸對響應(yīng)時間的影響指導(dǎo)設(shè)計高響應(yīng)速度的總溫探針.