時域有限差分法(FDTD算法)的基本原理及仿真

1、精選文檔時域有限差分法（FDTD算法）時域有限差分法是1966年K.S.Yee發(fā)表在AP上的一篇論文建立起來的，后被稱為Yee網(wǎng)格空間離散方式。這種方法通過將Maxwell旋度方程轉(zhuǎn)化為有限差分式而直接在時域求解, 通過建立時間離散的遞進序列, 在相互交織的網(wǎng)格空間中交替計算電場和磁場。FDTD算法的基本思想是把帶時間變量的Maxwell旋度方程轉(zhuǎn)化為差分形式，模擬出電子脈沖和抱負(fù)導(dǎo)體作用的時域響應(yīng)。需要考慮的三點是差分格式、解的穩(wěn)定性、吸取邊界條件。有限差分通常接受的步驟是：接受肯定的網(wǎng)格劃分方式離散化場域；對場內(nèi)的偏微分方程及各種邊界條件進行差分別散化處理，建立差分格式，得到差分方程組；結(jié)

2、合選定的代數(shù)方程組的解法，編制程序，求邊值問題的數(shù)值解。1. FDTD的基本原理FDTD方法由Maxwell旋度方程的微分形式動身，利用二階精度的中心差分近似，直接將微分運算轉(zhuǎn)換為差分運算，這樣達到了在肯定體積內(nèi)和一段時間上對連續(xù)電磁場數(shù)據(jù)的抽樣壓縮。Maxwell方程的旋度方程組為： （1）在直角坐標(biāo)系中，（1）式可化為如下六個標(biāo)量方程：， （2）上面的六個偏微分方程是FDTD算法的基礎(chǔ)。Yee首先在空間上建立矩形差分網(wǎng)格，在時刻時刻，F(xiàn)(x,y,z)可以寫成 （3）用中心差分取二階精度：對空間離散：對時間離散： （4）Yee把空間任一網(wǎng)格上的E和H的六個重量，如下圖放置:圖1 Yee氏網(wǎng)格

3、及其電磁場重量分布在FDTD中，空間上連續(xù)分布的電磁場物理量離散的空間排布如圖所示。由圖可見，電場和磁場重量在空間交叉放置，各重量的空間相對位置也適合于Maxwell方程的差分計算，能夠恰當(dāng)?shù)孛枋鲭姶艌龅膫鞑ヌ匦?。同時，電場和磁場在時間上交替抽樣，抽樣時間間隔相差半個時間步，使Maxwell旋度方程離散以后構(gòu)成顯式差分方程，從而可以在時間上迭代求解，而不需要進行矩陣求逆運算。因此，由給定相應(yīng)電磁問題的初始條件，F(xiàn)DTD就可以逐步推動地求得以后各個時刻空間電磁場的分布。依據(jù)這一原則可以寫出六個差分方程:（5）其余的也如法可以寫出，每個網(wǎng)格點上的個場分兩的新值依靠于該點在前一時間步長時刻的值機該點

4、四周的接近點上另一場量在早半個時間步長時的值。因此任一時刻可一次算出一個點，并行算法可計算出多個點。通過這些運算可以交替算出電場磁場在各個時間步的值。已知 0時刻空間各處的電磁場初始值依據(jù)上述FDTD差分方程組可得出計算電磁場的時域推動計算方法，如圖2所示。計算 時刻空間各處的磁場值 計算 時刻空間各處的電場值 循環(huán)n次圖2 FDTD在時域的交叉半步逐步推動計算2.數(shù)值穩(wěn)定性條件時間步長，空間步長，必需滿足肯定的關(guān)系，否則就使得數(shù)值表現(xiàn)不穩(wěn)定，表現(xiàn)為：隨著計算步數(shù)的增加，計算場量的數(shù)值會無限的增大，這種增大不是由于誤差積累造成的，而是由于電磁波的傳播關(guān)系被破壞造成的。所以，必需滿足肯定的關(guān)系以

5、保證穩(wěn)定性。Taflove等在1975年對Yee氏差分格式的穩(wěn)定性進行了爭辯，并導(dǎo)出了對時間步長的限制條件。數(shù)值解是否穩(wěn)定主要取決于時間步長與空間步長、的關(guān)系。對于非均勻媒質(zhì)構(gòu)成的計算空間選用如下的穩(wěn)定性條件： （6）若接受均勻立方體網(wǎng)格：， （7）而一般?。?，c為光速。當(dāng)，不相等時， （8）3.數(shù)值色散FDTD網(wǎng)格中，會導(dǎo)致數(shù)字波模在網(wǎng)格中發(fā)生轉(zhuǎn)變，這種轉(zhuǎn)變是由于計算網(wǎng)格本身引起的，而非物理因素，所以必需考慮。即在FDTD網(wǎng)格中，電磁波的相速與頻率有關(guān)，電磁波的相速度隨波長、傳播方向及變量離散化的狀況不同而轉(zhuǎn)變。色散將導(dǎo)致非物理因素引起的脈沖波形畸變、人為的各向異性和虛假折射等現(xiàn)象。明顯，色

6、散與空間、時間的離散間隔有關(guān)，如下式所示： （9）與數(shù)值色散關(guān)系相對應(yīng)，在無耗介質(zhì)中的單色平面波，色散解析關(guān)系是： （10）由式（9）可知，當(dāng)式（9）中的、均趨于零時，它就趨于式（10）。也就是說數(shù)值色散是由于用近似差分替代連續(xù)微分而引起的，而且在理論上可以減小到任意程度，只要此時時間步長和空間步長都足夠小。為獲得抱負(fù)的色散關(guān)系，問題空間分割應(yīng)依據(jù)小于正常網(wǎng)格的原則進行。一般選取的最大空間步長為，為所爭辯范圍內(nèi)電磁波的最小波長。由上分析說明，數(shù)值色散在用FDTD法分析電磁場傳播中的影響是不行能避開的，但我們可以盡可能的減小數(shù)值色散的影響?，F(xiàn)在適當(dāng)選取時間和空間步長，傳播方向，可以得到抱負(fù)狀況，如下所示：3-D方形網(wǎng)格：（數(shù)值穩(wěn)定的極限狀態(tài)，可得抱負(fù)色散關(guān)系）取波沿對角線傳播， （11）2-D方形網(wǎng)格:也是沿對角線傳播， （12）1-D網(wǎng)格： （13）4.吸取邊界條件在電磁場的輻射和散射問題中，邊界總是開放的，電磁場占據(jù)無限大空間，而計算機內(nèi)存是有限的，所以只能模擬有限空間。即：時域有限差分網(wǎng)格將在某處被截斷。這要求在網(wǎng)格截斷處不能引起波的明顯反射，因而對向外傳播的波而言，就像在無限大的空間傳播一樣，一種行之有效的方法是在截斷處設(shè)置一種吸取邊界條件。使傳播到截斷出的波被邊界吸取而不產(chǎn)生反射。下面只給出Engquist-Majda吸取邊界條件，接受Mur差分格式，其總體虛