厲夫兵，冷俊敏，李紅蓮

(北京信息科技大學(xué) 信息與通信工程學(xué)院，北京 100101)

誘餌球鏡面反射內(nèi)壁的輻射傳遞系數(shù)計(jì)算*

厲夫兵，冷俊敏，李紅蓮

(北京信息科技大學(xué) 信息與通信工程學(xué)院，北京 100101)

針對(duì)內(nèi)表面具備鏡面反射特性的空間誘餌球，基于蒙特卡羅方法建立漫發(fā)射模型，推導(dǎo)光線的最大反射次數(shù)，利用光線在球腔的傳播特性求解光線與球壁的多次反射交點(diǎn)，計(jì)算球面微元吸收的光束能量，獲取不同反射率情況下單一球面微元對(duì)全部球面微元的輻射傳遞系數(shù)。計(jì)算結(jié)果表明，單一微元對(duì)球面各微元的輻射傳遞系數(shù)中存在2個(gè)峰值，分別為微元對(duì)其自身的輻射傳遞系數(shù)，以及微元對(duì)其球?qū)ΨQ微元的輻射傳遞系數(shù)。隨著反射率的增大，微元對(duì)其自身的輻射傳遞系數(shù)逐漸接近于微元對(duì)其球?qū)ΨQ微元的輻射傳遞系數(shù)；當(dāng)反射率接近于1時(shí)，兩者的值近似相等。

輻射傳遞系數(shù)；輻射換熱；鏡面反射；蒙特卡羅；角系數(shù)；球

0 引言

在彈道目標(biāo)飛行中段釋放氣球誘餌的方法可以降低目標(biāo)的探測(cè)概率[1-2]，增加天基紅外監(jiān)視系統(tǒng)對(duì)目標(biāo)的識(shí)別與跟蹤難度[3-4]。公開(kāi)文獻(xiàn)表明，空間氣球多采用鋁膜-聚酯薄膜-鋁膜的超薄層結(jié)構(gòu)[5]。在日光照射下，氣球向陽(yáng)面與背陰面的溫差會(huì)加劇氣球內(nèi)部的輻射換熱，而輻射換熱又反過(guò)來(lái)影響著目標(biāo)的溫度分布。因此，提高輻射換熱的計(jì)算精度對(duì)提高目標(biāo)溫度的計(jì)算準(zhǔn)確度具有重要意義。

近年來(lái)，國(guó)內(nèi)外許多學(xué)者對(duì)影響彈頭、誘餌等空間目標(biāo)的溫度與紅外輻射特性的主要因素進(jìn)行了研究[6-8]，獲得了目標(biāo)在日照區(qū)及地球陰影區(qū)的整體溫度或溫度分布隨時(shí)間的變化[9-10]。傳統(tǒng)方法計(jì)算目標(biāo)內(nèi)部輻射換熱時(shí)，通常將目標(biāo)內(nèi)壁視為灰體面元[11-12]，利用輻射角系數(shù)建立傳熱方程，通過(guò)多次迭代求解灰體面元間的輻射換熱，整個(gè)求解過(guò)程較費(fèi)時(shí)間。張偉清等先利用輻射角系數(shù)求解灰體面元間的輻射傳遞系數(shù)[13]，再直接利用輻射傳遞系數(shù)求解灰體面元間的輻射換熱，避免了直接利用角系數(shù)求解輻射換熱時(shí)的迭代計(jì)算過(guò)程，加速了計(jì)算進(jìn)程。值得注意的是，當(dāng)氣球內(nèi)壁鋁膜足夠光滑時(shí)，氣球內(nèi)表面將具備良好的鏡面反射特性而非漫反射特性，此時(shí)利用基于面元灰體假設(shè)的計(jì)算方法將會(huì)給輻射換熱計(jì)算引入新的誤差，因此需要計(jì)算鏡反射面元間的輻射傳遞系數(shù)。

