離軌薄膜帆原子氧侵蝕的缺陷密度效應(yīng)數(shù)值模擬研究

傅宇蕾，惲衛(wèi)東，曹爭利，彭福軍，耿海峰，孫承月，吳宜勇

(1.上海宇航系統(tǒng)工程研究所上海市空間飛行器機構(gòu)重點實驗室，上海201108；2.哈爾濱工業(yè)大學(xué)材料科學(xué)與工程學(xué)院，哈爾濱150001)

1 引言

近年來，小衛(wèi)星的發(fā)展呈現(xiàn)井噴態(tài)勢，在2019年的航天發(fā)射總載荷中，立方星及重約10kg的微小衛(wèi)星占比高達37.6%。大多數(shù)小衛(wèi)星很難在短時間內(nèi)離軌并再入大氣層燒毀，從而成為長期駐留軌道、威脅其他航天器的空間碎片，給航天器的安全運行帶來了巨大威脅[1]。如圖1所示為近地軌道空間碎片分布圖[2]，對于空間碎片風(fēng)險，尤其要關(guān)注的軌道區(qū)域是800km附近軌道高度的太陽同步軌道。在這一軌道高度，碎片高度密集。另外，在正面相撞的情況下，產(chǎn)生的碰撞速度最大。因此發(fā)生在太陽同步軌道上的碰撞能量最大。這表明惡性碰撞最有可能發(fā)生在太陽同步軌道。如果碰撞的連鎖效應(yīng)開始，這一效應(yīng)極有可能發(fā)生在太陽同步軌道。根據(jù)機構(gòu)間空間碎片協(xié)調(diào)委員會 (IADC:Inter-Agency Space Debris Coordination Committee)2007年達成的國際指南要求，低軌碎片離軌時間不大于25年[3，4]。

增阻離軌薄膜帆是針對低地軌道衛(wèi)星空間碎片清除技術(shù)中可行性最高、經(jīng)濟性能最好的技術(shù)之一。離軌帆收攏小，展開面積大，在衛(wèi)星壽命結(jié)束后啟動，利用稀薄大氣阻力效應(yīng)進行降軌，大大縮短航天器在軌滯留時間，留出寶貴的軌道資源。如圖2所示為薄膜帆式離軌模塊概念圖。

離軌帆一般工作在800km以下的近地軌道，要長期暴露在高低溫、原子氧、紫外輻射、空間碎片等惡劣環(huán)境中，其中，原子氧對帆面侵蝕最為嚴(yán)重。離軌帆帆面薄膜材料的制備、帆面的折疊收攏以及發(fā)射過程中均不可避免地會產(chǎn)生孔隙、裂紋、折痕等缺陷。這些缺陷部位將成為原子氧掏蝕效應(yīng)的集中作用位置，嚴(yán)重時會導(dǎo)致帆面破碎，嚴(yán)重影響離軌帆工作性能。

圖1 近地軌道空間碎片分布圖Fig.1 Space debris distribution in low Earth orbit

圖2 薄膜帆式離軌模塊概念圖Fig.2 Conceptual map of membrane deorbit sail

國內(nèi)外學(xué)者提出了原子氧掏蝕模型[5—10]，研究了初始裂紋尺寸對侵蝕形貌的影響。在掏蝕模型中通過設(shè)定初始缺陷尺寸、原子氧能量和入射角度等參數(shù)，可以模擬材料缺陷處侵蝕形貌的演進過程。近年來，在傳統(tǒng)二維掏蝕模型的基礎(chǔ)上發(fā)展了三維模型，并且在參數(shù)選取和算法優(yōu)化方面也進行了改進[11，12]。然而，現(xiàn)有的模型仍然是針對單一缺陷進行掏蝕行為預(yù)測的，針對缺陷密度效應(yīng)和缺陷間交互影響作用的掏蝕行為研究較為欠缺。本文針對離軌帆帆面為研究對象，離軌帆帆面輕薄，對缺陷敏感，在折疊收攏以及發(fā)射過程中均不可避免地會產(chǎn)生較為密集的缺陷，缺陷密度與薄膜膜厚具有可比性，因此重點研究了帆面薄膜材料的原子氧掏蝕缺陷密度效應(yīng)。

2 離軌帆薄膜缺陷表征

