黃敏超，吳 尚

(國防科技大學(xué) 航天科學(xué)與工程學(xué)院， 湖南 長沙 410073)

?

太陽能熱推進(jìn)系統(tǒng)聚光器參數(shù)優(yōu)化*

黃敏超，吳 尚

(國防科技大學(xué) 航天科學(xué)與工程學(xué)院， 湖南 長沙 410073)

采用先進(jìn)遺傳算法對太陽能熱推進(jìn)系統(tǒng)一次聚光器參數(shù)進(jìn)行優(yōu)化，達(dá)到減小推進(jìn)系統(tǒng)質(zhì)量的目的。以聚光器太陽光收集效率和聚光器質(zhì)量為優(yōu)化目標(biāo)函數(shù)，建立聚光器工作的數(shù)學(xué)模型，并開展相關(guān)仿真研究。仿真結(jié)果表明，先進(jìn)遺傳算法可有效用于太陽光收集效率和聚光器質(zhì)量優(yōu)化分析。

聚光器；收集效率；質(zhì)量；優(yōu)化

太陽能熱推進(jìn)系統(tǒng)一次聚光器是整個系統(tǒng)的重要組成部分，決定著供給能量的大小，影響著整個系統(tǒng)的光熱轉(zhuǎn)換效率和推力器的性能。在太陽能熱推進(jìn)系統(tǒng)一次聚光器的設(shè)計方案中。理想的旋轉(zhuǎn)拋物面聚光器是太陽能熱微推進(jìn)系統(tǒng)一次聚光器的最佳選擇[1]。理想的旋轉(zhuǎn)拋物面聚光器具有良好的聚光特性,能把平行于光軸的入射光線匯聚在焦點(diǎn)上，具有功率高、聚光比高、質(zhì)量輕和體積小的優(yōu)點(diǎn)。旋轉(zhuǎn)拋物面聚光器的聚光原理及光路如圖1所示，聚光器的頂點(diǎn)在旋轉(zhuǎn)拋物面的正中心，平行于光軸入射到拋物線上的任一點(diǎn)的光線SM，其反射光線都匯聚于拋物線的焦點(diǎn)F上。旋轉(zhuǎn)拋物面是拋物線繞光軸旋轉(zhuǎn)一周得到的,因此拋物面上的點(diǎn)都具備同樣的聚光性質(zhì),平行于光軸入射到拋物面上的所有反射光線都將聚集在焦點(diǎn)上,這就是旋轉(zhuǎn)拋物面聚光器聚光的原理[2]。

在滿足微小衛(wèi)星的任務(wù)及動力配置基礎(chǔ)上，要求太陽能熱推力器推力F至少為100 mN，比沖Isp至少為800 s。以一次聚光器安裝的接受半角φ和焦距f為自變量，以光的收集效率和一次聚光器質(zhì)量為目標(biāo)函數(shù)，采用Isight軟件中的領(lǐng)域培植遺傳算法(Neighborhood Cultivation Genetic Algorithm, NCGA)對一次聚光器進(jìn)行優(yōu)化，得到適合于空間應(yīng)用的太陽能熱推進(jìn)系統(tǒng)一次聚光器結(jié)構(gòu)參數(shù)。

圖1 旋轉(zhuǎn)拋物面光路圖

1 聚光器性能評估參數(shù)

對于微小衛(wèi)星來說，其推進(jìn)系統(tǒng)主要用于其軌道保持和姿態(tài)控制，決定了其成本低、質(zhì)量小、體積小、性能高、研制周期短等特點(diǎn)。綜合考慮太陽能熱推進(jìn)聚光器系統(tǒng)性能評估，主要參數(shù)包括：光收集效率和質(zhì)量。太陽能熱微推進(jìn)聚光器的優(yōu)化是由硬件設(shè)備的經(jīng)濟(jì)性和所要求的收集能力共同決定的，即聚光器系統(tǒng)想要獲得較高的光收集效率，就需要增大一次聚光器的面積，但是一次聚光器面積的增加必然會導(dǎo)致聚光器系統(tǒng)質(zhì)量的增加，所以必須尋找這個折中點(diǎn)。本文采用先進(jìn)遺傳算法，在推進(jìn)系統(tǒng)能夠滿足衛(wèi)星設(shè)計功能的基礎(chǔ)上，對太陽能熱推進(jìn)系統(tǒng)參數(shù)進(jìn)行優(yōu)化分析，使光收集效率較高和聚光器系統(tǒng)質(zhì)量較小。

