王飛魏兵

(西安電子科技大學(xué)物理系,西安 710071)

(2012年11月26日收到;2013年1月22日收到修改稿)

1 引言

鐵氧體材料為磁色散介質(zhì),在有外加磁場(chǎng)時(shí)又是磁各向異性介質(zhì).近年來(lái),對(duì)其物理特性的研究已成為熱點(diǎn)之一,如文獻(xiàn)[1,2]分別研究了磁化鐵氧體材料的自旋電流及自旋極化電流的影響和作用.磁化鐵氧體材料的磁導(dǎo)系數(shù)是頻率的張量函數(shù),對(duì)含有磁化鐵氧體目標(biāo)的電磁散射進(jìn)行嚴(yán)格三維分析是非常困難的.本文主要研究應(yīng)用于磁各向異性色散介質(zhì)的電磁散射的時(shí)域有限差分(finite difference time domain,FDTD)算法.文獻(xiàn)中處理各向同性色散介質(zhì)電磁問題的FDTD方法有：遞歸卷積(recursive convolution,RC)法[3]、分段線性遞歸卷積(piecewise linear recursive convolution,PLRC)法[4]、電流密度卷積(current density convolution,JEC)法[5]、分段線性電流密度遞歸卷積(piecewise linear current density recursive convolution,PLJERC)法[6,7]、輔助方程 (auxiliary differential equation,ADE)法[8,9]、Z 變換 (Z-transform,ZT)法[10-12]、移位算子(shift operator,SO)法[13]和半解析遞推卷積(semianalytical recursive convolution,SARC)法[14]等.外加磁場(chǎng)時(shí)鐵氧體介質(zhì)具有磁各向異性色散特性,因此對(duì)磁化鐵氧體問題,需要將上述算法進(jìn)行推廣.Chen等[15,16]計(jì)算分析了處于自由空間和波導(dǎo)內(nèi)的飽和磁化鐵氧體目標(biāo)的環(huán)流特征.楊利霞等分別采用Pade′-FDTD法[17],RC-FDTD法[18]和直接差分法[19]計(jì)算了飽和磁化鐵氧體球的電磁散射.

Z變換法是數(shù)字信號(hào)處理中處理離散信號(hào)的一種常用的變換方法[20].在FDTD法中運(yùn)用Z變換的優(yōu)點(diǎn)是將與頻率相關(guān)的媒質(zhì)的介電常數(shù)由頻域變換到Z域,把時(shí)域中的卷積積分轉(zhuǎn)化為Z域內(nèi)的乘法,可以方便地求出媒質(zhì)的差分方程,避免了時(shí)域卷積運(yùn)算.本文根據(jù)飽和磁化鐵氧體介質(zhì)電磁特性,將Z變換法推廣應(yīng)用于任意磁化方向的鐵氧體目標(biāo)FDTD分析.考慮到在實(shí)驗(yàn)室參照系中飽和磁化鐵氧體磁導(dǎo)系數(shù)張量的形式和磁化方向密切相關(guān),首先利用坐標(biāo)系轉(zhuǎn)換矩陣給出實(shí)驗(yàn)室系中飽和磁化鐵氧體的頻域磁化率張量,再采用部分分式展開方法通過(guò)Z變換得到磁化率張量的Z域形式,最后將飽和磁化鐵氧體的頻域磁本構(gòu)關(guān)系轉(zhuǎn)換到Z域再根據(jù)移位定理轉(zhuǎn)換到離散時(shí)域,從而實(shí)現(xiàn)電磁場(chǎng)量B→H的迭代計(jì)算.該方法對(duì)于研究鐵氧體目標(biāo)處于任意方向外磁場(chǎng)情況下的電磁散射問題具有應(yīng)用價(jià)值.作為驗(yàn)證,計(jì)算了飽和磁化鐵氧體球的后向雷達(dá)散射截面(radar scattering cross-section,RCS),所得結(jié)果與文獻(xiàn)[21]結(jié)果一致.理論推導(dǎo)及算例表明該方法是正確有效的.文中時(shí)諧因子取exp(jωt).

