張磊
(北京跟蹤與通信技術研究所,北京 100094)
煙幕遮蔽是眾多隱身技術中的一種,已有數十年的軍事應用歷史。目前,紅外煙幕廣泛應用于大氣層內對紅外制導及紅外偵察設備的干擾,對于太空飛行器,煙幕遮蔽也是有效的隱身手段之一。太空飛行器的紅外輻射主要來自溫度引起的自發熱輻射,在空間中釋放煙幕,可以在其周圍形成一定濃度的粒子層,并且與飛行器伴飛較長時間。在此期間,飛行器的紅外輻射需要穿過粒子層才能被探測,由于粒子層具有對紅外光的吸收和散射能力,飛行器可被探測到的能量將減少,以達到隱身的效果[1-5]。
由于太空飛行器載荷有限,煙幕的拋撒方式、粒子的種類和濃度、粒子在空間中的擴散情況都對其紅外隱身效果有巨大影響。本文通過建立光線穿過粒子層的輻射傳輸模型,分析不同擴散情況、不同粒子,不同濃度下煙幕在不同波段的遮蔽效果,為煙幕遮蔽技術的發展提供基礎支撐。
球形粒子的消光因子Qe,散射因子Qs,吸收因子Qa根據MIE散射理論得到[6-9]:
(1)
(2)
Qa(m,χ)=Qe-Qa,
(3)
式中:m為粒子復折射率,m=n-ik,n為折射指數,k為吸收指數;χ為尺度參數;an,bn為MIE散射系數。
散射相函數Φ采用線性各向異性相函數
Φ(Θ)=1+A1(cosΘ),
(4)
式中:A1=1。
煙幕是具有一定厚度、一定濃度的粒子系,這里采用獨立散射假設計算粒子系輻射參數。
(5)
(6)
(7)
式中:β為消光系數;σs為散射系數;k為吸收系數;Ni為粒徑Di的粒子數密度,單位為m-3;Ce,i,Cs,i,Ca,i分別表示粒徑為Di的單粒子的消光截面、散射截面和吸收截面。
粒子屬于吸收、發射、散射性介質,其輻射傳輸方程為
(8)
式中:Iη為光譜輻射強度;s為傳輸距離;η為波數,Ib,η為光譜黑體輻射強度;Ωi為第i個立體角。方程采用有限體積法進行離散,共軛梯度法求解離散的輻射傳輸方程。
本文對單粒子輻射特性參數進行計算,并與文獻[10]中的結果進行對比,如圖1,2所示。其中粒子尺度χ=1,圖1為散射因子1a)為文獻結果,1b)為本文計算結果;圖2為吸收因子2a)為文獻結果,2b)為本文計算結果。計算結果與文獻中的結果吻合良好。

圖1 粒子散射因子Fig.1 Particle scattering factor

圖2 粒子吸收因子Fig.2 Particle absorption factor
根據典型煙幕拋撒狀態,設定太空飛行器周圍均勻散布著一個煙幕層[11-12],煙幕層呈球狀,厚度為1 m,煙幕層距飛行器40 m。如圖3所示煙幕半球截面示意圖。

圖3 煙幕系統建模示意圖Fig.3 Diagram of smoke modeling
煙幕層中粒子的復折射率為1+1i,數密度為1e12 m-3,粒子的直徑為1 μm,粒子自身的溫度為均勻的200 K[13-15]。在這種條件下,基于上述理論模型,計算了飛行器和煙幕層總的輻射輻射強度,表1給出了煙幕遮蔽前后飛行器輻射強度對比。

表1 煙幕遮蔽前后彈頭輻射強度對比Table 1 Comparison of radiation intensity withand without smoke
可以看出,經過煙霧的遮蔽,飛行器在中波紅外的輻射急劇減小,說明煙幕在中波紅外波段具有較強的紅外隱身能力。
煙幕從太空飛行器拋撒出來后會向空間不斷擴散,為了研究擴散過程對煙幕消光特性的影響,分析了煙幕拋撒過程中,不同時刻的輻射強度,設定煙幕按均勻速度向空間均勻擴散,煙幕的總質量保持不變。
圖4給出了3~5 μm波段彈頭輻射強度隨擴散距離的變化。彈頭及煙幕總體的紅外輻射都隨擴散距離的增加而增大。這是因為隨著擴散的增加,煙幕的密度下降,消光能力下降,而隨著煙幕投影面積的增加其自發輻射能量的比重在增加。