1.1 光的收集效率

太陽能熱推進(jìn)系統(tǒng)的光線的傳播過程如圖2所示，將收集到的太陽輻射能轉(zhuǎn)化為工質(zhì)氣體的動能，其過程首先要經(jīng)過一次聚光器的反射。因為對于聚光器不同的反射積層，不同的反射材料的反射率是不同的，對太陽輻射的反射并不能達(dá)到100%，一定會有一部分損失。通常采用現(xiàn)代鍍膜技術(shù)，反射率可高達(dá)90%～98%，稱為一次聚光器的反射率ηc1；經(jīng)一次聚光器改變傳播方向的光線匯聚到二次聚光器的球形入口處，此過程又會有一部分損失，即反射光到達(dá)二次聚光器入口表面的百分比一般在88%～96%之間，稱為主聚光器的截取因子ηi1；在二次聚光器內(nèi)全反射，最后被吸收器/推力室內(nèi)表面吸收的光線又會有一部分損失，吸收效率一般在87%～96%之間[3]，稱為折射式二次聚光器的效率ηc2；吸收器/推力室內(nèi)表面吸收涂層一般在0.8以上[4]，用ηα表示。

[5]，太陽光的收集效率Ceff為經(jīng)過上述過程損失所剩光能乘以系數(shù)cos4(φ/2)，其表達(dá)式為:

Ceff=ηc1·ηi1·ηc2·ηαcos4(φ/2)

(1)

圖2 太陽能熱微推進(jìn)系統(tǒng)光線傳播路線圖Fig.2 Light path diagram of solar thermal micro thrust system

1.2 一次聚光器質(zhì)量

一次聚光器主要由底板層、反射層和防護(hù)層組成[6-7]。聚光器的襯底材料選擇厚度為0.5 mm的聚亞酰胺薄膜，密度為1.44 g/cm3，此材料在航天工程中的多次使用證明其具有力學(xué)性能好、抗張強(qiáng)度高、不易產(chǎn)生褶皺和化學(xué)穩(wěn)定性好等優(yōu)點(diǎn)。反射膜選擇成膜性好、全太陽光譜范圍內(nèi)反射率高的鋁膜，厚度為10 μm，密度為2.7 g/cm3。為了提高反射率，在反射膜上面鍍一層厚度為120 nm的Al-Ag復(fù)合膜作為增反膜。為了防止反射膜及增反膜在太空中的蒸發(fā)，在增反膜上面鍍一層厚度為275 nm的SiO2薄膜，因為納米級的厚度太小，故忽略增反膜和防氧化膜的質(zhì)量，所以一次聚光器的質(zhì)量M=S(h1ρ1+h2ρ2)。

2 聚光器數(shù)學(xué)模型

2.1 基本模型

本文聚光器系統(tǒng)基準(zhǔn)構(gòu)型一次聚光器的聚光比為1500，入射半角Φ為45°，由此可以推導(dǎo)出一次聚光器拋物線的方程為y2=2.8z，光的收集效率Ceff.ref=52.72%，一次聚光器的質(zhì)量為Mref=0.82 kg。

在太陽能熱推進(jìn)系統(tǒng)中，收集到的太陽輻射最終將轉(zhuǎn)化為工質(zhì)氣體的動能，可寫出如式(2)所示關(guān)系式：

(2)

2.2 聚光器模型

在滿足微小衛(wèi)星任務(wù)需求的情況下，要求太陽能熱推進(jìn)系統(tǒng)的推力不小于0.1 N，比沖不小于800 s，建模如下：

選取目標(biāo)函數(shù)一次聚光器光的收集效率Ceff和質(zhì)量M，為簡化計算假設(shè)光收集效率中ηc1,ηi1,ηc2,ηα為固定值，其表達(dá)式分別為：

Ceff=0.723 6cos4(φ/2)

(3)

(4)

其中：φ是一次聚光器拋物面的位置角；f是一次聚光器的焦距；R為一次聚光器旋轉(zhuǎn)半徑，表達(dá)式為：

設(shè)計變量選取如下：

1)拋物面位置角φ：如圖1所示，φ表示入射光的反射光線與光軸的夾角，其范圍取20°～60°。

2)一次拋物面聚光器的焦距f與位置角φ和半徑R的關(guān)系密切，也是影響聚光器性能的主要因素，其范圍取0.1～1 m。