由于鏡面反射的存在，面元的輻射熱量不可避免地在球腔內(nèi)多次反射，因此優(yōu)先采用蒙特卡羅方法計(jì)算面元間的輻射傳遞系數(shù)。經(jīng)典的蒙特卡羅方法利用面元的發(fā)射率隨機(jī)吸射光線，通過(guò)統(tǒng)計(jì)面元吸收光線數(shù)目的方法計(jì)算輻射傳遞系數(shù)[14-15]，這種隨機(jī)吸收光線的方法在面元發(fā)射率ε較大時(shí)會(huì)導(dǎo)致反射光線數(shù)量的急劇減少。例如，當(dāng)系統(tǒng)發(fā)射率ε=0.9、面元發(fā)射光線數(shù)量N=10 000條時(shí)，經(jīng)過(guò)系統(tǒng)的1，2，3次反射后，剩余光線數(shù)量分別約為1 000，100，10。多次反射發(fā)生后，由于剩余光線的樣本數(shù)量太少，難以對(duì)某面元的反射光線在半空間的分布情況進(jìn)行分析。通過(guò)大量增加發(fā)射光線數(shù)量的方法雖然能夠緩解上述問(wèn)題，但同時(shí)也引入了大量的碰撞與求交計(jì)算，限制了蒙特卡羅方法的計(jì)算速度。

本文通過(guò)建立鏡反射內(nèi)壁面元的蒙特卡羅漫發(fā)射模型，根據(jù)光線的最大反射次數(shù)計(jì)算每次反射事件發(fā)生時(shí)光線損失的能量份額，利用光線在球腔內(nèi)的傳播特性直接求解光線與內(nèi)壁的多次反射交點(diǎn)，通過(guò)統(tǒng)計(jì)各次反射事件中面元吸收的光線能量份額計(jì)算面元間的輻射傳遞系數(shù)，在此基礎(chǔ)上分析鏡反射內(nèi)壁的輻射傳遞系數(shù)分布特性。

1 蒙特卡羅漫發(fā)射模型

1.1 目標(biāo)建模與微元剖分

氣球目標(biāo)被建模為半徑為R0的球面，并在天頂角和方位角方向進(jìn)行等間隔剖分，獲得一系列球面微元，并利用微元在天頂角、方位角方向的位序確定微元的位置。由球的性質(zhì)可知，內(nèi)壁微元表面任意位置點(diǎn)的法矢量為該點(diǎn)指向球心的單位矢量。

1.2 光線的隨機(jī)發(fā)射位置、方向與能量

在蒙特卡羅方法中，面元向半空間的輻射能量被視為從面元隨機(jī)位置、以隨機(jī)方向發(fā)射的光線的集合，每束光線都具備一定的能量[14-15]。假定面元向半空間發(fā)射的光線總數(shù)為N，則每束光線攜帶的初始能量為

(1)

式中：T為面元的瞬態(tài)溫度；ρ為面元的反射率；σ為Steven-Boltzmann常數(shù)。

在球坐標(biāo)系中，光線在球面微元內(nèi)表面的隨機(jī)發(fā)射位置可表示為

(2)

式中：θa，θb分別是球面微元天頂角的上、下界；φa，φb分別是微元方位角的上、下界；Rθ，Rφ是在[0，1]區(qū)間滿足均勻分布的隨機(jī)變量。

光線在隨機(jī)發(fā)射位置點(diǎn)上的隨機(jī)發(fā)射方向的求解過(guò)程可參考文獻(xiàn)[14]。通過(guò)分析面元向半空間范圍內(nèi)投射光束的二維概率密度，可推導(dǎo)光束的二維分布函數(shù)并得到光束的發(fā)射方向與隨機(jī)數(shù)的關(guān)系為

(3)

式中：θf(wàn)，φf(shuō)分別為發(fā)射光束的天頂角、方位角；Rθ，Rφ為[0，1]區(qū)間內(nèi)均勻分布的隨機(jī)數(shù)。

2 輻射傳遞系數(shù)的計(jì)算

2.1 光線在球腔中的傳播