2 飽和磁化鐵氧體Maxwell旋度方程和頻域本構(gòu)關(guān)系

飽和磁化鐵氧體介質(zhì)的Maxwell旋度方程及頻域本構(gòu)關(guān)系為

其中μr(ω)是飽和磁化鐵氧體相對(duì)磁導(dǎo)系數(shù)張量的頻域形式.按標(biāo)準(zhǔn)Yee元胞[22]對(duì)旋度方程(1)和(2)式差分離散便可實(shí)現(xiàn)H→E→B的FDTD迭代計(jì)算,而要實(shí)現(xiàn)B→H的迭代計(jì)算則需要對(duì)頻域本構(gòu)關(guān)系(3)式進(jìn)行特殊處理.(3)式可以寫為另一種形式[23]

其中μ∞是無(wú)限大頻率時(shí)的相對(duì)磁導(dǎo)系數(shù)張量,M為磁化矢量

其中 χm(ω)為頻域磁化率張量,設(shè) χmαβ(ω)是χm(ω)的元素 (α,β=x,y,z),有

由(3)—(5)式可知

3 頻域磁化率張量從外磁場(chǎng)坐標(biāo)系到實(shí)驗(yàn)室坐標(biāo)系的轉(zhuǎn)換

各向異性介質(zhì)的本構(gòu)參數(shù)與所選取的坐標(biāo)系有關(guān).飽和磁化鐵氧體的磁化率張量在外磁場(chǎng)坐標(biāo)系 (設(shè)為 x′y′z′系)中有簡(jiǎn)單的表達(dá)式,而電磁測(cè)量往往在實(shí)驗(yàn)室坐標(biāo)系(設(shè)為xyz系)下進(jìn)行.下面根據(jù)外磁場(chǎng)系和實(shí)驗(yàn)室系的轉(zhuǎn)換矩陣以及外磁場(chǎng)系中飽和磁化鐵氧體的相對(duì)磁導(dǎo)系數(shù)張量,得到實(shí)驗(yàn)室系中磁化率張量的表達(dá)式.

圖 1 外磁場(chǎng) (x′y′z′)系和實(shí)驗(yàn)室 (xyz)系

建立外磁場(chǎng)直角坐標(biāo)系：x′y′z′系,并令外磁場(chǎng)H0方向?yàn)閦′軸正向,即H0=z?′H0,z?′為z′軸正向單位矢,H0是H0的大小.設(shè)外磁場(chǎng)H0在實(shí)驗(yàn)室直角坐標(biāo)系(xyz系)中的方向角為θt,φt,如圖1所示.

設(shè) x′y′z′系中,飽和磁化鐵氧體的相對(duì)磁導(dǎo)系數(shù)張量為 μ′r(ω),根據(jù)文獻(xiàn) [23]有

式中

上式中ω0=γH0,γ=1.76×107rad/Oe·s為旋磁比,ωm=γ4πMs,Ms為飽和磁化率,α為阻尼因子.

式可知

設(shè)μ′∞為單位矩陣,即

由 (8)—(11)式可得到 x′y′z′系中

其中

可以證明(12)式中

其中

xyz系中的 μr(ω)可由 x′y′z′系中的 μ′r(ω)通過(guò)轉(zhuǎn)換計(jì)算得到[24]

U 為x′y′z′系到 xyz系的轉(zhuǎn)換矩陣,UT是 U 的轉(zhuǎn)置矩陣,且

將(10)式代入(16)式有

即有

根據(jù)(12),(14),(16)式可知,χm(ω)的各元素也可寫成jω的有理分?jǐn)?shù)形式

其中qn的值見(15)式,pαβn的值可由(15)和(19)式計(jì)算得到, 即

