孔蓓蓓 盛新慶

(北京理工大學信息與電子學院電磁仿真中心,北京 100081)

?

利用近場Gordon積分近似的矩量法-物理光學混合方法

孔蓓蓓 盛新慶

(北京理工大學信息與電子學院電磁仿真中心,北京 100081)

針對電大金屬目標的電磁計算,提出了一種新的矩量法-物理光學(Method of Moment-Physical Optics,MoM-PO)混合方法,以解決傳統(tǒng)的MoM-PO混合法中PO區(qū)域和MoM區(qū)域耦合項的計算時間過長問題.用Gordon積分計算目標的PO區(qū)域對MoM區(qū)域的作用,并加入近場近似處理.該方法避免了傳統(tǒng)的MoM-PO混合法中耦合項積分方程的計算,加大PO區(qū)域的剖分面元的大小,能夠有效地降低矩陣規(guī)模和未知數個數,因而降低了內存,減少了計算時間.數值算例結果表明,近場Gordon積分近似的MoM-PO混合方法能夠有效地減少耦合項的計算時間,并能達到理想的精度.

MoM;PO;混合方法;Gordon積分

DOI 10.13443/j.cjors.2015110603

引 言

計算電大導體目標的電磁散射是計算電磁學中一個重要的和具有挑戰(zhàn)性的問題[1].矩量法(Method of Moment,MoM)由于其精確性和高效性,成為計算電磁學的主要方法之一,在計算三維目標的散射問題上得到了廣泛應用.用矩量法求解表面積分方程時,需要采用適當的基函數(如RWG(Rao-Wilton-Glisson)基函數[2])對等效源進行離散,為了達到理想的精度,對三維目標的剖分通常采用較小的網格.隨著目標電尺寸的增大,矩量法的計算成本大大增加,無法滿足計算需求.

物理光學(Physical Optics,PO)方法[3]非常適合求解電大問題,然而,對于目標中較為復雜和細小的部分,物理光學法不能準確地近似表面電流,為此通常采用與矩量法相結合的方法[4-9].尤其對于天線的輻射問題,MoM-PO混合法十分有效.此外,由于PO方法對于目標陰影區(qū)無法準確地近似,對曲率較小的電大目標的散射,也可以采用MoM-PO混合法.MoM-PO混合法將目標分為PO區(qū)域和MoM區(qū)域,PO區(qū)域的等效電流用PO電流近似,MoM區(qū)域的等效電流仍然采用離散積分方程、求解矩陣的方法進行計算.但是,還需要計算PO區(qū)域的電流對MoM區(qū)域的作用,而這部分仍然需要較長的時間[10],因為兩個區(qū)域的基函數之間的耦合都需要進行計算.

為了減少耦合項的計算時間,本文提出了一種利用近場Gordon積分近似的MoM-PO混合方法.PO區(qū)域對MoM區(qū)域的作用直接通過求解PO電流在MoM區(qū)域的散射場得到,散射場的求解采用Gordon積分方法[11],同時加入近場近似[12]以減少近場散射場的相位誤差.由于該方法不需要用基函數表示PO電流,對PO區(qū)域可以采用較大的面元進行剖分,計算耦合項時也不需要求積分方程,從而大大減少了耦合項的計算時間.

1 利用近場Gordon積分近似的MoM-PO方法

1.1 傳統(tǒng)的MoM-PO混合法

對于一個平面波入射下的完美電導體(Perfect Electric Conductor,PEC)目標的散射問題,入射場用(Ei,Hi)來表示,以電場積分方程為例,可以建立如下方程[13]:

(1)

L算子定義為

(2)

(3)

r代表場點位置,r′代表源點位置.

用RWG基函數對面積分方程(1)進行離散,gi為RWG矢量基函數,并采用伽遼金匹配,可得到下面的離散方程:

(4)

將式(4)簡寫為[Pij]{J}={bi},其中J為待求未知數.傳統(tǒng)的MoM-PO混合法將積分區(qū)域分為兩部分:PO區(qū)域和MoM區(qū)域.PO區(qū)域為目標的光照區(qū),該區(qū)域的電流可以近似地表示為

(5)

NPO為PO區(qū)域未知數的個數．PO區(qū)域的J可由式(5)得到,從而,待求未知數僅為MoM區(qū)域的未知數,求解的離散方程為

(6)

(7)

1.2 利用近場Gordon積分近似的MoM-PO方法

MoM-PO混合法使未知數的個數和矩陣規(guī)模都得到了縮減,降低了內存并減少了計算時間.然而,隨著計算目標的增大,PO區(qū)域與MoM區(qū)域的耦合項b′的計算時間也大大增加.由式(7)可以看出,b′的計算時間主要受PO區(qū)域離散未知數的個數NPO以及每一個離散位置對MoM區(qū)域的作用(即積分項PijJj)的計算時間影響.為了減少耦合項的計算時間,我們提出了利用Gordon積分和近場近似的MoM-PO混合方法.