圖1所示為參考坐標(biāo)系OXYZ中球面微元A上光束的某次隨機(jī)發(fā)射位置p(θp，φp)。為計(jì)算方便，建立與參考坐標(biāo)系OXYZ重合的本體坐標(biāo)系Oxyz，并將本體坐標(biāo)系Oxyz先沿z軸逆時(shí)針旋轉(zhuǎn)角度φp，再沿y軸順時(shí)針旋轉(zhuǎn)角度π-θp，使該坐標(biāo)系的z軸負(fù)半軸恰好通過(guò)光線的發(fā)射位置p。因此，從本體坐標(biāo)系Oxyz向參考坐標(biāo)系OXYZ的坐標(biāo)轉(zhuǎn)換矩陣為

(4)

式中：Ly(β)，Lz(γ)分別為坐標(biāo)系沿y軸、z軸逆時(shí)針旋轉(zhuǎn)β，γ角度對(duì)應(yīng)的旋轉(zhuǎn)矩陣，即

圖1 面元隨機(jī)發(fā)射位置p在坐標(biāo)系OXYZ與Oxyz中的位置Fig.1 Random position ‘p’ on a spherical facet A in OXYZ coordinates and Oxyz coordinates

如圖2所示，在坐標(biāo)系Oxyz中，光線從點(diǎn) p[0，0，-R0] 以角度 (θf(wàn)，φf(shuō)) 發(fā)射。光線與球壁第1次相交于點(diǎn)d (1次反射點(diǎn))，并遵循鏡面反射定律繼續(xù)向前傳播。由于球內(nèi)表面上任意點(diǎn)的法矢量都指向球心，故過(guò)點(diǎn)d的反射光線必定在由球心、點(diǎn)p和點(diǎn)d所確定的大圓切平面內(nèi)。容易推得，發(fā)射光線的多次反射光線也全部落在該圓內(nèi)。從圖2中可以看出，發(fā)射光線與球壁的1次反射點(diǎn)d在大圓中的圓心角α1=π-2θf(wàn)，光線與球壁的k次反射點(diǎn)對(duì)應(yīng)的圓心角為

(5)

式中：%表示求余操作。

容易求出，k次反射點(diǎn)在坐標(biāo)系Oxyz中的天頂角θk和方位角φk為

(6)

將坐標(biāo)系Oxyz中的k次反射點(diǎn)坐標(biāo) (R0，θk，φk) 先轉(zhuǎn)換為直角坐標(biāo)，再利用式(4)的坐標(biāo)轉(zhuǎn)換矩陣L求反射點(diǎn)在坐標(biāo)系OXYZ中的直角坐標(biāo)，并將該直角坐標(biāo)轉(zhuǎn)換為球坐標(biāo)(R0，θr，φr)，最后利用反射點(diǎn)在坐標(biāo)系OXYZ中的天頂角θr和方位角φr求出反射點(diǎn)所在的球面微元。

圖2 光線在球腔中的多次反射 (坐標(biāo)系Oxyz)Fig.2 Multi-reflections of a light ray inside the sphere in Oxyz coordinates

2.2 輻射傳遞系數(shù)的計(jì)算

輻射傳遞系數(shù)是離開(kāi)一個(gè)表面到達(dá)另一個(gè)表面的能量份額，其中包含了從任意面元反射后到達(dá)目標(biāo)表面的能量。因此，球面微元的輻射傳遞系數(shù)需要考慮光線的多次反射貢獻(xiàn)。

本文中，單束光線攜帶的初始能量e0是被面元依吸收率逐漸吸收而衰減的。當(dāng)光線發(fā)生第1次反射時(shí)，反射光束的能量變?yōu)棣裡0，被交點(diǎn)所屬微元吸收的能量為(1-ρ)e0。當(dāng)光線發(fā)生第2次反射時(shí)，反射光束的能量變?yōu)棣?e0，被交點(diǎn)所屬微元吸收的能量為(1-ρ)ρe0。以此類推，第k次反射時(shí)，反射光束的能量變?yōu)棣裬e0，被交點(diǎn)所屬微元吸收的能量為

(7)

