許 勇,黃 勇,余永剛

(中國空氣動力研究與發(fā)展中心 計算空氣動力學(xué)研究所,四川 綿陽 621000)

電大尺寸目標(biāo)電磁散射的并行FVTD計算

許 勇*,黃 勇,余永剛

(中國空氣動力研究與發(fā)展中心 計算空氣動力學(xué)研究所,四川 綿陽 621000)

實際目標(biāo)的高頻電磁散射模擬屬于電大尺寸目標(biāo)電磁計算問題，其高精度數(shù)值計算通常伴隨著大存儲量和大計算量的沉重負(fù)擔(dān)。為保證高精度數(shù)值計算，構(gòu)建了直接求解電磁學(xué)麥克斯韋方程組，采用貼體多塊結(jié)構(gòu)網(wǎng)格的時域有限體積法(FVTD)解算器。為解決大規(guī)模網(wǎng)格帶來的大計算量問題，則采用MPI并行編程，進(jìn)行網(wǎng)格多進(jìn)程分割、負(fù)載平衡以及進(jìn)程通信設(shè)置和程序并行化處理，成功地對電大尺寸飛翼外形進(jìn)行了L波段雙站電磁散射場和雷達(dá)截面(RCS)計算。結(jié)果表明pmbRCS3d這一并行高精度電磁散射模擬軟件具有穩(wěn)定和魯棒特性，適合進(jìn)一步應(yīng)用于目標(biāo)更高頻段電磁散射計算。

麥克斯韋方程組；時域有限體積法(FVTD)；并行計算；L波段；雷達(dá)截面(RCS)

0 引 言

飛行器氣動/隱身優(yōu)化設(shè)計需要發(fā)展快速高效、高精度電磁數(shù)值方法，但是復(fù)雜外形目標(biāo)的電磁散射包含了鏡面散射、邊緣繞射、爬行波等多種復(fù)雜現(xiàn)象，準(zhǔn)確模擬非常困難。

高頻漸進(jìn)等傳統(tǒng)工程方法雖然快捷，能捕獲主要散射特征，但誤差較大，不能滿足隱身設(shè)計的精細(xì)要求。隨著計算機(jī)能力的提高，國內(nèi)外高精度數(shù)值方法獲得巨大發(fā)展，主要有計算微分方程的時域有限差分FDTD[1-3]方法，和計算積分方程的多層快多極子方法[4-7]。其中直接求解電磁學(xué)麥克斯韋方程組的時域微分方法，最早來源于Yee K S[1]提出的FDTD，它基于蛙跳格式，電場和磁場在網(wǎng)格交錯配置以彌補格式黏性，有易編程和整個散射空間可統(tǒng)一計算的優(yōu)點，但由于使用笛卡兒直角網(wǎng)格，模擬物面存在階梯效應(yīng)。20世紀(jì)90年代以來，基于流體力學(xué)的歐拉方程組和電磁學(xué)的麥克斯韋方程組具有共同的雙曲型數(shù)學(xué)特征，人們開始在計算電磁學(xué)(CEM)引入計算流體力學(xué)(CFD)技術(shù)[8-15]，比較典型的是采用貼體坐標(biāo)系的時域有限體積法(FVTD)。

FVTD解算器有以下特點：直接求解麥克斯韋方程組；在寬頻入射波的一個計算過程中可得到多個頻率的結(jié)果；能達(dá)到接近矩量法(MOM)的數(shù)值計算精度；存儲僅與未知量數(shù)目同數(shù)量級，能不需特殊處理地模擬散射、腔激勵等復(fù)雜現(xiàn)象。但因涉及波在空間長時間傳播計算，以及抑制耗散和色散要有一定的網(wǎng)格密度，網(wǎng)格規(guī)模會隨頻率增加。

實際飛機(jī)等目標(biāo)在高頻區(qū)為電大尺寸，例如B2隱形轟炸機(jī)，長度和翼展達(dá)數(shù)十米，對應(yīng)高頻區(qū)(X波段)波長最大到1000量級，每個波長取10～20個網(wǎng)格點，則整個規(guī)模高達(dá)10億量級。超大計算規(guī)模成為制約氣動/隱身優(yōu)化的瓶頸。目前解決途徑是在高性能計算機(jī)集群上，采用大規(guī)模并行計算技術(shù)。

為解決電大目標(biāo)高頻電磁散射大規(guī)模計算，本文通過空間網(wǎng)格塊分割、負(fù)載平衡、多進(jìn)程MPI通信等，建立了多塊結(jié)構(gòu)網(wǎng)格FVTD并行程序pmbRCS3d，在與金屬球RCS解析解驗證對比基礎(chǔ)上，計算了電大尺寸飛翼外形L波段雙站電磁散射場和雷達(dá)截面(RCS)，結(jié)果表明了并行計算程序的正確性和魯棒性，利于向更高頻段(X波段)擴(kuò)展。

