第四章 结 论
通过对时域差分法的学习而进行的关于二维FDTD的TE波仿真，巩固了迭代式的形式与相应的循环编程技巧。但也存在一定问题，由于在吸收上我的编程思路不是很规范，所以在吸收效果上与理论上的完全匹配层吸收效果还有差距，但是已经基本达到了我的预期目的。可能是我的分区方法不对的缘故，原因就是出在sigma_max的设置上面，在高斯制下长度单位变为cm，所以同样的空间步长在高斯制下数量级就扩大了100倍，而sigma_max的计算式里面分母含有空间步长，所以如果不加以修正，sigma_max在高斯制下就变成了国际制下的1%，太大的导电率和太小的导电率都不能实现较好的吸收效果，下阶段应将其做适当修正并且写粒子仿真时用高斯制Maxwell方程可以参考角点反射大，多设置几层PML试一试。
参考文献
[1]J.M.Jin, M.Zunoubi, K.C.Donepudi, W.C.Chew, Frequency-Domain and Time-Domain Finite-Element Solution of Maxwell’s Equations Using Spectral Lanczos Decomposition Method. Comput. Methods Appl. Mech. Engrg. 1999,169:279~296
[2]K.Dongsoo,H.B.Lee. Hybird Full-Wave Analysis of Via-hole Grounds Using Finit-Difference and Finite-Element Time-Domain Methods. IEEE Transactions on Microwave Theory and Techniques. 1997,45(12):2217~2222 [3] Kane S. Yee Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media[J].IEEE Trans Antennas Propagation, 1966,AP-14(3):302~307 [4] A.Taylor and M.E.Brodwin. Numerical solution of steady-state EM scattering problems using the time-Maxwell’s equations. [J] IEEE Trans.Microwave Theroy Tech.,Aug. 1975，MTT-23:623~630
[5] X.Min,W.Sun and K.M.Chen. Stability analysis of the finite difference time domain method applied to unbounded electromagnetic problem [J]. IEEE Antennas and Propagation Society International Symposium,Dallas:May 1990,(4):1640~1643
[6]A.Taflove. Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures [J]. Wave Motion,June 1988,10(6): 547~582
[7] I.S.Kim and W.J.R.Hoefer. Numerical dispersion characteristics and stability factor for the TD-FD method [J].Electronics Letters, July 1990, 26(7): 485~487
[8] K.L.shlager,J.G.Msloney, S.L.Ray,and etc.. Relative accuracy of several finite-difference time-domain methods in two and three dimensions [J]. IEEE Transactions on Antennas and Propagation,Dec.1993,AP-41(12)：1732~1737
[9] S.L.Ray. Numerical dispersion and stability characteristics of finite-difference time-domain methods on nonorthogonal meshes [J]. IEEE Transactions on Antennas and Propagation,Feb.1993,AP-41(2)：233~235 [10] H.Shi and J.L.Drewriak. Dispersive comparison for DSI- and tenor-based nonorthogonal FDTD [J]. IEEE Microwave and Guided Wave Letters, May 1996,6(5):193~195
[11] F.Xiao and H.Yabe. Numerical dispersion relation for FDTD method in general curvilinear coordinates [J]. IEEE Microwave and Guided Wave Letters, Feb.1997,7(2):48~50
[12] R.Holland and J.Williams．Total field versus scattered field finite difference codes:a comparative assessment[J].IEEE Transactions on Nuclear Science，Dec.1983，NS-30(6)：4583～4588
[13]G.Mur．Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations[J]．IEEE Transactions on Electromagnetic Compatibility，Nov.1981，EMC-23(4)：377～382
[14]B.Engquist and A.Majda. Absorbing boundary conditions for the numerical simulation of waves [J].Mathematics of the Computation，July 1977，31(139)：629～651
[15]J.P.Berenger．A perfectly matched layer for the absorption of electromagnetic waves[J]．Journal of Computational Physics，Oct.1994，114(2)：185～200
[16]J.P.Berenger．Perfectly matched layer for the FDTD solution of wave-structure interaction problem[J]．IEEE Transactions on Antennas Propagation，Jan.1996，AP-44(1)：110～117
[17]J.P.Berenger．Three-dimensional perfectly matched layer for the absorption of electromagnetic waves[J]．Journal of Computational Physics，Sept.1996，127(2)：363～379
[18]Z.S.Sacks，D.M.Kingsland，R.Lee and J.F.Lee．A perfectly matched anisotropic absorber for use as an absorbing boundary condition[J]．IEEE
Transactions on Antennas and Propagation，Dec.1995，AP-43(12)：1460～1463
[19]D.E.Merewether.Transient currents on a body of revolution by an electromagnetic pulse [J]. IEEE Transactions on Electromagnetic Compatibility，May 1971，EMC-13(2)：41～44
[20]R.Holland．THREDS: a finite-difference time-domain EMP code in 3D spherical coordinates [J]．IEEE Transactions on Nuclear Science，Dec.1983，NS-30(6)：4592～4595
[21]S.S.Zivanovic，K.S.Yee and K.K.Mei．A subgridding method for the time-domain finite-difference method to solve Maxwell’s equations [J]．IEEE Transactions on Microwave Theory and Techniques，Mar.1991，MTT-39(3)：471～479
[22]P.Thoma and T.Weiland．A consistent subgridding scheme for the finite difference time domain method[J]．International Journal of Numerical Modeling: Electronic Network Device and Fields，1996，9(5)：359～374 [23]K.S.Yee，D.Ingham and K.Shlager．Time-domain extrapolation to the far field based on FDTD calculations[J]．IEEE Transactions on Antennas and Propagation，Mar.1991，AP-39(3)：410～413
[24]R.J.Luebbers，K.S.Kunz，M.Schneider and etc.．A finite-difference time-domain near zone to far zone transformation[J]．IEEE Transactions on Antennas and Propagation，Apr.1991，AP-39(4)：429～433
[25]R.J.Luebbers，D.Ryan and J.Beggs．A two dimensional time-domain near-zone to far-zone transformation[J]．IEEE Transactions on Antennas and Propagation，July 1992，AP-40(7)：848～851
[26]M.Kragalalott，M.S.Kluskens and W.P.Pala．Time-domain fields exterior to a two dimensional FDTD space[J]．IEEE Transactions on Antennas and Propagation，Nov.1997，AP-45(11)：1655～1663
[27]冯慈璋，马西奎. 工程电磁场导论 [M]. 高等教育出版社 2003 [28] D.M.Sullivan．Electromagnetic Simulation Using the FDTD Method[M]．New York：IEEE Press，2000
[29]B.J.Zook. A Dispersive Outer Radiation Boundary Condition for FDTD Calculation. Proceedings of the 10th Annual Review of Progress in Applied Computational Electromagnetics. 1994,1:240～247
[30]Engquist B, Majda A. Absorbing Boundary Conditions for the Numerical Simulation of Waves. Mathematics of Computation. 1977,31(7):629～651 [31]T.G.Moore,J.G.Blaschak,A.Taflobe,G.A.Kriegsmann. Theory and Application of Radiation Boundary Operators.IEEE Trans. AP.1988,36(12)1797～1812