4 實(shí)驗(yàn)室系中極化率張量的Z域形式

由上可知,χm(ω)各個(gè)元素 χmαβ(ω)可寫成jω的有理分式的真分式,可采用部分分式展開法展開為部分分式進(jìn)而轉(zhuǎn)換到Z域.令

由(15),(20)式知A(jω)有兩個(gè)共軛復(fù)數(shù)單根(2單極點(diǎn)),即A(jω)可寫為

2個(gè)共軛復(fù)數(shù)單根為

根據(jù)部分分式展開法,χmαβ(ω)可展開為部分分式

其中系數(shù)

上式的Z域形式為

其中

設(shè)其為張量χ1和χ2的元素,在(29),(30)式中同時(shí)考慮了(24)式.

5 飽和磁化鐵氧體的Z變換——FDTD方法

飽和磁化鐵氧體頻域本構(gòu)關(guān)系(4)式和磁化矢量(5)式的Z域形式為

將(29),(30)式代入(32)式并整理得

將(33)式代入(31)式并整理得

根據(jù)Z域和離散時(shí)域?qū)?yīng)關(guān)系和移位定理[22],可以將以上Z域關(guān)系式過(guò)渡到離散時(shí)域,于是得到(33),(34)式對(duì)應(yīng)的離散時(shí)域形式為

綜上所述,可以得到飽和磁化鐵氧體的ZTFDTD計(jì)算步驟：

1)由B,M →H,用(36)式計(jì)算;

2)由H→M,用(35)式計(jì)算;

3)由H→E→B,用旋度方程(1),(2)式的差分離散式計(jì)算;

4)回到步驟1).

6 算例驗(yàn)證及數(shù)值結(jié)果

飽和磁化鐵氧體球半徑為d=1.5 cm,參數(shù)為

設(shè)x方向極化的平面高斯脈沖沿z方向入射,其中τ=34Δt和t0=0.8τ.

圖2(a)和(b)分別為該飽和磁化鐵氧體球的同極化和交叉極化后向RCS.圖中實(shí)線表示當(dāng)外加磁場(chǎng)沿z軸正向(θt=0,φt=0)時(shí),本文方法的計(jì)算結(jié)果,作為對(duì)比,圖中同時(shí)給出了文獻(xiàn)[21]的計(jì)算結(jié)果,在圖中用圓圈表示.由圖可見,兩者符合得非常好.

在圖2中同時(shí)給出了當(dāng)外加磁場(chǎng)沿x軸正向(θt=90°,φt=90°)時(shí)該球的后向 RCS,在圖 2 中用虛線表示.由圖2(a)可見,兩種外磁場(chǎng)方向情形下,該鐵氧體球的同極化后向RCS在頻率小于12GHz時(shí)幾乎無(wú)差別,在更高頻率時(shí)將略有不同;而又由圖2(b)所示,該鐵氧體球的交叉極化后向RCS對(duì)外磁場(chǎng)方向非常敏感,在該算例的兩種外磁場(chǎng)方向情形下,其差別很大.

7 結(jié)論

利用坐標(biāo)系轉(zhuǎn)換矩陣給出實(shí)驗(yàn)室系中飽和磁化鐵氧體的頻域磁化率張量,并采用部分分式展開方法通過(guò)Z變換得到磁化率張量的Z域形式,以Z變換為橋梁將飽和磁化鐵氧體的頻域磁本構(gòu)關(guān)系轉(zhuǎn)換到離散時(shí)域,實(shí)現(xiàn)了電磁場(chǎng)量B→H的迭代計(jì)算.從而給出了任意磁化方向鐵氧體電磁散射的ZT-FDTD算法,該方法對(duì)于研究鐵氧體目標(biāo)處于任意方向外磁場(chǎng)情況下的電磁散射問題具有應(yīng)用價(jià)值.作為驗(yàn)證,計(jì)算了飽和磁化鐵氧體球的后向RCS,所得結(jié)果與文獻(xiàn)結(jié)果一致.理論推導(dǎo)及算例表明該方法正確可行,推導(dǎo)簡(jiǎn)單,概念簡(jiǎn)明.