(8)

(9)

圖1 三角面元的符號

耦合項b′為PO區(qū)域的等效電流JPO在MoM區(qū)域上的散射場與試函數的內積:

b′=∫gi·Esds.

(10)

與式(7)相比,式(10)不需要計算積分方程,從而不需要用基函數表示PO電流,對PO區(qū)域的剖分可以選擇較大的面元,大大減少了計算時間.

另外,我們將多層快速多極子技術(Multilevel Fast Multipole Algorithm,MLFMA)應用于MoM中,加速對MoM區(qū)域矩陣填充和未知數的求解,對于計算電大尺寸的目標,可以節(jié)約計算時間.在接下來的數值算例中,混合法的MoM區(qū)域未知數的計算都采用MLFMA技術進行加速.

2 數值算例

在本節(jié)中,我們分別用MLFMA、傳統(tǒng)混合法以及本文中提出的近場Gordon積分近似的混合方法計算了兩個算例.根據入射場的入射方向,將目標的光照區(qū)作為PO區(qū)域,陰影區(qū)作為MoM區(qū)域.對于近場Gordon積分近似的混合方法,我們采用平均邊長為0.5λ的三角面元對PO區(qū)域進行剖分,MoM區(qū)域的剖分大小為0.1λ,兩個區(qū)域之間采用共形網格.MLFMA和傳統(tǒng)混合法均采用邊長為0.1λ的網格剖分.

首先計算了一個底面積為8λ×8λ,高為5λ的PEC立方體.平面波沿著-Z軸入射,VV極化下的雙站RCS計算結果如圖2所示.表1統(tǒng)計了不同方法計算的未知數個數以及耦合項的計算時間,為了比較兩種混合法的計算精度,引入RCS的均方根誤差:

(11)

式中:σM代表MLFMA的計算結果,作為精度比較的參考值;σH代表混合法的計算結果;N代表觀察角度的個數.從仿真結果可以看出,近場Gordon積分近似的混合方法計算的RCS結果能夠較為準確地描繪RCS的變化趨勢,與傳統(tǒng)的混合法相比,耦合項的計算時間大大減少,減少了約99%.這是由于在計算耦合項的過程中,我們的方法不需要計算積分方程,每個亮區(qū)面元的散射場可以由面元的頂點位置解析計算出來,從而節(jié)省了時間.加大亮區(qū)面元的剖分,可以大大減少亮區(qū)面元散射場計算個數,使耦合項計算時間減少.另外,由于兩個區(qū)域之間采用共形網格,使得MoM區(qū)域的未知數個數有所下降.

(a) 模型

(b) 雙站RCS圖2 大小為8λ×8λ×5λ的立方體模型及其RCS仿真

未知數個數耦合項計算時間/sεRMS/dB近場Gordon近似混合法414531374.7393傳統(tǒng)混合法77104150252.7498

由于本文所提方法不僅考慮了亮區(qū)電流的等效,還在計算耦合項時采用了遠場散射場的計算方法,比起傳統(tǒng)的用RWG將PO電流用基函數表示的方法誤差有所增大,對于圖2中的算例,我們的方法誤差上升了2 dB.

同樣計算了一個半徑為3λ的PEC半球,平面波垂直于半球平面入射.計算結果如圖3和表2所示,與傳統(tǒng)方法相比,MoM區(qū)域未知數個數幾乎相同,耦合項的計算時間明顯減少,εRMS上升了2.6 dB.

(a) 模型

(b) 雙站RCS圖3 半徑為3λ的半球模型及其RCS仿真

未知數個數耦合項計算時間/sεRMS/dB近場Gordon近似混合法201805764.5167傳統(tǒng)混合法2030017091.9029

3 結 論

本文提出了一種新的MoM-PO混合方法,用Gordon積分計算散射體亮區(qū)對暗區(qū)的作用,并加入近場近似處理.該方法避免了耦合項計算過程中積分方程的計算,因而可以在PO區(qū)域采用較大的面元進行剖分,有效地降低了耦合項的計算時間.從仿真結果來看,近場Gordon積分近似的MoM-PO混合方法能夠準確地計算目標的散射,與傳統(tǒng)的混合法相比,誤差增加范圍不超過3 dB.

[1] 潘小敏, 盛新慶. 電特大復雜目標電磁特性的高效精確并行計算[J].電波科學學報, 2008, 23(5): 888-893.

PAN X M, SHENG X Q. Efficient and accurate parallel computation of electromagnetic scattering by extremely large targets [J]. Chinese journal of radio science, 2008, 23(5): 888-893. (in Chinese)

[2] RAO S M, WILTON D R, GLISSON A W. Electromagnetic scattering by surfaces of arbitrary shape[J]. IEEE transactions on antennas and propagation, 1982, 30(3): 409-412.

