耿海，孫明明，羅俊華，劉家濤

蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)試驗(yàn)室 蘭州 730000

離子推力器基于其高比沖、長(zhǎng)壽命及高效率等優(yōu)點(diǎn)，在航天器在軌推進(jìn)領(lǐng)域得到了廣泛應(yīng)用[1-2]。

離子推力器在軌飛行前，需在地面進(jìn)行一系列測(cè)試及考核。由于其只能在真空環(huán)境下工作，因此必須建造地面真空環(huán)境設(shè)備以開展相關(guān)的測(cè)試和驗(yàn)證。電推進(jìn)真空環(huán)境設(shè)備為封閉空間，且通過(guò)各類真空泵維持設(shè)備內(nèi)真空度始終處于較高水平(一般要求真空度≤10-3Pa)。同時(shí)，由于推力器在工作時(shí)噴射出的高速、高溫且呈電中性狀態(tài)的氙氣羽流，若直接轟擊至設(shè)備壁面，會(huì)造成設(shè)備壁面出現(xiàn)離子刻蝕現(xiàn)象。因此，在電推進(jìn)真空測(cè)試設(shè)備內(nèi)須配置濺射靶以降低羽流中高能離子造成的壁面刻蝕影響，且濺射靶需設(shè)計(jì)為至少對(duì)推力器羽流具有一次光學(xué)屏蔽能力[3]。2019年中國(guó)空間技術(shù)研究院KM4真空環(huán)境設(shè)備內(nèi)開展30 cm離子推力器聯(lián)試試驗(yàn)，從結(jié)果來(lái)看，推力器在穩(wěn)定運(yùn)行一段時(shí)間后，濺射靶和柵極溫度均快速升高，且柵極截獲電流增大，與之對(duì)應(yīng)的則是推力器的工作性能(基于電參數(shù)得到的推力、比沖等)發(fā)生了衰退。分析其原因認(rèn)為，由于濺射靶在推力器羽流影響下溫度快速升高，并通過(guò)熱輻射引起柵極表面溫度升高，近一步造成柵極間距變化，從而引起束流聚焦?fàn)顟B(tài)以及整體工作性能的變化。因此，研究濺射靶對(duì)推力器影響，對(duì)保證推力器地面測(cè)試試驗(yàn)的順利開展以及分析在軌和地面條件下推力器性能不一致性均具有重要意義。

