二维各向异性磁等离子体的无条件稳定ADE-CNAD-FDTD算法

详细信息    查看全文 | 推荐本文 |

英文篇名：Unconditionally stable auxiliary differential equation Crank-Nicolson-approximate-decoupling FDTD algorithm for 2-D anisotropic magnetized plasma 作者：李建雄 ; 庄永佳 ; 李现国 英文作者：Li Jianxiong;Zhuang Yongjia;Li Xianguo;School of Electronics and Information Engineering,Tianjin Polytechnic University;Tianjin Key Laboratory of Optoelectronic Detection Technology and Systems; 关键词：辅助微分方程 ; Crank-Nicolson ; approximate-decoupling算法 ; 时域有限差分 ; 磁等离子体 英文关键词：auxiliary differential equation(ADE);;Crank-Nicolson-approximate-decoupling(CNAD);;finite-difference time-domain(FDTD);;magnetized plasma 中文刊名：QJGY 英文刊名：High Power Laser and Particle Beams 机构：天津工业大学电子与信息工程学院;天津市光电检测技术和系统重点实验室; 出版日期：2017-12-22 10:58 出版单位：强激光与粒子束 年：2018 期：v.30;No.247 基金：国家自然科学基金项目(61372011) 语种：中文; 页：QJGY201801006 页数：5 CN：01 ISSN：51-1311/O4 分类号：22-26

摘要

针对二维各向异性磁等离子体提出一种有效的无条件稳定算法,新算法结合了辅助微分方程(ADE)方法与Crank-Nicolson approximate-decoupling(CNAD)时域有限差分算法仿真各向异性磁等离子体介质。传统的ADE-FDTD方法应用在一维各向异性色散介质具有较高的精度和效率,将提出的新算法ADECNAD-FDTD应用到二维各向异性磁等离子体介质中不仅解决了电磁波在具有各向异性和频率色散特性介质中传播的仿真难题,而且去除了CFL稳定性条件。该算法在保留了原有的精度情况下大幅度地提高了计算效率并成为无条件稳定的形式。给出一个算例证明该算法的有效性,通过模拟电磁波在磁等离子体中的传播,仿真结果与传统的ADE-FDTD算法对比,证实了该算法的高效率、无条件稳定性和高精度。
        An effective unconditionally stable implementation of the auxiliary differential equation Crank-Nicolson-approximate-decoupling finite-difference time-domain(ADE-CNAD-FDTD)algorithm for 2-D anisotropic magnetized plasma is proposed.The conventional ADE-FDTD method for 1-D anisotropic dispersive media has high efficiency and accuracy.This paper extends this method to 2-D anisotropic magnetized plasma with the CNAD scheme.The proposed formulations not only solves the problem that incorporates both anisotropy and frequency dispersion at the same time,but also eliminates the Courant-Friedrich-Levy(CFL)stability constraint.A numerical example has been carried out to validate the proposed formulations in the 2-D FDTD domain composed of anisotropic magnetized plasma.The results prove that the proposed formulations significantly save time and perform stably with acceptable accuracy.

引文

[1]Luebbers R J,Hunsberger F,Kunz K S.A frequency dependent time domain formulation for transient propagation in plasma[J].IEEE Trans Antennas and Propagation,1991,39(1):29-34.
    [2]刘少斌,莫锦军,袁乃昌.各向异性磁等离子体的辅助方程FDTD方法[J].物理学报,2004,53(7):2233-2236.(Liu Shaobin,Mo Jinjun,Yuan Naichang.An auxiliary differential equation FDTD method for anisotropic magnetized plasma.Acta Physica Sinica,2004,53(7):2233-2236)
    [3]Xu L J,Yuan N C.JEC-FDTD for 2-D conducting cylinder coated by anisotropic magnetized plasma[J].IEEE Microwave and Wireless Components Letters,2005,15(12):892-894.
    [4]Liu S,Zhong S,Liu S B.Piecewise linear recursive convolution FDTD method for magnetized plasmas[J].Journal of Systems Engineering and Electronics,2006,17(2):290-295.
    [5]杨宏伟,袁洪,陈如山,等.各向异性磁化等离子体的SO-FDTD方法[J].物理学报,2007,56(3):1443-1446.(Yang Hongwei,Yuan Hong,Chen Rushan,et al.SO-FDTD analysis of anisotropic magnetized plasma.Acta Physica Sinica,2007,56(3):1443-1446)
    [6]Taflove A,Hagness S C.Computational electrodynamics:The finite-difference time-domain method[M].3rd ed.Boston:Artech House,2005.
    [7]Namiki T.A new FDTD algorithm based on alternating-direction implicit method[J].IEEE Trans Microwave Theory and Techniques,1999,47(10):2003-2007.
    [8]Sun G,Trueman C W.Unconditionally stable Crank-Nicolson scheme for solving two-dimensional Maxwell's equations[J].Electronics Letters,2003,39(7):595-597.
    [9]Lee J,Fornberg B.A split step approach for the 3-D Maxwell’s equations[J].Journal of Computational and Applied Mathematics,2003,158(2):485-505.
    [10]Shibayama J,Muraki M.Efficient implicit FDTD algorithm based on locally one-dimensional scheme[J].Electronics Letters,2005,41(19):1046-1047.
    [11]Sun G,Trueman C W.Approximate Crank-Nicolson schemes for the 2-D finite-difference time-domain method for TEz waves[J].IEEE Trans Antennas and Propagation,2004,52(11):2963-2972.