給定閾值η，當(dāng)?shù)趉次反射發(fā)生時(shí)，若反射光束的能量ρke0與初始能量e0之比低于閾值η時(shí)視為光束能量被完全吸收，即ρk≤η，則最小反射次數(shù)nr的取值可表示為

(8)

式中：「?表示上取整。

綜上所述，對(duì)于球面微元的任意一束初始發(fā)射光線，利用式(5)～(7)可一次求得其nr個(gè)反射點(diǎn)的所屬微元及微元吸收能量。對(duì)微元的全部N束發(fā)射光線進(jìn)行遍歷，統(tǒng)計(jì)每束光線的全部nr個(gè)反射點(diǎn)的所屬微元，并計(jì)算微元吸收能量，可求得該微元對(duì)其它微元的輻射傳遞系數(shù)為

(9)

式中：Ej為微元i的輻射能量中最終被微元j吸收的部分，包括微元i對(duì)微元j的直接入射能量和被其它面元反射后到達(dá)微元j的能量。

3 數(shù)值計(jì)算結(jié)果與分析

本例的參數(shù)設(shè)置如下。氣球半徑R0=1 m，內(nèi)壁鏡面鋁膜的反射率為ρ，球面沿天頂角剖分18份 (dθ=10°)，沿方位角方向剖分36份 (dφ=10°)。對(duì)特定的球面微元，隨機(jī)發(fā)射光線的數(shù)目N=100萬(wàn)條。光束被多次吸收、反射后，當(dāng)反射光束的能量與初始能量之比低于閾值η=0.01時(shí)，輻射能量被視為全部吸收，光線的多次反射計(jì)算過(guò)程結(jié)束。

3.1 輻射角系數(shù)

當(dāng)球內(nèi)壁為黑體時(shí) (ε=1，ρ=0)，輻射傳遞系數(shù)等于輻射角系數(shù)。球內(nèi)壁任意兩球面微元m，n之間的角系數(shù)理論值為[15]

(10)

式中：Sn為微元n的面積；θa，θb分別為微元n的天頂角上、下界；φa，φb分別為微元的方位角上、下界。

由式(10)可知，球面微元間的角系數(shù)僅與微元n的面積有關(guān)，與發(fā)射微元m的位置和面積無(wú)關(guān)。

圖3a)為球面發(fā)射微元A(天頂角、方位角中心為 (135°，115°)) 對(duì)球內(nèi)壁全部微元的角系數(shù)計(jì)算值二維展開(kāi)圖，垂直方向?yàn)榍騼?nèi)壁微元的天頂角，水平方向?yàn)榉轿唤?。圖3b)為3個(gè)特定的方位角下 (對(duì)應(yīng)于圖3a)方框) 的輻射角系數(shù)計(jì)算值與理論值對(duì)比，圖3c)為計(jì)算值與理論值的相對(duì)誤差。

圖3 球面微元的角系數(shù)理論值與計(jì)算值Fig.3 Theoretical and calculated shape factors for sphere facets

從圖3a)可以看出，該微元對(duì)球面各微元的輻射角系數(shù)在方位角方向保持相對(duì)穩(wěn)定，僅在天頂角方向變化。從圖3b)～c)可以看出，輻射角系數(shù)計(jì)算值與理論值很接近。在本例的參數(shù)設(shè)置下，相對(duì)誤差大部分處于±5%的范圍內(nèi)，最大相對(duì)誤差不超過(guò)10%。由理論式(10)可知，當(dāng)球面沿方位角方向均勻剖分時(shí)，φb，φa為常數(shù)，角系數(shù)與方位角無(wú)關(guān)，由此證明了本文方法計(jì)算輻射角系數(shù)的正確性。

3.2 輻射傳遞系數(shù)