從目前國(guó)內(nèi)外的研究情況來(lái)看，研究方向主要集中在設(shè)備內(nèi)部的推力器羽流分布特性，對(duì)于設(shè)備內(nèi)引入濺射靶后的羽流特性，及濺射靶對(duì)推力器的影響研究甚少。國(guó)外前期對(duì)設(shè)備內(nèi)推力器羽流特性進(jìn)行了典型研究。Goebel[4]指出，封閉真空設(shè)備內(nèi)的中性氣體主要由3部分組成，即推力器羽流、中和器逸出氣體及設(shè)備壁面出氣。其中推力器羽流是主要組成部分，并且羽流可以近似等效為從與柵面尺寸相當(dāng)?shù)膱A形出口，并以與出口垂直的方向噴射。Soulas[5]對(duì)NEXT離子推力器在設(shè)備內(nèi)的中性羽流分布進(jìn)行了模擬，結(jié)果顯示真空設(shè)備下游區(qū)域的背景氣體壓強(qiáng)主要由推力器噴射的中性羽流與設(shè)備壁面發(fā)生的碰撞而形成。Crofton[6]測(cè)量了真空設(shè)備內(nèi)的T5離子推力器的噴射羽流，發(fā)現(xiàn)中和器羽流相比推力器羽流可忽略。Sengupta[7-8]測(cè)量了NSTAR離子推力器在地面測(cè)試設(shè)備內(nèi)工作時(shí)的中性氣體密度和等離子體參數(shù)，測(cè)量結(jié)果顯示柵極出口附近具有較高的氣體密度，而沿推力器中心線向設(shè)備底部(濺射靶)方向，氣體密度呈現(xiàn)單調(diào)遞減。Katz[9]模擬了NSTAR離子推力器在地面測(cè)試設(shè)備內(nèi)部的氣體密度分布，計(jì)算結(jié)果顯示，離子推力器工作時(shí)其內(nèi)部壓強(qiáng)在1×10-2Pa量級(jí)或者更低，并且采用Knudesen流(即分子流)模型可以描述推力器放電室內(nèi)的中性氣體密度。除此之外，中性原子與設(shè)備壁面的碰撞過(guò)程呈余弦分布狀態(tài)，碰撞將使得羽流氣體溫度與壁面溫度達(dá)到一致。Walker[10]對(duì)50 kW霍爾推力器在真空系統(tǒng)內(nèi)的羽流特性以及背景壓強(qiáng)分布進(jìn)行了模擬。結(jié)果顯示，推力器出口(柵極面)處的壓強(qiáng)較大，而隨著羽流在真空設(shè)備內(nèi)的擴(kuò)散過(guò)程，在距離推力器出口下游約2 m左右的位置，推力器羽流壓強(qiáng)逐漸降低至稍高于設(shè)備背景壓強(qiáng)，在靠近設(shè)備壁面處，羽流壓強(qiáng)基本與設(shè)備背景壓強(qiáng)一致。Korkut[11]對(duì)真空設(shè)備內(nèi)離子推力器中性羽流以及氣體反流過(guò)程進(jìn)行了三維模擬，模擬結(jié)果顯示反流氣體增強(qiáng)了推力器內(nèi)部的放電效應(yīng)，并且此時(shí)產(chǎn)生的CEX離子密度相比沒有反流氣體影響時(shí)要高出一個(gè)量級(jí)。除此之外模擬結(jié)果顯示，推力器羽流中的中性原子速度基本上是固定常數(shù)。

綜上，目前關(guān)于設(shè)備內(nèi)濺射靶對(duì)羽流分布以及濺射靶對(duì)推力器的影響研究較少，且為更深入分析前述2019年KM4真空設(shè)備內(nèi)的推力器性能衰退原因，本文計(jì)劃采用有限元分析和試驗(yàn)驗(yàn)證的方法，首先對(duì)地面真空設(shè)備內(nèi)的離子推力器羽流分布進(jìn)行模擬，獲得由于羽流對(duì)濺射靶的能量沉積效應(yīng)所造成的濺射靶溫度變化，在此基礎(chǔ)上進(jìn)一步分析濺射靶溫度升高對(duì)離子推力器柵極溫度以及熱形變位移所造成的影響，最后開展試驗(yàn)驗(yàn)證，并進(jìn)行誤差分析。

1 推力器羽流模型