[3] 陳博韜, 雷振亞, 謝擁軍,等. 基于改進物理光學法的電大目標雙站RCS的預估[J]. 電波科學學報, 2010, 25(5): 960-965.

CHEN B T, LEI Z Y, XIE Y J, et al. Bistatic RCS Prediction with improved PO method for electrically large targets[J]. Chinese journal of radio science, 2010, 25(5): 960-965. (in Chinese)

[4] Medgyesi-Mitschang L, WANG D S. Hybrid solutions for scattering from perfectly conducting bodies of revolution[J]. IEEE transactions on antennas and propagation, 1983, 31(4): 570-583.

[5] JAKOBUS U, LANDSTORFER F M. Improved PO-MM hybrid formulation for scattering from three-dimensional perfectly conducting bodies of arbitrary shape[J]. IEEE transactions on antennas and propagation, 1995, 43(2): 162-169.

[6] JAKOBUS U, LANDSTORFER F M. Improvement of the PO-MoM hybrid method by accounting for effects of perfectly conducting wedges[J]. IEEE transactions on antennas and propagation, 1995, 43(10): 1123-1129.

[7] WEI X C, LI E P. Wide-band EMC analysis of on-platform antennas using impedance-matrix interpolation with the moment method-physical optics method[J]. IEEE transactions on electromagnetic compatibility, 2003, 45(3): 552-556.

[8] CHEN M, ZHANG Y, ZHAO X W, et al. Analysis of antenna around NURBS surface with hybrid MoM-PO technique[J]. IEEE transactions on antennas and propagation, 2007, 55(2): 407-413.

[9] 龔主前, 朱國強, 秦家銀.MM-PO混合法中一致性劈面繞射電流研究[J].電波科學學報, 2008, 23(3): 576-579.

GONG Z Q, ZHU G Q, QIN J Y. Uniform expressions for the diffracted currents on the wedge surface for the MM-PO technique[J]. Chinese journal of radio science, 2008, 23(3): 576-579. (in Chinese)

[10]LIU Z L, WANG C F. Efficient iterative Meth-od of moments-physical optics hybrid technique for electrically large objects[J]. IEEE transactions on antennas and propagation, 2012, 60(7): 3520-3525.

[11]GORDON W B. Far field approximations to the kirchhoff-helmholtz representation of scattered fields[J]. IEEE transactions on antennas and propagation, 1975, 23(7): 590-592.

[12]LEGAULT S R. Refining Physical Optics for Near-field Computations[J]. Electronics Letters, 2004,40(1)：71-72.

[13]盛新慶. 計算電磁學要論[M]. 北京: 科學出版社, 2004,16-72.

孔蓓蓓 (1990-),女,甘肅人,北京理工大學博士研究生,研究方向為電磁場與微波技術.

盛新慶 (1968-),男,安徽人,北京理工大學教授,博士生導師,2001年度中國科學院“百人計劃”入選者.2004年度教育部長江學者特聘教授.主要從事計算電磁學、目標電磁特性與探測技術、隱身目標分析與設計、天線分析與設計、電磁環(huán)境預測技術等方面的研究.

Efficient hybrid MoM-PO technique with near-field Gordon integral

KONG Beibei SHENG Xinqing

(CenterforElectromagneticSimulation,SchoolofInformationandElectronics,BeijingInstituteofTechnology,Beijing100081,China)

For the electromagnetic calculation of electrically large conducting object, a new hybrid method of moments-physical optics(MoM-PO) technique is presented to overcome the time-consuming in calculating the coupling term of PO region and MoM region with the increase of the target size by using the conventional hybrid MoM-PO method. The near-field scattering is calculated by Gordon integral to avoid the calculation of the integral equation of the coupling term of the conventional method and to increase the mesh size of PO region, which provides the possibility to reduce the computational time and computer memory. The results of numerical examples indicate the effectiveness of the new method. Compared with the traditional hybrid MoM-PO method, the computational time of the coupling term is reduced. Meanwhile, appropriate accuracy can be achieved.

MoM; PO; hybrid method; gordon integral

10.13443/j.cjors.2015110603

2015-11-06

國家自然科學基金(61371002)； 111引智計劃(B14010)

O441.4

A

1008-9268(2016)04-0786-05

孔蓓蓓, 盛新慶. 利用近場Gordon積分近似的矩量法-物理光學混合方法[J]. 電波科學學報，2016,31(4):786-790.

KONG B B, SHENG X Q. Efficient hybrid MoM-PO technique with near-field Gordon integral[J]. Chinese journal of radio science，2016,31(4):786-790. (in Chinese). DOI: 10.13443/j.cjors.2015110603

聯系人： 孔蓓蓓 E-mail: kongbeibei@bit.edu.cn