1 數(shù)值方法

1.1 控制方程

麥克斯韋方程組的兩個時變旋度微分方程，法拉第電磁感應(yīng)定律和安培環(huán)路定理，在自由空間、直角坐標(biāo)系下可表示為守恒形式：

式中:B是磁感應(yīng)強度矢量，H是磁場強度，D是電位移矢量，E是電場強度矢量。

對于復(fù)雜外形物體，經(jīng)坐標(biāo)變換，

得到曲線坐標(biāo)系下麥克斯韋方程組守恒形：

式中，V是坐標(biāo)變換的雅可比矩陣行列式值。

1.2 積分方法

采用散射場變量既可保持守恒形式不變，又能避免傳播入射場的數(shù)值計算誤差，入射場解析給出。有限體積法能保持整個網(wǎng)格空間的通量守恒，對式(4)在網(wǎng)格單元作時、空積分，得到數(shù)值離散形式：

有限體積法的關(guān)鍵是精確模擬原變量Q在單元分界面的狀態(tài)變量，以得到相應(yīng)精確的分界面流通量F，本文F采用STEGER-WARMING分裂：

S、S-為相似矩陣，Λ+、Λ-分別為正負(fù)特征值構(gòu)成的對角矩陣，QL、QR代表分界面處左右狀態(tài)變量，采用MUSCL格式插值獲得，時間推進(jìn)則采用4步龍格-庫塔法。

1.3 入射電磁波波形

入射電磁波選取簡諧波(A0是波幅，c為波速)：

1.4 邊界條件

完全導(dǎo)電壁處的邊界條件：

截斷的網(wǎng)格空間外邊界相容條件：

1.5 近場到遠(yuǎn)場的轉(zhuǎn)換

穩(wěn)定的時變場在頻域收斂后，采用格林定理：

R是觀察點，R′是虛擬積分面上點矢量，r是觀察方向單位矢量，A是虛擬積分面，k是自由空間波數(shù)，j是虛數(shù)單位。

1.6 雷達(dá)散射截面計算

1.7 并行計算

電大尺寸電磁散射伴隨大規(guī)模三維網(wǎng)格，須用多進(jìn)程并行計算，包括：(1) 多進(jìn)程分割和負(fù)載平衡；(2) 程序并行化處理，采用MPI接口進(jìn)行通信。

2 計算結(jié)果與分析

2.1 并行化多塊FVTD程序pmbRCS3d驗證

多塊網(wǎng)格FVTD程序mbRCS3d，并行化擴(kuò)充為pmbRCS3d。首先選擇帶解析解的金屬球(ka=10，其中k為波數(shù)，a是半徑)為驗證目標(biāo)。

圖1是網(wǎng)格塊經(jīng)8個進(jìn)程分割后表面網(wǎng)格。圖2是表面散射場等值線云圖。圖3是表面誘導(dǎo)電流分布，結(jié)果顯示前向散射最強。圖4是雙站RCS數(shù)值與精確解析級數(shù)解比較，吻合很好，說明整個并行FVTD算法具有很好的精度。

2.2 電大尺寸實際外形RCS計算

針對飛翼外形(類似X-47)的L波段(f=1.2 GHz)電磁散射問題，長度11.64 m，翼展18.9 m，長度對應(yīng)電尺寸為47個波長，翼展對應(yīng)76個波長。采用36塊網(wǎng)格離散整個計算空間，網(wǎng)格點16 426 816個。網(wǎng)格密度根據(jù)電磁波頻率選取，要求每個坐標(biāo)方向上一個波長(250 mm)內(nèi)至少布15～20個網(wǎng)格點，散射體表面網(wǎng)格加密，遠(yuǎn)離物面區(qū)域適當(dāng)疏松，以便減少計算量，圖5是重新分塊后的物面網(wǎng)格邊界。在集群上使用16個進(jìn)程計算，耗時420 h。

該大網(wǎng)格量電磁場并行FVTD計算，得到穩(wěn)定振蕩的簡諧場和逐步收斂的頻域散射場，從圖6和圖7散射場云圖能看出散射波耦合強烈區(qū)域。圖8的結(jié)果(其中φ為與+X軸所成方位角)表明：前向RCS較大，后向RCS較小在-35 dB附近，有4個較強散射方位(分別對應(yīng)4個垂直于機(jī)身前緣和機(jī)翼前緣的方位)。

3 結(jié) 論

