李西園，侯雅琴，高慶華，張麗娜，王 晶,3

（1. 北京衛(wèi)星環(huán)境工程研究所，北京 100094； 2. 北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院，北京 100191；3. 北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點實驗室，北京 100094）

0 引言

火星表面大氣壓力約為700 Pa，主要成分為CO2，溫度-120～20 ℃，風(fēng)速 0～15 m/s。極低氣壓環(huán)境下存在的對流換熱使火星表面巡視器面臨著與地球和一般航天器軌道空間不同的環(huán)境。為了在極端環(huán)境下達到熱模型修正、早期故障篩除、性能測試等目的，火星巡視器的地面熱試驗一般需要極端環(huán)境疊加風(fēng)速模擬[1-2]，即在中型空間環(huán)境模擬設(shè)備內(nèi)部，通過調(diào)溫?zé)岢聊M冷黑背景，調(diào)溫冷板模擬巡視器的地表傳熱，壓力控制系統(tǒng)模擬1400 Pa左右壓力環(huán)境，通過風(fēng)機模擬最大15 m/s的風(fēng)速，并需要對試驗中的風(fēng)速進行實時測量。

目前工業(yè)上常用的風(fēng)速測量手段主要包括動壓式、熱式和超聲式等。動壓式風(fēng)速測量主要包括風(fēng)杯式、皮托管和五孔探針等形式，其中皮托管和五孔探針已經(jīng)被廣泛應(yīng)用于低氣壓風(fēng)速測量領(lǐng)域。在火星表面壓力下，0～15 m/s風(fēng)速帶來的動壓僅有0.01～10 Pa，盡管現(xiàn)有微壓傳感器產(chǎn)品的分辨率可以達到10-3Pa，但其體積、質(zhì)量往往較大，且無法在低氣壓下使用，需要使用管路在容器外進行測量，而多點風(fēng)速、風(fēng)向的測量將會使系統(tǒng)變得非常復(fù)雜[3]。熱式風(fēng)速測量包括熱線、熱球和熱膜等形式，其工作原理是可根據(jù)不同風(fēng)速下傳感器探頭換熱量的不同計算風(fēng)速[4]。在“Netlander”“海盜號”“小獵犬二號”等火星探測器上均搭載有熱式風(fēng)速儀，部分還具有通過尾跡線測量風(fēng)向的能力[5-6]。熱式風(fēng)速儀結(jié)構(gòu)簡單緊湊、質(zhì)量小、功耗低、不需要移動部件，但其缺點也很明顯：1）傳感器對標定的依賴性大，無法通過分析手段獲得信號?風(fēng)速的關(guān)聯(lián)式；2）隨著時間推移，傳感器阻值等特性變化會引起標定漂移；3）空氣溫度變化會導(dǎo)致信號的變化，故傳感器的測量精度非常依賴于氣體溫度的測量。目前火星探測器使用的熱式風(fēng)速儀均為定制產(chǎn)品。超聲風(fēng)速測量利用了聲波信號在介質(zhì)中傳播速度僅和氣體成分及溫度相關(guān)，而和壓力無關(guān)的特性。相比熱式風(fēng)速儀，超聲風(fēng)速儀不需要每個都標定，其主要缺點在于超聲波在低氣壓下衰減較大，現(xiàn)有產(chǎn)品往往很難從噪聲信號中篩取超聲換能器發(fā)射的信號，需要設(shè)計全新的換能器[7-8]，同時設(shè)備體積較大，難以滿足多點風(fēng)速測量的要求。

綜上所述，目前在極低氣壓下，風(fēng)速測量的主要途徑為定制型風(fēng)速傳感器。本文主要針對在低氣壓下使用現(xiàn)有工業(yè)用恒熱流熱球風(fēng)速傳感器的適用性，采用無量綱數(shù)分析方法建立熱球風(fēng)速傳感器表面的傳熱模型，分析傳感器在不同壓力下的響應(yīng)特性，并搭建試驗系統(tǒng)進行測試驗證。

1 熱球風(fēng)速儀對流換熱模型

1.1 仿真模型

熱球風(fēng)速探頭的核心為陶瓷封裝的微小球體，直徑約為0.6 mm，球體內(nèi)部包括熱電偶的熱端和加熱絲，外部包括熱電偶的冷端，圖1所示為本文使用的熱球風(fēng)速探頭的結(jié)構(gòu)。