離軌帆的缺陷主要在制備、帆面折疊收攏、發(fā)射過程中形成。需要對這些典型缺陷特征進行幾何測量分析，通過使用掃描電子顯微鏡精確表征缺陷形貌和尺度，為原子氧掏蝕行為預(yù)測過程提供建模依據(jù)。如圖3(a)所示，作為對照組的完整薄膜，可以看到鍍鋁薄膜表面致密完整，均勻性較好，但是，由于薄膜制備過程中不可避免地會摻雜雜質(zhì)，膜表面存在亞微米尺度不均勻分布的凸點，經(jīng)EDS面掃分析成分后，發(fā)現(xiàn)其硅含量較高，初步推斷為薄膜制備過程中的硅膠類雜質(zhì)。如圖3(b)所示為離軌帆帆面劃痕，劃痕是由于不規(guī)范的折疊工藝或制樣工藝引起的，對劃痕尺度的衡量意義不大，應(yīng)該盡量規(guī)范操作，避免鋁膜表面不必要的傷害。圖3(c)和 (d)所示分別為單次折疊后反面 (凹面)和正面 (凸面)的典型形貌。凹面鍍鋁層在壓應(yīng)力作用下形成條狀剝蝕，條間距為2～3μm量級，條狀裂紋寬度約為500nm～1μm量級。凸面鍍鋁層在拉應(yīng)力作用下形成鱗片狀剝蝕，典型鱗片尺度為2μm×2μm ～5μm ×5μm 量級， 鱗狀裂紋寬度約為500nm～1μm量級。

圖3 離軌帆帆面缺陷表征Fig.3 SEM graph of membrane defects

3 離軌帆離軌過程原子氧環(huán)境分析

離軌帆收攏小，展開大，其工作原理是利用薄膜結(jié)構(gòu)大展收比的特點，在衛(wèi)星壽命結(jié)束后啟動，利用稀薄大氣阻力效應(yīng)進行降軌，大大縮短航天器在軌滯留時間，留出寶貴的軌道資源。其一般工作于800km以下的低地軌道，其中，原子氧侵蝕對離軌帆帆面材料的影響最為嚴(yán)重，需要對離軌過程中原子氧累積通量進行分析。以金牛座離軌帆為例進行計算，其于2019年9月12日發(fā)射進入740km的太陽同步軌道。如圖4所示為金牛座離軌帆在軌運行圖。

圖4 金牛座離軌帆在軌運行照片F(xiàn)ig.4 On-orbit operation of Taurus CubeSat deorbit sail

圖5 金牛座離軌帆在軌原子氧累積通量分析Fig.5 Atomic fluence of Taurus deorbit drag sail

根據(jù)軌道分析軟件算出的離軌過程中軌道隨時間變化情況，估計離軌帆全壽命周期原子氧累積通量。估算過程中太陽活動指數(shù)取Ap(180) =150，地磁指數(shù)取F10.7(前一天) =150，F(xiàn)10.7(81天平均) =150。如圖5所示為金牛座離軌帆在軌原子氧累積通量分析圖。分析可知，金牛座離軌帆在阻力最大的飛行姿態(tài)下，經(jīng)歷的原子氧通量最小為4.67×1020atoms/cm2。在沒有任何防護鍍層的情況下，取離軌帆帆面聚合物侵蝕率 (Erosion Rate)3.4×10-24cm3/atom[13]，帆面材料在離軌過程中剝蝕厚度高達16μm。

4 原子氧蒙特卡洛侵蝕模型的建模與研究

4.1 問題簡化

在低地軌道中，原子氧與材料的掏蝕效應(yīng)，可以在統(tǒng)計意義上等效于一定數(shù)量的粒子與航天器表面材料相互作用的過程。本文采用蒙特卡洛仿真法，該方法是一種通過隨機變量進行統(tǒng)計試驗，從而解決工程實際問題的一種數(shù)值模擬算法。由地面模擬試驗和空間暴露實驗可知，原子氧對碳氫化合物材料的侵蝕主要集中在保護層的缺陷部位。如圖6所示，入射的原子氧原子以一定的概率與聚合物薄膜發(fā)射反應(yīng)。未反應(yīng)的原子發(fā)生多次鏡面反射或漫反射，直到它最終與材料發(fā)生反應(yīng)或者從缺陷中逃逸出掏蝕空腔。

仿真過程中進行了如下假設(shè):

(1)原子氧掏蝕效應(yīng)被限制在材料表面的垂直法線平面上。

(2)原子氧與鍍層不發(fā)生反應(yīng)或結(jié)合。

(3)原子氧間互相碰撞后不會重新組合成氧分子。

圖6 原子氧與薄膜材料作用示意圖Fig.6 Schematic diagram of atomic oxygen interaction with membrane