圖4 3~5 μm波段彈頭輻射強度隨擴散距離的變化Fig.4 Variation of radiation intensity of warhead with diffusion distance in 3~5 μm
粒子光學屬性、粒徑、煙幕厚度、煙幕密度等參數是影響煙幕隱身效果的主要因素,下面分析這些參數對目標中波紅外的隱身效果的影響。
復折射率是粒子自身的光學性質,其中復折射率的實部(折射指數n)影響粒子的散射能力,復折射率的虛部(吸收指數k)反應粒子的吸收能力。總體輻射強度是綜合粒子復折射率、粒子尺度、濃度的復雜函數如圖5所示,粒子的折射指數對總體輻射的影響并不大。而吸收指數對總體輻射強度有較大影響,并且影響并不是線性的,根據粗略計算,輻射的最大值出現在k=1附近。

圖5 3~5 μm波段輻射隨復折射率的變化規律Fig.5 Variation of radiation with complex refractive index in 3~5 μm
圖6給出了輻射隨粒徑的變化規律,結果表明,總輻射隨粒徑的變化規律較為復雜,在0.1~10 μm的范圍內,在粒徑等于1 μm左右時具有最佳隱身效果,之后隨著粒徑的增加隱身效果有少許下降,但在粒徑<1 μm時,隱身效果驟減。另外,在中波紅外隱身效果最佳的粒徑尺度是會隨著粒子光學性質(復折射率)而發生變化的。

圖6 3~5 μm波段輻射隨粒徑的變化規律Fig.6 Variation of radiation with particle size in 3~5 μm
圖7給出了輻射隨粒子數密度的變化規律,隨著煙幕中粒子數密度的增大,煙幕吸收散射的衰減作用就越強。但當粒子數密度達到一定數量時,煙幕粒子的隱身能力將不會變化。

圖7 3~5 μm波段輻射隨粒子數密度的變化規律Fig.7 Variation of radiation with particle number density in 3~5 μm
圖8為輻射隨煙幕厚度的變化規律,煙幕的隱身能力,在確定的邊界條件下,對煙幕層的厚度并不敏感。隨著厚度的增加彈頭和煙幕的總體輻射強度略有增大,但變化較小。

圖8 3~5 μm波段輻射隨煙幕厚度的變化規律Fig.8 Variation of radiation with smoke thickness in 3~5 μm
在本文中,基于MIE散射理論求解粒子系輻射傳輸模型,針對太空飛行器被煙幕包裹的場景進行了初步輻射計算。基于這些理論分析,在設定條件下,得到了以下幾點初步結論:
(1) 在特定的波段內,煙幕對太空飛行器能夠起到較大隱身作用。在本文設定的邊界條件中,煙幕在3~5 μm波段紅外隱身效果較好,在最好的情況下,90%的飛行器自發輻射都無法穿過煙幕層。
(2) 煙幕拋撒出來之后,隨著擴散進行,煙幕的隱身效能呈下降趨勢,因此煙幕內維持多久的隱身效能取決于煙幕擴散的速度,也就是煙幕拋撒的技術。
(3) 按照初步計算,粒子的折射指數對煙幕隱身效果影響較小,吸收指數對煙幕隱身效果影響較大。通常粒子的折射指數比較穩定,吸收指數隨相態變化而變化,可以通過改變粒子的相態改變吸收指數,已達到較好的效果。
(4) 粒子的尺度對隱身效果影響較大。煙幕光學性質會影響最佳粒徑的選擇,因此在選擇粒徑和煙幕材料是要結合考慮。
(5) 煙幕的密度存在飽和值,超出飽和值后隱身效果不會增加;煙幕厚度對隱身效果沒有明顯改善;因此考慮飛行器載荷限制,理論上應該存在最佳的煙幕拋灑量。