圖1 熱球風(fēng)速探頭結(jié)構(gòu)Fig. 1 Structure of hot-bulb anemometer probe

風(fēng)速測量中，加熱絲上的恒定電流產(chǎn)生恒定的功率，此時，熱電偶冷端溫度可以認為是流體溫度，而熱端溫度為熱球溫度。當(dāng)不同速度流體流過熱球表面時，對流換熱系數(shù)的不同導(dǎo)致了球表溫度的不同，繼而在熱電偶冷、熱端之間形成熱電勢差；通過毫伏變送器可以對熱電勢差進行測量、傳輸，根據(jù)傳感器的標定數(shù)據(jù)可換算出當(dāng)前的風(fēng)速。

1.2 無量綱數(shù)模型建立

傳感器敏感頭與周圍環(huán)境的換熱可以分為輻射換熱、對流換熱及通過引線導(dǎo)熱3種，當(dāng)忽略沿細線上的導(dǎo)熱時，在熱平衡狀態(tài)下探頭的加熱量等于探頭的對外換熱量之和，

式中：Qpower為傳感器敏感頭熱耗，W；Qconv為通過對流散失的熱量，W；Qrad為通過輻射散失的熱量，W；h為對流換熱系數(shù)，W/(m2·K)；S為熱球表面積，m2；ε為熱球表面發(fā)射率；σ為斯忒藩?玻耳茲曼常量，5.67×10-8W/(m2·K4)；Ts、Te分別為熱球表面溫度和來流氣體溫度?？紤]到敏感頭實際與設(shè)計尺寸的偏差及引線上損失的熱耗時，式(1)可表達為

式中：C1為敏感頭上電阻在整個加熱回路電阻中的占比；C2為敏感頭實際表面積與設(shè)計表面積的比值。式(2)可以繼續(xù)簡化為

其中C為常數(shù)，包括了C1和C2兩項常數(shù)的影響，通過數(shù)據(jù)擬合得到。

式(1)～(3)中的對流換熱系數(shù)h可由努賽爾數(shù)Nu表示，

式中：λ為傳熱系數(shù)，W/(m·K)，隨氣體壓力變化不大；l為特征長度，m；Nu一般可以表示為雷諾數(shù)Re的函數(shù)。球體表面流動的雷諾數(shù)可以表示為

式中：v為氣體流速，m/s；l為特征長度，m，取0.000 6 m；υ為氣體的運動黏度，m2/s。υ隨氣體密度ρ變化，故將其轉(zhuǎn)換為動力黏度μ和ρ的比值，其中，動力黏度μ是不隨ρ變化的，氣體密度ρ可由理想氣體狀態(tài)方程，通過氣體壓力、溫度計算得到。根據(jù)計算，常壓、0～15 m/s風(fēng)速下，Re約為20～500；而在 1400 Pa環(huán)境下，Re僅為 0.1～2。針對球體表面強制對流Nu?Re經(jīng)驗關(guān)聯(lián)式，Kramers、Whitaker、Yuge、Vilet、Raithby 等學(xué)者通過試驗，給出了適用于不同范圍的表達式[9]，但大多適用于中高雷諾數(shù)，在低雷諾數(shù)（0.1～10）范圍偏差較大，不同經(jīng)驗公式間偏差可達20%[10]。根據(jù)統(tǒng)計，Whitaker關(guān)聯(lián)式

在低雷諾數(shù)時與試驗結(jié)果最為接近，其中Pr為流體普朗特數(shù)。

若式(4)～(6)中的定性溫度均選取熱球表面溫度Ts與來流氣體溫度Te的平均值，傳熱系數(shù)λ、動力黏度μ、普朗特數(shù)Pr均表示為定性溫度的函數(shù)，則當(dāng)來流氣體的溫度、速度已知時，可通過式(4)～(6)計算敏感頭表面對流換熱系數(shù)h。再將各參數(shù)代入式(3)后可求解非線性方程獲得熱球風(fēng)速傳感器的表面溫度。當(dāng)熱電偶的輸出信號ΔU與溫差ΔT滿足ΔU=k×ΔT關(guān)系時，即可由ΔT計算得到熱電偶輸出信號ΔU。本文試驗所用熱電偶的k約為0.041 mV/℃。

