三维柱坐标系中PML—FDTD算法及其奇点处理方法

总结

在三维圆柱坐标系下给出了FDTD有限差分格式并编程实现，并给出了对于圆柱坐标系带来的奇点问题解决方法，推导出对于解决磁场奇点问题的公式。利用辅助方程法重新推导出了PML中差分公式，通过编写程序，仿真测试。结果表明了本文处理奇点的可行性，正确性。

参考文献

[1] 葛德彪. 电磁波时域有限差分方法[M].西安：西安电子科技大学出版社，2011

[2] YEE K S. Numerical solution of initial bounary value problem involving Maxwell’s equations in isotropic media [ J ]. IEEE Trans Antennas and Propagation，1966，14：302 - 307.

[3] BERENGER J P. A perfectly matched layer for the absorption of electro-magnetic vaves [ J ]. Computational Physics，1994，114（2）：185 - 200.

[4] ZHAO L，CANGELLAR IS A C. GT - PML：Generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-do-main grids [ J ]. IEEE Trans Microwave Theory and Techniques，1996，44（12）：2555 - 2563.

[5] GHEONJ IAN A，JOBAVA R. Non2uniform conforming mesh generator for FDTD scheme in 3D cylindrical coor-dinate system [A]. IEEE D IPED Proceedings[C]. [ s.l. ]：[ s. n. ]，2000：41 - 44

[6] Y.Chen，R.Mittra，and P.Harms，“Finite-difference time-domain algorithm for solving Maxwell’s equations in rotationally symmetric geometries，” IEEE Trans. Microwave Theory Tech.，vol.44，pp.832-839，June 1996.

[7] D.W.Prather and S.Shi，“Formulation and application of the finite-difference time-domain method for the analysis of axially symmetric diffractive optical elements，” J.Opt. Soc.Amer.A.，vol. 16，pp.1131-1141，1999.

[8] F. Liu，S. Crozier，“An FDTD model for calculation of gradient-induced eddy currents in MRI system，” IEEE Trans Appl Supercond.，vol. 14，pp. 1983-1989，2004.

[9] Zilvinas Kancleris，“Handling of Singularity in Finite-Difference Time-Domain Procedure for Solving Maxwell’s Equations in Cylindrical Coordinate System，” IEEE Trans. Antennas Propagat.，vol. 56，pp. 610-613，2008.

[10] Jieru Chi，Feng Liu，Ling Xia，Tingting Shao，David G. Mason，“An Improved Cylindrical FDTD Method and its Application to Field-tissue Interaction Study in MRI，”IEEE trans on antennas and propagation，vol. 47，no. 2，rebruary 2011

推荐访问:奇点 坐标系 算法 方法 FDTD