通過散射電磁場的FVTD解算器和基于MPI的并行編程，一定程度解決了電大尺寸目標(biāo)高頻散射大計算量難題。計算和驗證表明：該系統(tǒng)能提供高精度的電磁場和RCS分布，為氣動/隱身優(yōu)化提供精細(xì)工具,并有利于向更高頻段擴(kuò)展。

[1]Yee K S.Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media[J].IEEE Trans.Antennas and Propagation,1966,14(3):302-307.

[2]Holland R.Finite-difference solution of Maxwell’s equations in generalized coordinates[J].IEEE Trans.Nuclear Science,2007,30(6):4589-4591.

[3]Fusco M.FDTD algorithm in curvilinear coordinates[J].IEEE Transactions on Antennas and Propagation,1990,38(1):76-89.

[4]Song J M,Chew W C.Fast multipole method solution of combined field integral equation[J].IEEE Trans.on Antennas and Propagation,1994,7(16):760-765.

[5]Nachman A.A brief perspective on computational electromagnetic[J].Journal of Computational Physics,1996,125:237-239.

[6]V elamparambil S,Chew W C,Song J M.10 million unknowns:is it that big?[J].IEEE Antennas and Propagation Magazine,2003,45(2):43-58.

[7]Pan X M,Sheng X Q.A sophisticated parallel MLFMA for scattering by extremely large targets[J].IEEE Antennas and Propagation Magazine,2008,50(3):129-138.

[8]Shankar V,Hill W,Mohammadian A H.A CFD-based finite-volume procedure for computational electromagnetics inter-disiplinary applications of CFD methods[R].A1AA CP 89-1987.

[9]Shankar V,Hill W,Mohammadian A H.A time-domain differential solver for electromagnetic scattering problems[J].Proc.IEEE,1989,77(5):709-721.

[10]Ahuja V,Long L N.A parallel finite-volume Runge-Kutta algorithm for electromagnetic scattering[J].Journal of Computational Physics,1997,137:299-320.

[11]Kabakian A V.A three-dimensional spectral collocation time- domain solver for electromagnetic wave scattering[R].AIAA 98-0980.

[12]Shang J S,Gaitonded.Characteristic-based time-dependent Maxwell equations solver on a general curvilinear frame[R].AIAA 93-3178,1993.

[13]許勇,樂嘉陵.基于CFD的電磁散射對數(shù)值模擬[J].空氣動力學(xué)學(xué)報,2004,22(2):185-189.

[14]Shang J S.Time-domain electromagnetic scattering simulations on multicomputers[J].Journal of Computational Physics,1996,128:381-390.

[15]Shang J S,Gaitonde D.Characteristic-based time-dependent Maxwell equations solvers on a general curvilinear frame[J].AIAA Journal,1995,33(3):491-498.

ParellelFVTDcomputationforelectromagneticscatteringofelectricallylargeobjects

XU Yong*,HUANG Yong,YU Yonggang

(ComputationalAerodynamicsInstituteofChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China)

Numerical simulation of real military objects′ electromagnetic scattering at high frequency bands is ordinary corresponding to the electromagnetic field calculation for electrically large objects.This simulation can be computationally expensive with large storage and calculation works when high precision numerical algorithm is used.To ensure high computational precision,we programed a finite volume time domain (FVTD) solver pmbRCS3d based on body-fitted multi-blocks grids,which computes the Maxwell equations of electromagnetism directly.In order to handle the massive calculation work due to huge number of grids,a parallel MPI programing technique was used together with careful design of the grid partition,load balance,data communication,and parallelization.The pmbRCS3d software was applied successfully to calculate the electromagnetic scattering field and the bi-station RCS of an L-band electrically large flying-wing shape.The results show that the developed pmbRCS3d software has stable and robust characteristic and can be used for the RCS calculation of higher frequency band objects.

Maxwell equations; finite volume time domain method; parallel computing; L frequency band; radar cross section

0258-1825(2017)06-0797-04

O441.4

A

10.7638/kqdlxxb-2015.0071

2015-06-03;

2015-07-24

裝備預(yù)先研究項目(51313010303)

許勇*(1971-),男，四川彭州人，副研究員，研究方向：計算電磁學(xué).E-mail：stephen000@sina.com

許勇,黃勇,余永剛.電大尺寸目標(biāo)電磁散射的并行FVTD計算[J].空氣動力學(xué)學(xué)報,2017,35(6):797-800.

10.7638/kqdlxxb-2015.0071 XU Y,HUANG Y,YU Y G.Parellel FVTD computation for electromagnetic scattering of electrically large objects[J].Acta Aerodynamica Sinica,2017,35(6):797-800.