1.3 自然對流與強制對流強度分析

熱式風(fēng)速儀傳感器探頭表面存在自然對流和強制對流。對于恒熱流風(fēng)速探頭而言，自然對流隨著探頭的表面溫度而變化，會給測量帶來額外的誤差，故測量時應(yīng)當(dāng)確保探頭表面的強制對流占據(jù)主導(dǎo)地位?？筛鶕?jù)無量綱數(shù)Gr/Re2來判斷強制對流或自然對流的主導(dǎo)地位：Gr/Re2＜0.01時強制對流占據(jù)主導(dǎo)地位，Gr/Re2＞10時自然對流占據(jù)主導(dǎo)地位，Gr/Re2在0.01～10之間時則需要同時考慮這2種對流因素的影響。

式中：Gr為格拉曉夫數(shù)；g為重力加速度，m/s2；ΔT為熱球與環(huán)境氣體間的溫差，K；av為體脹系數(shù)，K-1，其表達式為

其中：ρ為氣體密度，kg/m3；P為壓力，Pa；T為定性溫度，K；R為摩爾氣體常數(shù)?？梢婓w脹系數(shù)與壓力無關(guān)，即衡量自然對流與強制對流模式占比的無量綱數(shù)僅和重力加速度、定性溫度、來流氣體與熱球表面的溫差、特征長度、來流速度有關(guān)，其中重力加速度、特征長度、來流速度為給定值，溫差可通過1.2節(jié)方法求解獲得，定性溫度T取來流氣體與熱球表面溫度的平均值。

2 試驗驗證

2.1 試驗系統(tǒng)及方法

為了對無量綱數(shù)模型進行驗證，需要在地面建立低氣壓風(fēng)速模擬手段。目前對于1400 Pa左右低氣壓下的風(fēng)速模擬有旋轉(zhuǎn)、閉式風(fēng)洞和開式風(fēng)洞3種基本方法[3]。其中旋轉(zhuǎn)法的基本原理如圖2所示，旋轉(zhuǎn)的風(fēng)速儀穿過流體，流體介質(zhì)不運動，風(fēng)速儀與流體介質(zhì)之間存在相對運動。一般是把懸臂和轉(zhuǎn)臺放置于空間環(huán)境模擬容器中，容器中充入與使用環(huán)境具有相同壓力、溫度的氣體。旋轉(zhuǎn)風(fēng)速模擬系統(tǒng)具有成本低，結(jié)構(gòu)簡單，容器內(nèi)壓力、氣體成分穩(wěn)定的優(yōu)點。

圖2 基于旋轉(zhuǎn)的風(fēng)速模擬系統(tǒng)組成Fig. 2 Wind speed simulation system based on rotation

本文基于中型空間環(huán)境模擬設(shè)備搭建了旋轉(zhuǎn)風(fēng)速模擬系統(tǒng)，見圖3。熱球風(fēng)速探頭固定于懸臂端部，已知懸臂半徑和懸臂轉(zhuǎn)速時可知當(dāng)前風(fēng)速，懸臂上安裝的12位精度毫伏信號采集器可直接工作于真空環(huán)境下，采集分辨率約為0.01 mV，轉(zhuǎn)換為485信號后由轉(zhuǎn)臺上的電刷滑環(huán)導(dǎo)出，并傳送至容器外的計算機存儲。傳感器的恒定熱流施加使用了容器外的程控電源，功率恒為0.08 W。

圖3 風(fēng)速模擬系統(tǒng)實物Fig. 3 Photo of the wind speed simulation system

2.2 試驗結(jié)果與分析

分別對常壓、4×104Pa、1400 Pa下的風(fēng)速傳感探頭輸出進行了測試。在試驗中背景溫度與流體溫度相同，均為室溫。

2.2.1 常壓下