(4)原子氧入射位置在缺陷上方部位隨機生成，空間中由于飛行器的運動原子氧速度在Maxwell分布的基礎(chǔ)上附加了飛行器的速度，其運動方向也由于附加了飛行器速度不再均勻而呈角分布，正面撞擊時主要集中在材料的負法線附近。

(5)原子氧與材料的熱同化過程在仿真中以熱同化概率進行等效。單個原子氧產(chǎn)生的效應(yīng)可以看成是原子氧逐個作用效應(yīng)的累積?？臻g中原子氧濃度很低，原子氧的平均自由程很大。原子氧在掏蝕空洞中反射所經(jīng)歷的路程與平均自由程相比很小，在此過程中與其他原子氧碰撞的機會很少，所以原子氧相互碰撞對基蝕空洞發(fā)展的影響可以忽略。

(6)反應(yīng)的原子氧以兩種方式離開材料表面:

鏡面反射:原子氧發(fā)生彈性碰撞，此時原子氧離開表面時能量保持不變，因而下次碰撞時的反應(yīng)概率不變。

漫反射:原子氧發(fā)生非彈性碰撞，漫反射的原子氧以Maxwell速度分布離開材料表面，散射角度服從余弦分布。

4.2 數(shù)值模擬算法

本文使用MATLAB平臺開發(fā)了蒙特卡洛原子氧掏蝕仿真模擬模型。入射原子氧撞擊到鍍層上或通過缺陷打到聚合物基體上，根據(jù)反應(yīng)概率判斷是否發(fā)生反應(yīng)，若發(fā)生反應(yīng)則將格點從陣列中去掉；若不發(fā)生反應(yīng)則原子氧從撞擊點發(fā)生鏡面反射或漫反射。反射后的原子氧繼續(xù)運動，如果碰到新的格點便重復(fù)上述過程。當(dāng)原子氧發(fā)生反應(yīng)或從缺陷處逸出到空間，上一個原子氧作用的模擬便完成了，可以開始下一個原子氧的模擬。仿真算法流程如圖7所示。采用多維數(shù)組存儲仿真中的格點。

圖7 原子氧掏蝕仿真流程圖Fig.7 AO undercutting simulation algorithm

4.3 關(guān)鍵參數(shù)

如表1所示為仿真用到的反應(yīng)參數(shù)。

表1 蒙特卡洛仿真參數(shù)[5]Table 1 Parameters for AO interaction with deorbit sail

(1)原子氧個數(shù)。蒙特卡洛仿真法采用有限多個仿真分子代替大量的真實分子，合理的仿真分子數(shù)的選取要兼顧計算效率及統(tǒng)計得到的宏觀物理量的真實性。本文在數(shù)值模擬過程中，首先根據(jù)LDEF空間飛行試驗確定大缺陷近似材料無保護層時的模擬分子數(shù)，然后參照該數(shù)據(jù)對數(shù)值模擬的結(jié)果進行修正，進而確定小缺陷時模擬分子數(shù)[15]。如式 (1)、式 (2)所示:

式中，M為大缺陷仿真分子數(shù)，atoms；F為模擬原子氧累積通量，atoms/cm2；L為缺陷長度，cm/cell；Wc為大缺陷寬度，cm；λ為原子氧對離軌帆材料侵蝕系數(shù)，cm3/atom；H為大缺陷網(wǎng)格數(shù)，cell；W為小缺陷網(wǎng)格數(shù)，cell。

(2)反應(yīng)概率。原子氧與聚合物材料的初始反應(yīng)概率如式 (3)[14]所示:

式中，k為系數(shù)，k=0.1165；E0為反應(yīng)活化能，eV；E為原子氧撞擊能量，eV。

后續(xù)反應(yīng)概率取決于原子氧能量和入射角度，如式 (4)[14]所示:

式中，A為系數(shù)，A=0.3596；Pini為初始反應(yīng)概率；E為原子氧撞擊能量，eV；θ為碰撞時入射原子氧與材料法線間夾角，deg。

5 仿真校驗與結(jié)果討論

將NASA長期暴露試驗的數(shù)據(jù)與數(shù)值模擬結(jié)果進行對比，從而檢驗?zāi)Ｐ偷挠行?。模擬的基本數(shù)據(jù)是通過大量飛行試驗的結(jié)果綜合所得的。圖8所示為根據(jù)美國長期暴露試驗 (Long Duration Exposure Facility，LDEF)，掃描電鏡測試結(jié)果畫出的掏蝕空腔輪廓。選取Kapton材料為例，原子氧通量為5.77×1021atoms/cm2，缺陷寬度為2μm進行模擬，從兩幅圖的對比可以看出，無論是從空腔的輪廓形貌還是從掏蝕的深度、寬度和角度考慮，本文的仿真方法都較好地模擬了空間飛行時原子氧對 Kapton材料的作用過程。