[1]He PB,Liu WM 2005 Phys.Rev.B 72 064410

[2]Li Z D,Liang JQ,Li L,Liu WM 2004 Phys.Rev.E 69 066611

[3]Luebbers R J,Hunsberger F,Kunz K S 1991 IEEE Trans.Antennas Propag.39 29

[4]Kelley DF,Lubbers RJ1996 IEEETrans.Antenn.Propag.44 792

[5]Chen Q,Katsurai M,Aoyagi PH 1998 IEEE Trans.Antenn.Propag.46 1739

[6]Liu SB,Mo JJ,Yuan N C 2004 Acta Phys.Sin.53 778(in Chinese)[劉少斌,莫錦軍,袁乃昌2004物理學(xué)報(bào)53 778]

[7]Xu L J,Yuan N C 2005 IEEEMicrow.Wirel.Compon.Lett.15 277

[8]Nickisch L J,Franke PM 1992 IEEEAntennas Propagat.Mag.34 33

[9]Takayama Y,Klaus W 1994 IEEEMicrow.Wirel.Compon.Lett.12 102

[10]Sullivan D M 1992 IEEETrans.Antenn.Propag.40 1223

[11]Sullivan D M 1995 IEEETrans.Antenn.Propag.43 676

[12]Sullivan D M 1996 IEEETrans.Antenn.Propag.44 28

[13]Ge D B,Wu Y L,Zhu X Q 2003 Chin.J.Radio Sci.18 359(in Chinese)[葛德彪,吳躍麗,朱湘琴2003電波科學(xué)學(xué)報(bào)18 359]

[14]Zhang Y Q,Ge D B 2009 Acta Phys.Sin.58 4573(in Chinese)[張玉強(qiáng),葛德彪2009物理學(xué)報(bào)58 4573]

[15]Chen R S,Edward K N Yung,Xie Z M,Wang Y 1998 International Journal of Infrared and Millimeter Waves19 329

[16]Edward K N Yung,Chen RS,Wu K 1998 IEEProc.H 145 433

[17]Yang L X,Ge D B 2006 Acta Phys.Sin.55 1751(in Chinese)[楊利霞,葛德彪2006物理學(xué)報(bào)55 1751]

[18]Yang L X,Ge D B,Wang G,Yan S 2007 Acta Phys.Sin.56 6937(in Chinese)[楊利霞,葛德彪,王剛,閻述2007物理學(xué)報(bào)56 6937]

[19]Yang L X,Ge D B,Zhao Y H,Wang G,Yan S 2008 Acta Phys.Sin.57 2936(in Chinese)[楊利霞,葛德彪,趙躍華,王剛,閻述2008物理學(xué)報(bào)57 2936]

[20]Ludeman L C 1986 Fundamentals of Digital Signal Processing(New York：Harper and Row Publishers)

[21]Schuster J,Lubbers R 1996 IEEE Antenn.Propag.Soc.AP-SInt.Symp.(Digest)3 1648

[22]Ge D B,Yan Y B 2011 Finite-Difference Time-Domain Method for Electromagnetic Waves(3th Ed.)(Xi’an：Xidian University Press)p11,272(in Chinese)[葛德彪,閆玉波2011電磁波時(shí)域有限差分法(第三版)(西安：西安電子科技大學(xué)出版社)第11,272頁(yè)]

[23]Bi DX 1985Electromagnetic Field Theory(Beijing：Publishing House of Electronics Industry)(in Chinese)[畢德顯1985電磁場(chǎng)理論(北京：電子工業(yè)出版社)]

[24]Kong JA 2002 Electromagnetic Wave Theory EMWPublishing(Beijing：Higher Education Press)