獲得推力器在真空設(shè)備內(nèi)的羽流特性，是計(jì)算濺射靶能量沉積的前提，因此首先需建立起推力器羽流模型。圖1(a)為針對(duì)30 cm離子推力器性能測(cè)試及壽命考核而研制的專用地面測(cè)試設(shè)備，設(shè)備分為主艙和副艙，其中副艙尺寸為Φ2 m×2 m，主艙尺寸為Φ4.5 m×10 m，推力器被安裝于副艙，推力器軸線與設(shè)備軸線方向一致且柵面近似與副艙截面平行。設(shè)備兩側(cè)等間距總共安裝有6臺(tái)低溫泵，低溫泵總抽速約為260 kL/s，可在推力器處于額定工作條件下提供1×10-4Pa的真空度水平。設(shè)備主艙底部安裝有與主艙截面大小近似相等的采用TC-4鈦合金材料制成的濺射靶，濺射靶為百葉窗型結(jié)構(gòu)，具有18個(gè)等間距且均平行葉片，葉片傾角設(shè)計(jì)為45°以確保推力器羽流碰撞至濺射靶表面不會(huì)直接返回至推力器。其次，考慮到推力器羽流帶有較高的熱量和動(dòng)能，因此濺射靶葉片背面焊接有液氮盤管以降低推力器羽流造成的影響。在測(cè)試期間，設(shè)備內(nèi)部安裝有多種傳感器用于監(jiān)控試驗(yàn)過(guò)程中的溫度、真空度及束流大小的變化。圖1(b)為根據(jù)電推進(jìn)專用地面測(cè)試設(shè)備真實(shí)尺寸所建立的推力器羽流Fluent分析模型，由于該模型主要用于分析設(shè)備內(nèi)部推力器羽流分布，因此在模型建立過(guò)程中對(duì)結(jié)構(gòu)進(jìn)行了相應(yīng)簡(jiǎn)化，僅保留了關(guān)鍵的氣體入口、出口和濺射靶，推力器則以簡(jiǎn)單的柱段-錐段結(jié)構(gòu)代替，而對(duì)于濺射靶對(duì)離子推力器的影響，則必須新建完整的推力器結(jié)構(gòu)開展相應(yīng)分析。其次，在建立羽流模型中，對(duì)模型不同的區(qū)域選擇了不同的網(wǎng)格尺寸，尤其是對(duì)于氣體入口和出口、濺射靶等，通過(guò)增加節(jié)點(diǎn)數(shù)和控制單元尺寸，以最大程度保證計(jì)算精度。而對(duì)于羽流擴(kuò)散區(qū)域則可以適當(dāng)選擇較大的網(wǎng)格尺寸，以降低計(jì)算迭代時(shí)長(zhǎng)。其中，推力器出口面的最小單元尺寸為1.69×10-3m，是推力器口徑(約0.34 m)的1/200，出口面的最小單元尺寸為2.21×10-3m，約是低溫泵口徑(0.45 m)的1/200。

圖1 地面試驗(yàn)設(shè)備及有限元分析模型Fig.1 Ground test facility and FEM model of vacuum chamber

對(duì)于圖1(b)的邊界條件設(shè)置，首先需判斷設(shè)備內(nèi)推力器羽流的流體狀態(tài)，因此需計(jì)算設(shè)備內(nèi)Xe原子平均自由程和Knudson數(shù)：

(1)

式中：λ為Xe原子平均自由程；ra為Xe原子直徑(3.6×10-10m)；k為波爾茲曼常數(shù)；d為設(shè)備直徑；P為設(shè)備內(nèi)部平均壓強(qiáng)(取30 cm離子推力器實(shí)際工作中的最惡劣真空度5×10-3Pa，實(shí)際真空度要小于該值，按此計(jì)算的自由程λ為最小值)；T為設(shè)備內(nèi)Xe氣體溫度(取為600 K，近似為放電室溫度)?；谑?1)計(jì)算得到Knudson數(shù)大于10，設(shè)備內(nèi)氣體處于分子流狀態(tài)，且由于Xe原子平均自由程較大，因此推力器羽流基本為定向運(yùn)動(dòng)。基于計(jì)算結(jié)果和國(guó)外結(jié)論[9]，流體模型采用不可壓縮的Knudsen流，且分析過(guò)程同時(shí)考慮了能量交換過(guò)程。

其他邊界條件設(shè)置包括，推力器柵極面被設(shè)置為壓強(qiáng)入口(pressure inlet)，并且推力器羽流近似為從與柵面等大小的圓盤垂直噴射[4]。根據(jù)前期試驗(yàn)測(cè)量結(jié)果[12]，柵面壓強(qiáng)設(shè)置為0.01 Pa并假設(shè)柵面各處壓強(qiáng)均一致。低溫泵入口面設(shè)置為速度出口(velocity outlet)，由于單臺(tái)低溫泵口徑為0.45 m，對(duì)Xe氣抽速為4 600 L/s，計(jì)算得到圖1(b)所示模型的出口速度為28.7 m/s。設(shè)備壁面設(shè)置為出氣邊界(outflow boundary)以模擬設(shè)備材料氣體釋放過(guò)程(根據(jù)辜學(xué)茂等[13]的研究結(jié)論，設(shè)置為0.01 PaL/s)。同時(shí)，壓強(qiáng)入口的溫度設(shè)置為450℃(與柵極實(shí)測(cè)溫度相等[14])，流體的動(dòng)能和能量交換采用二階迎風(fēng)算法(second order upwind scheme)，默認(rèn)迭代次數(shù)為200次，并且迭代次數(shù)可隨著計(jì)算結(jié)果的收斂與否進(jìn)行增加。