圖4a)～c)分別為鏡面反射率ρ=0.40，0.70，0.99時(shí) (反射次數(shù)nr=6，13，459)，發(fā)射微元A(135°，115°) 對(duì)全體微元的輻射傳遞系數(shù)矩陣。從圖中可以看出，與圖3a)的黑體面元不同，鏡面反射時(shí)，輻射傳遞系數(shù)沿方位角方向也發(fā)生變化，并呈現(xiàn)出幾個(gè)特性。①隨著反射率的增大，輻射傳遞系數(shù)存在兩個(gè)較為明顯的峰值區(qū)域，一個(gè)恰為微元A(135°，115°)自身，另一個(gè)是其球心對(duì)稱微元B(45°，295°) (注：球面兩點(diǎn)a，b關(guān)于球心對(duì)稱時(shí)，其天頂角滿足θa+θb=180°，方位角滿足 |φa-φb|=180°)。②當(dāng)反射率較低時(shí)，微元A的輻射傳遞系數(shù)小于微元B，隨著反射率的增加，其輻射傳遞系數(shù)逐漸接近微元B。當(dāng)反射率接近于1時(shí)，兩者的值近似相等。

圖4 發(fā)射微元A對(duì)球面微元的輻射傳遞系數(shù)Fig.4 Radiative heat transfer coefficient between a manua-lly selected facet A and the whole spherical facet

圖5a)～b)為反射率ρ=0.99時(shí) (nr=459)，2個(gè)不同位置的發(fā)射微元A(5°，55°)，(85°，85°) 對(duì)全體球面微元的輻射傳遞系數(shù)矩陣。從圖5中可以看出，與圖4c)類似，輻射傳遞系數(shù)矩陣中的兩個(gè)峰值所在的微元在空間上關(guān)于球心對(duì)稱。由此可以推知，當(dāng)反射率較大時(shí)，任意位置的發(fā)射微元對(duì)全體球面微元的輻射傳遞系數(shù)均存在2個(gè)峰值，且分別位于微元自身及其球心對(duì)稱微元。

接下來(lái)首先分析圖5中輻射傳遞系數(shù)存在2個(gè)關(guān)于球心對(duì)稱的峰值的原因。在圖2所示的坐標(biāo)系Oxyz中，假設(shè)微元A，B分別為發(fā)射微元及其球心對(duì)稱微元。當(dāng)反射率ρ增大時(shí)，單束光線的反射次數(shù)將增加；當(dāng)ρ接近于1時(shí)，反射次數(shù)nr將非常大，各次反射點(diǎn)將散布于大圓圓周上。雖然單一光線的多次反射點(diǎn)在大圓圓周上的散布位置是確定的；但對(duì)于發(fā)射天頂角θf(wàn)∈[0，π/2]的大量隨機(jī)發(fā)射光線而言，其多次反射點(diǎn)落在大圓圓周各處的概率將近似相等。由于大圓切割球心對(duì)稱微元A，B形成的弧長(zhǎng)相等，因此，在方位角φ=φf(shuō)上，反射點(diǎn)落在微元A上弧線的概率等于其落在微元B上弧線的概率。特別地，由圖2可知，由于任意方位角φ∈[0，2π)的大圓都經(jīng)過(guò)微元A，B，因此反射點(diǎn)落在微元A，B上的概率遠(yuǎn)高于其他微元。即反射率接近于1時(shí)，微元A，B的輻射傳遞系數(shù)近似相等且遠(yuǎn)高于其他球面微元。

圖5 ρ=0.99時(shí)不同位置的發(fā)射微元A對(duì)應(yīng)的輻射傳遞系數(shù)矩陣Fig.5 Radiative heat transfer coefficient matrix corresponding to two different manually selected facets A when reflectivity ρ is 0.99

對(duì)圖4a)～c)中微元A的輻射傳遞系數(shù)隨反射率ρ的增大而逐漸逼近微元B的原因進(jìn)行分析時(shí)，需要對(duì)微元A，B吸收的光線能量成份進(jìn)行研究。通過(guò)對(duì)微元A的每條發(fā)射光線進(jìn)行追蹤，記錄微元上的各個(gè)反射點(diǎn)都源于光線的第幾次反射行為，然后統(tǒng)計(jì)微元上各次反射行為發(fā)生的次數(shù)。圖6所示為在圖4參數(shù)設(shè)置下，反射率ρ=0.40，0.70，0.83時(shí)，微元A，B上光線的第k次反射行為 (k=1，2，…，nr) 發(fā)生的次數(shù)統(tǒng)計(jì)。

