DGTD算法在電磁問題中的應用與分析.pdf

1、現(xiàn)在電磁工程的計算大都采用傳統(tǒng)的有限元法(Finite Element Method, FEM)或時域有限差分法(Finite Difference Time Domain, FDTD)等方法，然而這些方法在求解大規(guī)模多尺度的電磁問題上還面臨著很多挑戰(zhàn)。由于越來越多的電磁問題趨于大規(guī)模復雜化，因此選擇一種高效的算法具有很重要的意義。本文研究的時域間斷伽遼金方法(Discontinuous Galerkin Time Domain, DG

2、TD)，它是結(jié)合FEM和有限體積法(Finite Volume Method, FVM)優(yōu)點的算法，具有非結(jié)構(gòu)化網(wǎng)格離散優(yōu)勢，具有高階基函數(shù)特性的高計算精度，具備內(nèi)在的并行性，更適用于復雜多尺度的電磁問題計算。本文研究DGTD算法在電磁問題中的應用與分析，主要是針對一維、二維時域Maxwell方程的DGTD法，討論分析不同階數(shù)基函數(shù)下的穩(wěn)定性，以及不同網(wǎng)格密度和不同階數(shù)對結(jié)果精度的影響。
　　本文首先介紹了傳統(tǒng)的FEM以及時域有限

3、體積法(Finite Volume Time Domain, FVTD)的基本原理以及關鍵技術(shù)，包括 Galerkin法和數(shù)值通量的概念，以此為基礎引出DGTD法的原理及其關鍵技術(shù)。然后討論了DGTD算法的基本原理及關鍵技術(shù)，包括半離散公式的導出，質(zhì)量矩陣和剛度矩陣的計算，基函數(shù)的選擇，單元內(nèi)離散點的選擇，數(shù)值通量的概念以及時間偏導數(shù)的處理。接著討論一維、二維時域Maxwell方程的DGTD法及其關鍵技術(shù)。DGTD算法采用Legendr

4、e多項式作為單元高階基函數(shù)，并采用Legendre-Gauss-Lobatto離散點作為單元內(nèi)積分點，時間偏導的處理采用了4階低存儲Runge-Kutta時間方案。
　　最后對一維、二維時域電磁問題進行DGTD計算和分析。先通過Matlab語言編程實現(xiàn)一維算例，驗證了 DGTD法的可行性以及程序的正確性。然后又利用Gmsh軟件繪圖繪制出二維矩形波導和圓形波導并剖分網(wǎng)格，通過程序讀取網(wǎng)格數(shù)據(jù)到DGTD算法模塊；計算出波導截面的電磁場