2 設(shè)備內(nèi)羽流模擬結(jié)果

根據(jù)模型及邊界條件設(shè)置，模擬得到設(shè)備內(nèi)部推力器羽流特性如圖2所示。需注意的是，圖2(a)所示結(jié)果為僅考慮了濺射靶表面溫度對(duì)羽流溫度影響，未考慮羽流溫度對(duì)濺射靶溫度的作用，主要由于在建模過(guò)程中濺射靶僅以表面形式給出以模擬流-固耦合過(guò)程，未建立實(shí)體模型。濺射靶對(duì)推力器和空間的熱輻射以及自身的溫度傳導(dǎo)效應(yīng)需在下一小節(jié)中進(jìn)行討論。

從圖2(a)所示結(jié)果來(lái)看，推力器羽流溫度在真空艙軸線方向上呈現(xiàn)明顯的一致性。這主要由于真空艙內(nèi)氣體密度較低，羽流呈定向分子流狀態(tài)，羽流內(nèi)部幾乎不會(huì)發(fā)生粒子間的相互碰撞，即羽流在擴(kuò)散過(guò)程中幾乎沒有能量損失，而是僅與設(shè)備壁面發(fā)生碰撞，因此呈現(xiàn)圖2(a)所示結(jié)果。圖2(b)結(jié)果顯示，真空艙內(nèi)大部分的氣體壓強(qiáng)在2×10-3～6×10-3Pa之間，但濺射靶后部的氣體壓強(qiáng)要明顯高于真空艙中段區(qū)域的壓強(qiáng)，這是由于濺射靶阻擋了大部分反流氣體，導(dǎo)致濺射靶后部的氣體密度遠(yuǎn)高于處于分子流狀態(tài)的羽流擴(kuò)散區(qū)域(即設(shè)備中段)，并且由于濺射靶的45°傾角結(jié)構(gòu)，造成濺射靶后部靠設(shè)備底面區(qū)域的氣體壓強(qiáng)要高于頂部區(qū)域。圖2給出的模擬結(jié)果一方面顯示了分子流狀態(tài)下的推力器羽流溫度以及壓強(qiáng)分布特性，另一方面分析結(jié)果將作為關(guān)鍵的參數(shù)設(shè)置以開展濺射靶和推力器的相互影響分析。

3 濺射靶與推力器模型

圖1(b)所示模型僅針對(duì)推力器羽流特性分析，但對(duì)于濺射靶與推力器的相互影響分析，僅依靠圖1(b)模型是遠(yuǎn)不足的，需進(jìn)一步建立詳細(xì)的推力器模型和濺射靶模型。為獲得推力器羽流對(duì)濺射靶溫度的影響，以及濺射靶溫度升高后對(duì)推力器本體的影響，本文建立了詳細(xì)的推力器和濺射靶模型，如圖3所示。模型采用ANSYS建立，圖3(a)為僅顯示推力器和濺射靶后的分析模型，圖3(b)為劃分網(wǎng)格后的推力器模型。同時(shí)為了與真實(shí)條件保持一致，濺射靶葉片背面建立了液氮盤管模型，用于模擬實(shí)際工作中，液氮管對(duì)濺射靶的降溫過(guò)程。

圖3 推力器及濺射靶有限元分析模型Fig.3 FEM model of ion thruster and the sputtering target