圖6 不同反射率時(shí)面元A，B上各次反射發(fā)生的次數(shù)Fig.6 Statistics of multi-reflection of rays on facet A and its spherically symmetric facet B

當(dāng)η=0.01，ρ=0.4時(shí)，由式(8)可知光線發(fā)生6次反射即被完全吸收。如圖6a)所示，當(dāng)k=1時(shí)微元A，B上光線的第1次反射行為發(fā)生的次數(shù)近似相等，即微元A發(fā)出的光線中，與微元A，B直接發(fā)生碰撞的光線數(shù)量近似相等。由于微元A，B是關(guān)于球心對(duì)稱的微元，面積相等，因此光線與兩面元發(fā)生碰撞的概率相等，這與節(jié)3.1中輻射角系數(shù)的計(jì)算結(jié)論是一致的。

k=2表示光線的第2次反射行為，即發(fā)射微元A發(fā)出的光線經(jīng)過(guò)球面微元的1次鏡面反射后再與微元A或B相交。由圖2可知，光線由微元A發(fā)出后，只有那些能夠與微元B(或其附近微元) 發(fā)生直接碰撞的光線才有可能經(jīng)1次反射后與微元A相交，此部分光線數(shù)量較少。但是，光線由微元A發(fā)出后，所有能夠與“赤道”微元 (與Oxy平面相交的球面微元) 發(fā)生直接碰撞的光線經(jīng)過(guò)1次鏡面反射后，都有可能與微元B相交。因此，2次反射點(diǎn)落在微元B上的光線數(shù)量要遠(yuǎn)大于落在微元A上的光線數(shù)量，即微元B吸收的光線的2次碰撞能量遠(yuǎn)大于微元A，直接導(dǎo)致此時(shí)微元B的輻射傳遞系數(shù)較微元A大。

由圖6a)～c)可以看出，當(dāng)反射率ρ增大時(shí)，反射次數(shù)增多，且光線的第k次反射行為 (k≥3) 發(fā)生在微元A，B上的數(shù)目趨于相等。此時(shí)，雖然微元B吸收的光線的2次碰撞能量仍較微元A高，但相對(duì)于總能量而言，其所占的能量份額隨反射率的增大而逐漸減小。因此，如圖4所示，隨著反射率的增大，微元A的輻射傳遞系數(shù)逐漸接近于微元B。

在本例中，光線被多次反射后，其攜帶的能量會(huì)逐漸降低，若反射光射的能量與初始能量之比低于閾值η=0.01，則光線能量被視為全部吸收，即光束被吸收的能量份額為0.99。當(dāng)該閾值減小時(shí)，計(jì)算精度將會(huì)提高。圖7所示為η=0.000 1、其他設(shè)置與圖5相同時(shí)，得到的輻射傳遞系數(shù)矩陣。

圖7 η=0.000 1，ρ=0.99時(shí)不同位置的微元A對(duì)應(yīng)的輻射傳遞系數(shù)矩陣 Fig.7 Radiative heat transfer coefficient matrix corresponding to two different manually selected facets when η=0.000 1，ρ is 0.99

從圖7與圖5可以看出計(jì)算結(jié)果非常相近。這是因?yàn)楫?dāng)閾值η=0.000 1時(shí)，光束被吸收的能量份額為0.999 9；而當(dāng)η=0.01時(shí)，光束被吸收的能量份額為0.99。由此可見(jiàn)，兩者的能量變化不足1%，因此圖7得到的結(jié)果與圖5類似，但數(shù)據(jù)精度更高。

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