約束條件：g1表示推力至少為0.1 N，g2表示比沖至少為800 s。

得到一次聚光器光的收集效率Ceff隨接受半角φ和焦距f變化的響應(yīng)曲面如圖3所示，質(zhì)量M隨接受半角φ和焦距f變化的響應(yīng)曲面如圖4所示。

圖3 Ceff隨接受半角φ和焦距f變化的響應(yīng)圖Fig.3 Response diagram about Ceff with half-angle φ and focus f

圖4 M隨接受半角φ和焦距f變化的響應(yīng)圖Fig.4 Response diagram about mass M with half-angle φ and focus f

3 一次聚光器優(yōu)化分析

遺傳算法(Genetic Algorithm，GA)是模擬生物在自然環(huán)境中的遺傳和進(jìn)化過程而形成的一種自適應(yīng)全局優(yōu)化概率搜索算法，由美國Michigan大學(xué)的holland教授于1975年首次提出，源于達(dá)爾文的進(jìn)化論、孟德爾的群體遺傳學(xué)說和魏茨曼的五種選擇說。其基本思想是模擬自然界遺傳機(jī)制和生物進(jìn)化論而形成的一種過程搜索最優(yōu)解的算法。

圖5 遺傳算法流程圖Fig.5 Genetic algorithm flow chart

遺傳算法是模擬生物在自然環(huán)境中優(yōu)勝劣汰、適者生存的遺傳進(jìn)化過程而形成的一種具有自適應(yīng)能力的全局性惡概率搜索算法。它是從代表問題可能潛在階級的一個種群開始，首先將表現(xiàn)型映射到基因型及編碼，從而將解空間映射到編碼空間，每個編碼對應(yīng)問題的一個解，稱為染色體或個體。初始種群產(chǎn)生之后，按照適者生存和優(yōu)勝劣汰的原理，逐代演化生成出越來越好的近似解。在每一代，根據(jù)問題域中個體的適應(yīng)度大小選擇個體，并借助自然遺傳學(xué)的遺傳算子進(jìn)行組合交叉和變異，產(chǎn)生出代表新的解集的種群。這個過程使種群像自然進(jìn)化一樣，后代種群比前代更加適應(yīng)于環(huán)境，末代種群中的最優(yōu)個體經(jīng)過解碼可以作為問題近似最優(yōu)解，遺傳算法流程如圖5所示。目前在遺傳算法應(yīng)用中存在收斂速度慢、局部搜索能力差和早熟等問題，本文應(yīng)用改進(jìn)的遺傳算法NCGA和NSGA-II對問題進(jìn)行優(yōu)化分析。

本文研究的多目標(biāo)遺傳算法NCGA是由最早的GA發(fā)展而來，它視各目標(biāo)同等重要，通過排序后分組進(jìn)行交叉的方法實現(xiàn)“相鄰繁殖”的機(jī)制，從而使接近于Pareto前沿的解進(jìn)行交叉繁殖的概率增大，加速了計算收斂過程。NCGA步驟如下：

1)初始化：令t=0，設(shè)置第一代個體P0，種群數(shù)為N，計算個體對應(yīng)的適應(yīng)值函數(shù)，記為A。

2)令t=t+1，Pt=At-1。

3)排序：個體Pt按向聚集的目標(biāo)值的方向進(jìn)行排序。

4)分組：個體Pt根據(jù)上述排序分為若干組，每組由兩個個體組成。

5)交叉和變異：在每一組中執(zhí)行交叉和變異操作，由兩個父代個體產(chǎn)生兩個子代個體，同時父代個體被刪除。

6)重組：所有子代個體組成一組新的Pt。

7)更新：將At-1與Pt組合，按環(huán)境選擇(environment selection)機(jī)制，從2N個個體中，選出其中的N個個體。

8)終止：如果滿足終止條件，則終止優(yōu)化程序，否則返回至第2步。

根據(jù)工程經(jīng)驗，需要既考慮高聚光效率，又盡量減小一次聚光器質(zhì)量，因此設(shè)定目標(biāo)函數(shù)為Y=λ1Ceff-λ2M(其中λ1,λ2為聚光效率與聚光器質(zhì)量的加權(quán)值，λ1=0.6,λ2=0.4)。采用NCGA多目標(biāo)遺傳算法進(jìn)行優(yōu)化，種群大小取10，遺傳代數(shù)取30。結(jié)果如表1所示。