對(duì)于推力器內(nèi)部熱邊界條件的設(shè)置，本課題組前期已經(jīng)針對(duì)30 cm離子推力器的熱分析模型，模型材料屬性等效，工作時(shí)的熱邊界設(shè)置以及熱特性模擬開展過(guò)大量研究[14]，因此本文不再進(jìn)行敘述，僅對(duì)濺射靶與推力器間的熱輻射邊界進(jìn)行說(shuō)明。當(dāng)考慮濺射靶影響時(shí)，推力器自身內(nèi)部輻射關(guān)系保持不變，主要增加了推力器減速柵對(duì)濺射靶的輻射關(guān)系，其中減速柵和濺射靶的熱輻射系數(shù)分別為0.5和0.8(均為測(cè)量值)。同時(shí)，考慮到推力器側(cè)面與真空測(cè)試設(shè)備壁面存在輻射換熱，故設(shè)置推力器側(cè)面對(duì)設(shè)備壁面的輻射關(guān)系，二者發(fā)射系數(shù)分別為0.7和0.2(測(cè)量值)，設(shè)備壁面設(shè)置為恒溫-120℃(測(cè)量值)。當(dāng)不考慮濺射靶影響時(shí)，分析模型中去除濺射靶結(jié)構(gòu)，推力器減速柵直接與設(shè)備壁面發(fā)生熱交換關(guān)系，設(shè)備壁面熱輻射系數(shù)和溫度邊界同樣取為0.2和-120℃，其余輻射關(guān)系均保留不變。

獲得濺射靶溫度升高對(duì)推力器工作平衡溫度的影響，關(guān)鍵在于準(zhǔn)確模擬濺射靶在受高溫、高動(dòng)能羽流粒子沖擊時(shí)的溫度變化情況。推力器羽流可分為近場(chǎng)區(qū)和遠(yuǎn)場(chǎng)區(qū)，如圖4所示。其中近場(chǎng)區(qū)(約為2倍推力器半徑)存在強(qiáng)烈的電離效應(yīng)，以電荷交換碰撞、電離碰撞及離子復(fù)合碰撞為主，該區(qū)域具有高密度等離子體。而進(jìn)入遠(yuǎn)場(chǎng)區(qū)，粒子間碰撞以離子復(fù)合為主，且95%以上Xe離子均被還原為Xe原子，該區(qū)域的推力器羽流基本不再具有帶電粒子，可看做高速中性Xe氣流轟擊至低溫濺射靶，因此主要研究遠(yuǎn)場(chǎng)區(qū)推力器羽流對(duì)濺射靶的能量沉積效應(yīng)。

由于推力器羽流為定向分子流狀態(tài)，因此遠(yuǎn)場(chǎng)區(qū)以中性Xe原子和低溫濺射靶的碰撞為主。為計(jì)算高速Xe氣羽流轟擊至濺射靶所產(chǎn)生的能量沉積，首先假設(shè)設(shè)備內(nèi)氣體為單一氣體(Xe氣)，忽略其他如氧氣或水汽成分(從本課題組前期采用四極質(zhì)譜計(jì)對(duì)推力器工作時(shí)的設(shè)備內(nèi)氣體成分測(cè)量結(jié)果來(lái)看，氙氣分壓力比其他雜質(zhì)氣體分壓力要高出兩個(gè)量級(jí))。而對(duì)于中性Xe原子與濺射靶的碰撞，首先考慮其吸附方式。其中化學(xué)吸附主要發(fā)生在固體-氣體之間，且只有在具有化學(xué)活性的氣體與固體表面才能發(fā)生(如氧氣和金屬表面)，原子基本只能吸附單層，且吸附粒子與固體間存在電子交換。而物理吸附是分子間普遍存在的Van Der Waals力，任何表面對(duì)任何氣體均會(huì)發(fā)生物理吸附，固體表面可形成多層吸附，也可看做是凝結(jié)過(guò)程，并且當(dāng)氣體-固體間存在越高的溫差時(shí)，物理吸附更為明顯，而溫度越高，物理吸附幾率α越低。由此可見，由于中性Xe原子屬于惰性氣體，且Xe原子與濺射靶(-196℃)存在巨大溫差，因此濺射靶對(duì)Xe原子以物理吸附為主。根據(jù)相關(guān)資料[15]，當(dāng)溫度從170 K變化至70 K時(shí)，鈦材料對(duì)Xe原子的物理吸附幾率在0.86～1之間變化，根據(jù)圖2(a)得到的羽流溫度模擬結(jié)果，本文將濺射靶對(duì)Xe原子的吸附幾率設(shè)置為0.9。

圖4 真空艙內(nèi)中性羽流引出過(guò)程Fig.4 Extraction process of neutral plume in vacuum chamber