本文利用蒙特卡羅方法計(jì)算了球內(nèi)壁的輻射角系數(shù)及不同鏡面反射率情況下的輻射傳遞系數(shù)。輻射角系數(shù)的計(jì)算值與理論值一致。輻射傳遞系數(shù)的結(jié)果表明，首先，微元對(duì)全部球面微元的輻射傳遞系數(shù)矩陣中存在2個(gè)峰值，即微元對(duì)其自身、以及微元對(duì)其球心對(duì)稱微元的輻射傳遞系數(shù)。其次，隨著反射率的增大，微元對(duì)其自身的輻射傳遞系數(shù)逐漸接近于微元對(duì)其球心對(duì)稱微元的輻射傳遞系數(shù)。當(dāng)反射率接近于1時(shí)，兩者的輻射傳遞系數(shù)值近似相等。由此可知，對(duì)于誘餌球鏡反射內(nèi)壁微元而言，當(dāng)某微元的輻射出射能量一定時(shí)，反射率越大，輻射能量越傾向于向該微元及其球?qū)ΨQ微元集中。本文的結(jié)論有助于分析空間誘餌球向陽(yáng)面高溫微元與背陰面低溫微元之間的輻射換熱與溫度變化過(guò)程。

[1] 酈蘇丹，任萱.大氣層外誘餌釋放研究[J].宇航學(xué)報(bào)，2001，22(2)：100-104. LI Su-dan，REN Xuan.Research of Releasing Decoy Outside Atmosphere[J].Journal of Astronautics，2001，22(2)：100-104.

[2] 林兩魁，謝愷，徐暉，等.中段彈道目標(biāo)群的紅外成像仿真研究[J].紅外與毫米波學(xué)報(bào)，2009，28(3)：218-223. LIN Liang-kui，XIE Kai，XU Hui，et al.Research on Infrared Imaging Simulation of Midcourse Ballistic Object Target Complex[J].Journal of Infrared and Millimeter Waves，2009，28(3)：218-223.

[3] 王碩，張奕群，孫冰巖.紅外點(diǎn)目標(biāo)跟蹤方法綜述[J].現(xiàn)代防御技術(shù)，2016，44(2)：124-134. WANG Shuo，ZHANG Yi-qun，SUN Bing-yan.Review of Point Target Tracking Methods Based on Infrared Sensor[J].Modern Defence Technology，2016，44(2)：124-134.

[4] 浦甲倫，崔乃剛，郭繼峰.天基紅外預(yù)警衛(wèi)星系統(tǒng)及其探測(cè)能力分析[J].現(xiàn)代防御技術(shù)，2008，36(4)：68-72. PU Jia-lun，CUI Nai-gang，GUO Ji-feng.Space-Based Infrared System and the Analysis of Its Detecting Capability[J].Modern Defence Technology，2008，36(4)：68-72.

[5] SESSLER A M，CORNWALL J M，DIETZ B，et al.Countermeasures：a Technical Evaluation of the Operational Effectiveness of the Planned US National Missile Defense System[R].Cambridge，MA：Union of Concerned Scientists and MIT Security Studies Program，April 2000.

[6] 姚連興，侯秋萍，羅繼強(qiáng).彈道導(dǎo)彈中段目標(biāo)表面溫度與紅外突防研究[J].航天電子對(duì)抗，2005，21(2)：5-16. YAO Lian-xing，HOU Qiu-ping，LUO Ji-qiang.Research on Surface Temperature and Infrared Signature of the Ballistic Warhead in Midcourse[J].Aerospace Electronic Warefare，2005，21(2)：5-16.

[7] 汪洪源，陳赟.天基空間目標(biāo)紅外動(dòng)態(tài)輻射特性建模與仿真[J].紅外與激光工程，2016，45(5)：0504002. WANG Hong-yuan，CHEN Yun.Modeling and Simulation of Infrared Dynamic Characteristics of Space-Based Targets[J].Infrared and Laser Engineering，2016，45(5)：0504002.

[8] 敬韓博，王英瑞.臨近空間高超聲速目標(biāo)天基紅外探測(cè)技術(shù)研究[J].現(xiàn)代防御技術(shù)，2016，44(6)：80-84. JING Han-bo，WANG Ying-rui.Space-Based Infrared Detection for Near Space Hyptrsonic Targets[J].Modern Defence Technology，2016，44(6)：80-84.