表1 NCGA優(yōu)化結(jié)果

表1中每一組均為通過遺傳算法計算得到在一定入射角和焦距的條件下的最優(yōu)解(最小質(zhì)量和最高吸收效率)，改變?nèi)肷浣桥c焦距后得到10組數(shù)據(jù)，結(jié)果表明經(jīng)過優(yōu)化后的聚光器系統(tǒng)，光的收集效率較優(yōu)化前大幅提高，聚光系統(tǒng)的質(zhì)量也大幅減小，達(dá)到了優(yōu)化目的。

但是，當(dāng)入射角φ太小時，一次聚光器的聚光比太大，焦斑的理論溫度將達(dá)到5000 K以上，任何材料承受不了這樣的高溫。取聚光器系統(tǒng)安裝半角φ=28°，焦距f=0.775 m,此時光的收集效率Ceff=64.1%，質(zhì)量M=0.359 kg。

4 結(jié)論

以太陽能熱推進(jìn)系統(tǒng)一次聚光器為研究對象，開展了對一次聚光器結(jié)構(gòu)參數(shù)的優(yōu)化及性能研究。主要工作與結(jié)果如下：

1)建立了一次聚光器的性能評估參數(shù)，分析了焦距和接收半角對光的收集效率及聚光系統(tǒng)質(zhì)量的影響；

2)對聚光器系統(tǒng)進(jìn)行建模，在基準(zhǔn)構(gòu)型基礎(chǔ)上采用NCGA進(jìn)行優(yōu)化，優(yōu)化后光的收集效率提高了21.59%，質(zhì)量減輕了0.461 kg，因此采用的優(yōu)化策略有效。

[1] 馬寶珊. 世界太陽能、風(fēng)能、生物能的發(fā)展現(xiàn)狀和前景[J]. 黑龍江電力技術(shù)， 1997, 19(6): 380-384.

MA Baoshan. The development status and prospects of solar energy, wind energy and bio-energy in the world[J]. Heilongjiang Electric Power， 1997, 19(6): 380-384.(in Chinese)

[2] Feuermann D, Gordon J M, Ries H. High-flux solar concentration with imaging designs[J]. Solar Energy, 1999, 65(2): 83-89.

[3] Wong W A, Geng S M, Castle C H, et al. Design, fabrication and test of a high efficiency refractive secondary concentrator for solar applications[C]//Proceedings of Energy Conversion Engineering Conference and Exhibit, 2000: 1117-1125.

[4] Nakamura T, Sullivan D, Mcclanahan J, et al. Solar thermal propulsion for small spacecraft[C]// Proceedings of AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, 2013.

[5] Feuermann D, Gordon J M. Solar fiber-optic mini-dishes: a new approach to the efficient collection of sunlight[J]. Solar Energy, 1999, 65(3): 159-170.

[6] 李金麗， 劉全校， 許文才. 光學(xué)薄膜及其發(fā)展現(xiàn)狀[J]. 化工新型材料, 2012, 40(2): 14-15.

LI Jinli, LIU Quanxiao, XU Wencai. The development of optical thin film[J].New Material, 2012, 40(2): 14-15.(in Chinese)

[7] 趙秀琴. 增透膜和增反膜[J]. 太原師范學(xué)院學(xué)報(自然科學(xué)版), 2003, 2(4): 42-45.

ZHAO Xiuqin.AR coating and high reflecting film[J].Journal of Taiyuan Teachers College(Natural Science Edition), 2003, 2(4): 42-45.(in Chinese)

Condenser parameters optimization in solar thermal propulsion system

HUANG Minchao, WU Shang

(College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China)

An advanced genetic algorithm was used to make an optimization analysis on the solar thermal propulsion system condenser parameters for reducing the mass of system. Considering the sunlight collection efficiency and the mass of condenser as the optimization objective function, a mathematical model of condenser work was established and several related simulation researches were carried out. The simulation results show that the advanced genetic algorithm can be effectively used in the optimization analysis of the sunlight collection efficiency and the condenser mass.

condenser; collection efficiency; mass; optimization

??(

)

10.11887/j.cn.201605015

http://journal.nudt.edu.cn

2015-12-25

國家部委級基金資助項目(9140A20100413KG01293)

黃敏超(1969—)，男，重慶榮昌人，副教授，博士，碩士生導(dǎo)師，E-mail:h_mchao01@qq.com

V439

A

1001-2486(2016)05-094-05