由于遠(yuǎn)場(chǎng)區(qū)的Xe原子具有能量主要為動(dòng)能(Xe原子溫度約為600 K，熱能相比動(dòng)能基本可忽略)，其速度可表述為v=Ispg，因此單個(gè)Xe原子具有的動(dòng)能Wa為：

Wa=mav2/2=ma(Ispg)2/2

(2)

式中：Isp為推力器比沖；g和ma分別為重力加速度以及Xe原子質(zhì)量。單位時(shí)間碰撞至濺射靶單位面積的Xe原子數(shù)量N為：

(3)

式中：pf和Tf分別為遠(yuǎn)場(chǎng)區(qū)靠近濺射靶區(qū)域的氣體壓強(qiáng)和溫度；μ為氙原子相對(duì)質(zhì)量，根據(jù)圖2分析結(jié)果，將pf和Tf設(shè)置為6×10-3Pa及600 K，聯(lián)立式(2)(3)得到Wa=1.34×10-16J，N=4.7×1015/(cm2·s)。由于推力器羽流與濺射靶的碰撞面積近似為圓形[3]，根據(jù)推力器束流角度關(guān)系(束流半角為15°)，可得到濺射面積S為：

S=π(Ltan15°)2

(4)

式中：L為推力器柵面至濺射靶的距離。由于推力器羽流分布滿足高斯分布特性[16]，但隨著離出口面的距離增加，高斯分布呈現(xiàn)平緩特征，因此基本可認(rèn)為在L>1.5 m的羽流密度屬于均勻分布。基于式(4)得到濺射靶的羽流碰撞面積為5 075 cm2，因此濺射靶的熱邊界條件為：

Q=αNWaS

(5)

式中：Q為羽流沉積能量，單位為W；吸附幾率α為0.9。得到濺射靶熱流邊界(heat flow)為2 876 W，考慮施加區(qū)域面積后，可得到推力器羽流沉積在濺射靶上的熱流密度(heat flux)條件，如圖5(a)所示，同時(shí)考慮到濺射靶自身輻射以及其上安裝有液氮盤管冷卻系統(tǒng)，因此濺射靶設(shè)置對(duì)空間輻射以及溫度邊界(見圖5(b))。

圖5 濺射靶熱通量及冷卻管溫度設(shè)置Fig.5 Heat flux of the sputtering target and temperature setting of cooled tube

4 濺射靶對(duì)柵極的影響

根據(jù)推力器羽流沉積在濺射靶上的熱流密度以及其他邊界條件，基于ANSYS模擬得到濺射靶熱平衡后的溫度分布如圖6(a)所示。同時(shí)，考慮到濺射靶溫度升高后會(huì)對(duì)柵極帶來(lái)一定熱輻射影響，在模型中設(shè)置了濺射靶表面對(duì)推力器柵面的空間熱輻射關(guān)系。其中，濺射靶表面發(fā)射率取為0.8，柵極表面為0.5(Mo材料表面輕微氧化)，給出溫度平衡后的加速柵溫度如圖6(b)所示。

圖6 濺射靶與加速柵溫度分布Fig.6 Temperature distribution of the sputtering target and the accelerator grid

從圖6(a)所示結(jié)果來(lái)看，濺射靶在推力器羽流能量沉積的影響下，中心區(qū)域的溫度明顯高于其他區(qū)域，由于濺射靶葉片自身具有較高的輻射系數(shù)且葉片背部均焊接有液氮盤管，因此中心區(qū)域向周圍的熱傳導(dǎo)過(guò)程受明顯制約，高溫區(qū)域幾乎只限制在熱流密度加載區(qū)域。圖6(b)所示加速柵中心溫度為352℃，相比前期在沒有濺射靶影響下模擬得到的326℃的分析結(jié)果[14]，溫度提高了26℃，同時(shí)，屏柵中心溫度為440℃(限于篇幅未給出)，相比忽略濺射靶影響時(shí)模擬得到的407℃提高了33℃[14]。由此可見，濺射靶對(duì)柵極組件平衡溫度具有一定影響，但對(duì)于柵極來(lái)說(shuō)，關(guān)注點(diǎn)主要在于熱態(tài)下的柵極間距變化，因此基于濺射靶影響下的柵極溫度分析結(jié)果，對(duì)柵極間距變化進(jìn)一步開展分析。