[9] ZHU D Q，SHEN W T，CAI G B，et al.Numerical Simulation and Experimental Study of Factors Influencing the Optical Characteristics of a Spatial Target[J].Applied Thermal Engineering，2013(50)：749-762.

[10] 張駿，楊華，凌永順，等.彈道導(dǎo)彈中段彈頭表面溫度場(chǎng)分布理論分析[J].紅外與激光工程，2005，34(5)：583-586. ZHANG Jun，YANG Hua，LING Yong-shun，et al.Theoretical Analysis of Temperature Field on the Surface of Ballistic Missile Warhead in Midcourse[J].Infrared and Laser Engineering，2005，34(5)：583-586.

[11] LI F B，XU X J.Modeling Time-Evolving Infrared Characteristics for Space Object with Mircromotions[J].IEEE Transactions on Aerospace and Electronic Systems，2012，48(4)：3567-3577.

[12] 厲夫兵，許小劍.計(jì)算超薄多層介質(zhì)誘餌溫度的快速算法[J].紅外與激光工程，2011，40(6)：986-991.

LI Fu-bing，XU Xiao-jian.Fast Algorithm for Temperature Calculation of Multi-Layered Ultra-Thin Decoy[J].Infrared and Laser Engineering，2011，40(6)：986-991.

[13] 張偉清，宣益民，韓玉閣.單元表面間輻射傳遞系數(shù)的新型計(jì)算方法[J].宇航學(xué)報(bào)，2005，26(1)：77-81. ZHANG Wei-qing，XUAN Yi-min，HAN Yu-ge.A New Method of Radiative Transfer Coefficient between Unit Surfaces[J].Journal of Astronautics，2005，26(1)：77-81.

[14] HOWELL J R，PERLMUTTER M.Monte Carlo Solution of Thermal Transfer Through Radiant Media between Gray Walls[J].Journal of Heat Transfer，1964，86(1)：116-122.

[15] 卞伯繪.輻射換熱的分析與計(jì)算[M].北京：清華大學(xué)出版社，1988：81-83，144-145. BIAN Be-hui.Analysis and Calculation of Radiative Heat Transfer[M].Beijing:Tsinghua Universtiy Press,1988:81-83,144-145.

Radiative Heat Transfer Coefficient Calculation of Inflatable Sphere with Specular Inner Surface

LI Fu-bing，LENG Jun-min，LI Hong-lian

(Beijing Information Science & Technology University，School of Information & Communication Engineering，Beijing 100101，China)

Monte Carlo technique is used to calculate the radiative heat transfer coefficient for an inflatable sphere with specular reflection inner surface. Light rays are emitted diffusely from a spherical facet and being traced. The maximum number of times a ray can be reflected is derived, and the multi-intersections of a ray with the inner side of spherical surface are obtained. Radiative heat transfer coefficients between a manually selected facet and the whole spherical facet are calculated at different reflectivities. The results indicate that there are two peaks in the calculated coefficients, the smaller one is the coefficient between the selected facet and itself (i.e., the self-to-self coefficient), and the larger one is between the selected facet and its spherically symmetric facet (i.e., the self-to-spherical-symmetry coefficient). The self-to-self radiative heat transfer coefficient gradually approaches and finally equals the self-to-spherical-symmetry coefficient as the reflectivity increases to 1.

radiation transfer coefficient；radiation heat transfer；specular reflection；Monte Carlo；shape factor；sphere

2016-12-01；

2017-02-11

北京市自然科學(xué)基金資助項(xiàng)目(3164043)

厲夫兵(1982-)，男，山東日照人。講師，博士，主要從事目標(biāo)紅外輻射特性建模研究。

10.3969/j.issn.1009-086x.2017.03.030

TK124；TP391.9

A

1009-086X(2017)-03-0193-07

通信地址：100101 北京市朝陽(yáng)區(qū)北四環(huán)中路35號(hào)北京信息科技大學(xué)信息與通信工程學(xué)院

E-mail：lifubing@bistu.edu.cn