圖8 LDEF試驗原子氧掏蝕形貌對比圖Fig.8 AO undercutting profile for LDEF and Monte Carlo simulation results

另外，缺陷底部呈現(xiàn)出兩邊深，中間相對淺的倒三角鋸齒狀形貌。原子氧碰撞速度由其初始速度與飛船運動速度合成，初次入射的原子氧主要集中在材料界面負法線兩側(cè)呈角分布，又很難集中于底部中心區(qū)域。如圖9所示為原子氧初始碰撞角度統(tǒng)計圖。

圖9 原子氧初始碰撞角度統(tǒng)計圖Fig.9 AO incidence angle probability distribution

由于離軌帆帆面折疊收攏以及發(fā)射過程中均不可避免地會產(chǎn)生較為密集的缺陷，缺陷密度與薄膜膜厚具有可比性，因此需對帆面薄膜材料的原子氧掏蝕缺陷密度效應(yīng)進行研究。如圖10(a)、(b)所示分別為原子氧通量為4.67×1020atoms/cm2下500nm和1μm缺陷的掏蝕形貌。掏蝕深度分別為1.2μm和1.4μm，寬度分別為600nm和1.1μm。圖10(c)所示為1μm和500nm的兩個缺陷，缺陷間隔1μm，可以看出在原子氧通量為4.67×1020atoms/cm2的情況下，缺陷間交互效應(yīng)并不明顯。

圖10 原子氧掏蝕情況仿真結(jié)果圖Fig.10 Atomic oxygen erosion simulation results

為了進一步研究原子氧通量對缺陷交互效應(yīng)的影響，如圖11所示分別為原子氧累積通量為4.67 ×1020atoms/cm2、 9.34 ×1020atoms/cm2、2.95×1021atoms/cm2的情況下原子氧對薄膜帆掏蝕形貌的仿真圖。由仿真分析可知，隨著原子氧通量的增加，缺陷交互效應(yīng)越來越明顯。當(dāng)原子氧累積通量達到9.34×1020atoms/cm2時，雖然通量只增加了兩倍，但掏蝕深度增加了近10倍。當(dāng)原子氧累積通量達到2.95×1021atoms/cm2的情況下，離軌帆薄膜材料已經(jīng)被完全穿透，帆面會造成嚴(yán)重損傷，影響其功能的實現(xiàn)。

6 結(jié)論

(1)本文所提出的蒙特卡洛數(shù)值仿真法較為有效，可以作為離軌帆原子氧剝蝕情況預(yù)測，仿真結(jié)果與LDEF空間試驗結(jié)果有較好的一致性。

圖11 不同原子氧通量下薄膜掏蝕情況仿真結(jié)果圖Fig.11 Atomic oxygen undercutting simulation under different fluence

(2)離軌帆在折疊收攏過程中不可避免地會產(chǎn)生折痕缺陷，凹面鍍鋁層在壓應(yīng)力作用下形成條狀剝蝕，條間距為2～3μm量級，條狀裂紋寬度約為500nm～1μm量級。凸面鍍鋁層在拉應(yīng)力作用下形成鱗片狀剝蝕，典型鱗片尺度為2μm×2μm～5μm ×5μm 量級， 鱗狀裂紋寬度約為500nm～1μm量級。這些缺陷部位將成為原子氧集中作用的脆弱部位。

(3)通過軌道分析工具包計算出金牛座離軌帆的離軌時間最短約為37個月，在太陽活動強度和地磁強度取平均值，離軌時間取最小值的情況下，經(jīng)歷的原子氧累積通量約為4.67×1020atoms/cm2。

(4)離軌帆鍍鋁層折疊后會導(dǎo)致裂紋缺陷，在缺陷部位原子氧掏蝕效應(yīng)嚴(yán)重。通過仿真可知，即使是有缺陷的鍍鋁膜也比無鍍鋁層的純薄膜具有更好的原子氧耐受性。在理想的情況下，500nm和1μm的裂紋缺陷，在離軌過程中掏蝕深度和寬度分別為約1.2μm和1.4μm，寬度約為600nm和1.1μm，帆面可以保持基本完整。