屏柵和加速柵的初始冷態(tài)間距為0.85～0.95 mm，圖7給出了無(wú)濺射靶影響以及濺射靶影響下的推力器屏柵和加速柵中心的溫差變化和柵極中心的相對(duì)位移，即柵極熱態(tài)間距的縮小量。圖7(a)所示結(jié)果表明，考慮濺射靶的影響后，屏柵和加速柵相比無(wú)濺射靶情況下的最高溫度有所增大。從7(b)的結(jié)果來(lái)看，當(dāng)考慮濺射靶的影響后，屏柵和加速柵的間距縮小量從原來(lái)的0.560 mm(忽略濺射靶影響)增加到0.585 mm，若以柵極最小初始間距0.85 mm來(lái)考慮，無(wú)濺射靶和有濺射靶影響下的柵極間距分別穩(wěn)定在0.290 mm和0.265 mm。而前期試驗(yàn)表明[1]，屏柵和加速柵的安全間距為0.26 mm，因此在濺射靶的影響下，屏柵和加速柵間距已接近安全距離的臨界值。

圖7 屏柵和加速柵溫度差及熱態(tài)間距變化Fig.7 Temperature difference and hot gap variation between sc. grid and acc. grid

為了驗(yàn)證真空設(shè)備內(nèi)推力器羽流特性以及濺射靶對(duì)推力器影響的模擬結(jié)果的準(zhǔn)確性，本課題組開展了相關(guān)測(cè)試試驗(yàn)。其中，對(duì)于推力器羽流特性模擬結(jié)果的驗(yàn)證主要是測(cè)試羽流溫度和設(shè)備內(nèi)不同區(qū)域的真空度，以完成對(duì)圖2的驗(yàn)證。羽流溫度的測(cè)量采用10只鎧裝熱電偶完成，熱電偶按主艙尺寸等間距的貼裝于設(shè)備內(nèi)壁頂面，并且熱電偶測(cè)試端被人為彎折后，以距離主艙頂面30 cm的距離且垂直于頂面的方式，測(cè)量推力器羽流溫度以避免緊貼主艙壁面測(cè)量造成不準(zhǔn)確，試驗(yàn)和模擬結(jié)果的比對(duì)如圖8(a)所示。設(shè)備內(nèi)不同區(qū)域真空度的模擬結(jié)果比對(duì)是通過(guò)設(shè)備上安裝的電離真空規(guī)以完成。30 cm離子推力器地面測(cè)試設(shè)備上在其前端、中端以及靠近濺射靶位置上均安裝有真空規(guī)以監(jiān)控試驗(yàn)中的真空度變化，試驗(yàn)和模擬結(jié)果的比對(duì)如圖8(b)所示。對(duì)于濺射靶對(duì)推力器溫度影響的模擬結(jié)果驗(yàn)證，目前由于沒有直接測(cè)量柵極表面溫度的測(cè)試手段，僅對(duì)試驗(yàn)過(guò)程中的柵極安裝環(huán)部位進(jìn)行了溫度測(cè)量，模擬時(shí)設(shè)置的相應(yīng)邊界條件(如壁面溫度、發(fā)射率等)大部分均來(lái)自于測(cè)試值，以保證測(cè)試與模擬結(jié)果的可比性。試驗(yàn)結(jié)果顯示屏柵和加速柵邊緣處的熱平衡溫度分別為385℃和364℃[14]，與仿真結(jié)果得到的411℃和342℃(見圖6(b))相比誤差均為6%。