圖4和圖5分別為常壓下輸出信號與無量綱數(shù)方法預(yù)測的信號對比（其中擬合常數(shù)C≈0.067）及相應(yīng)風(fēng)速下的Gr/Re2，其中受限于懸臂在常壓下的空氣阻力，僅模擬了0～11 m/s風(fēng)速下的傳感器輸出。由圖可見，在常壓下，恒熱流式熱球風(fēng)速儀的無量綱數(shù)模型與實際輸出信號的吻合度非常好，偏差在0.05 mV以內(nèi)。同時，無量綱數(shù)Gr/Re2均＜0.01，即熱球風(fēng)速儀表面強制對流占主導(dǎo)地位。

圖4 常壓下探頭輸出與無量綱分析對比Fig. 4 Comparison between probe output and dimensional analysis under ambient pressure

圖5 常壓不同風(fēng)速下的Gr/Re2Fig. 5 Gr/Re2 at different wind speeds under ambient pressure

2.2.2 低壓下

繪制4×104Pa、1400 Pa下傳感器輸出與無量綱數(shù)模型分析結(jié)果的對比（見圖6）。在4×104Pa下，二者數(shù)值較為接近，但在1400 Pa壓力下，分析值與試驗輸出存在較大的偏差，這與文獻[11]中針對熱線風(fēng)速儀的試驗結(jié)果是一致的，即在極低氣壓下，Nu?Re關(guān)聯(lián)式會隨壓力的降低產(chǎn)生較大的偏差。因此，對于低氣壓下的熱式風(fēng)速儀的適用性研究，需要結(jié)合試驗數(shù)據(jù)進行修正。

圖6 4×104 Pa與1400 Pa下探頭輸出與無量綱分析對比Fig. 6 Comparison of probe output and Gr/Re2 analysis under pressures of 4×104 Pa and 1400 Pa

同時，由圖6可見，在1400 Pa下，傳感器在不同風(fēng)速下的靈敏度約為0.1～0.2 mV/(m·s-1)；在低氣壓下進行風(fēng)速測量試驗時，毫伏信號的數(shù)據(jù)采集器亦應(yīng)當(dāng)滿足相應(yīng)的精度要求。此外，由于輻射換熱的背景溫度相關(guān)性，標定環(huán)境應(yīng)該與使用環(huán)境具有接近的輻射換熱背景，或在傳感器表面安裝高反射比的風(fēng)道以降低輻射換熱的影響。

根據(jù)試驗結(jié)果計算熱球風(fēng)速探頭表面的Gr/Re2，繪制1400 Pa下的Gr/Re2曲線如圖7所示。當(dāng)風(fēng)速達到1 m/s以上時，無量綱數(shù)Gr/Re2恒＜0.01，即強制對流在總對流換熱中仍占據(jù)主導(dǎo)地位，在這個范圍內(nèi)自然對流不會給測量帶來額外的誤差。

圖7 1400 Pa，不同風(fēng)速下的Gr/Re2Fig. 7 Gr/Re2 at different wind speeds under the pressure of 1400 Pa

3 結(jié)論

針對低氣壓下的風(fēng)速測量問題，本文使用無量綱數(shù)分析方法建立了恒熱流熱球風(fēng)速探頭的風(fēng)速?輸出信號模型，基于中型空間環(huán)境模擬試驗設(shè)備搭建了旋轉(zhuǎn)式低氣壓風(fēng)速標定系統(tǒng)，分別在常壓、4×104Pa、1400 Pa的環(huán)境下進行了測試。得到了如下結(jié)論：

1）常壓下的試驗表明，本文建立的無量綱數(shù)分析模型與試驗測試取得了一致的結(jié)果，證明了該模型對于評估傳感器輸出的有效性。

2）4×104Pa下的試驗表明，無量綱數(shù)模型與試驗結(jié)果偏差主要集中在低風(fēng)速范圍，在風(fēng)速為2 m/s以上時，其偏差在10%以內(nèi)，該模型可用于4×104Pa以上環(huán)境的熱球風(fēng)速傳感器輸出信號預(yù)測。

3）1400 Pa下的試驗表明，無量綱數(shù)分析模型與試驗結(jié)果存在25%左右的偏差，對于極低氣壓下的熱球風(fēng)速儀適用性分析，需要結(jié)合試驗手段進行修正；同時，0.6 mm探頭表面仍以強制對流換熱為主，探頭的靈敏度為0.1～0.2 mV/(m·s-1)，使用更高精度的采集儀器可以滿足1400 Pa下的風(fēng)速測量需求。