從圖8的比對(duì)結(jié)果來(lái)看，仿真結(jié)果與試驗(yàn)結(jié)果比對(duì)誤差較小，且仿真與試驗(yàn)結(jié)果的趨勢(shì)相同，均顯示出處于分子流狀態(tài)下的推力器羽流在噴射和擴(kuò)散過(guò)程中壓強(qiáng)和溫度幾乎沒有較大變化的特性，但在靠近濺射靶區(qū)域壓強(qiáng)升高而溫度快速降低的特點(diǎn)，說(shuō)明百葉窗型濺射靶對(duì)羽流起到了“隔斷”效應(yīng)，導(dǎo)致主艙底部的氣體壓強(qiáng)明顯高于前端和中段區(qū)域，并且反流氣體大部分被限制在濺射靶后部，有利于降低反流氣體對(duì)推力器造成的放電損耗升高等不利影響。圖8(a)所示濺射靶后部羽流氣體的溫度計(jì)算值與試驗(yàn)值比對(duì)誤差在50℃左右，這主要由于在仿真過(guò)程中，未考慮羽流溫度對(duì)濺射靶的影響，而實(shí)際由于羽流粒子轟擊至濺射靶造成的能量沉積效應(yīng)，羽流溫度要明顯高于模擬值。

圖8 羽流溫度及壓強(qiáng)的仿真與試驗(yàn)結(jié)果比對(duì)Fig.8 The comparison results of the plume temperature and pressure

5 結(jié)束語(yǔ)

本文采用有限元模擬方法分別對(duì)真空艙內(nèi)推力器羽流特性和濺射靶對(duì)推力器柵極的影響進(jìn)行了分析模擬，并進(jìn)行了試驗(yàn)驗(yàn)證，得到以下結(jié)論：

1)離子推力器在真空設(shè)備內(nèi)的羽流可采用定向分子流模型進(jìn)行描述，羽流在真空艙內(nèi)的擴(kuò)散過(guò)程中幾乎沒有能量損失，羽流溫度在碰撞至濺射靶之前基本不會(huì)發(fā)生變化。

2)對(duì)于30 cm離子推力器，其工作時(shí)真空艙內(nèi)大部分區(qū)域的氣體壓強(qiáng)在2×10-3～6×10-3Pa，由于濺射靶阻擋了大部分反流氣體，導(dǎo)致濺射靶后部的氣體密度遠(yuǎn)高于處于分子流狀態(tài)的羽流擴(kuò)散區(qū)域。

3)30 cm離子推力器噴射羽流對(duì)濺射靶造成一定的能量沉積效應(yīng)，導(dǎo)致濺射靶中心區(qū)域的溫度明顯高于其他區(qū)域，且高溫區(qū)域幾乎只限制在熱流密度加載區(qū)域。在濺射靶影響下，推力器加速柵和屏柵中心溫度分別為352℃和440℃，邊緣溫度分別為342℃和411℃。屏柵和加速柵的間距縮小量從原來(lái)的0.560 mm增加到0.585 mm，此時(shí)柵極間距已接近安全距離的臨界值。

4)試驗(yàn)結(jié)果顯示，屏柵和加速柵邊緣處的熱平衡溫度分別為385℃和364℃，與仿真結(jié)果得到的411℃和342℃相比誤差均為6%，且試驗(yàn)測(cè)量得到的設(shè)備內(nèi)氣體壓強(qiáng)和推力器羽流溫度與仿真結(jié)果比對(duì)誤差較小，同時(shí)試驗(yàn)與仿真結(jié)果具有相同的趨勢(shì)。濺射靶后部羽流氣體的試驗(yàn)值高于計(jì)算值約50℃，主要由于在仿真過(guò)程中，忽略了羽流粒子的能量沉積效應(yīng)。

從本文結(jié)論來(lái)看，推力器羽流對(duì)濺射靶的能量沉積效應(yīng)，導(dǎo)致濺射靶溫度快速升高是后續(xù)濺射靶設(shè)計(jì)過(guò)程中需考慮的問題。而從濺射靶對(duì)推力器本體的影響結(jié)果來(lái)看，造成了柵極溫度的小幅度升高和柵極間距的惡化。在后續(xù)試驗(yàn)過(guò)程中，應(yīng)盡量控制濺射靶溫度，以消除由于濺射靶影響帶來(lái)的柵極熱態(tài)間距變化。