柱坐标系“分解”思想WLP-FDTD的PML实现  

PML Implementation of WLP-FDTD Algorithm Based on“Decomposition”Ideology in Cylindrical Coordinate System

作　　者：朱大伟 陈海林 徐佳琛 周晓晓 徐博奥 ZHU Dawei;CHEN Hailin;XU Jiachen;ZHOU Xiaoxiao;XU Boao(School of Electrical Engineering,Nanjing Vocational University of Industry Technology,Nanjing,Jiangsu 210023,China;National Key Laboratory on Electromagnetic Environmental Effects and Electro-optical Engineering,PLA Army Engineering University,Nanjing,Jiangsu 210007,China;School of Electronic and Optical Engineering,Nanjing University of Science and Technology,Nanjing,Jiangsu 210094,China;Unit of PLA,Beijing 100000,China)

机构地区：[1]南京工业职业技术大学电气工程学院,江苏南京210023 [2]陆军工程大学电磁环境效应与光电工程国家级重点实验室,江苏南京210007 [3]南京理工大学电子工程与光电技术学院,江苏南京210094 [4]中国人民解放军96901部队,北京100000

出　　处：《计算物理》2024年第4期440-452,共13页Chinese Journal of Computational Physics

基　　金：国家自然科学基金(61701221);南京工业职业技术大学引进人才科研启动基金(YK21-02-01)资助项目。

摘　　要：从柱坐标系下传统时域有限差分(FDTD)的基本方程出发,解决传统加权Laguerre多项式(WLP)-FDTD算法在计算上存在内存消耗大和计算效率低的问题。本文的“分解”思想分为两部分,首先,在频域上对电磁场方程进行第一次分解,并代入PML参数,将分解后的时域方程转换至Laguerre域,使原三维双向求解问题转换为二维单向规模来求解,降低计算的内存消耗。其次,通过LU分解对求解规模降低后的Laguerre域系数矩阵实施第二次分解,实现规避大型稀疏矩阵求解的第一步,进而通过追赶法求解三对角电磁场方程,以提高计算效率,带动内存消耗的降低。算例证明:对比传统WLP-FDTD法,本文算法可在内存消耗上降低57%,计算效率提高49%左右,且在不丢失精度的情况下,具有较好的电磁波吸收效果,误差反射系数可达-70 dB。Based on the basic equations of conventional finite-difference-time-domian(FDTD)algorithm in cylindrical coordinate system,in order to solve the problems of large memory consumption and low computational efficiency in computation of conventional weighted-Laguerre-polynomial(WLP)-FDTD algorithm.The idea of“Decomposition”in this paper is divided into two parts.Firstly,the electromagnetic field equation is decomposed for the first time in the frequency domain,and PML parameters are substituted,and the decomposed time-domain equation is converted to the Laguerre domain,so that the original three-dimensional bidirectional solution problem is converted into two-dimensional one-way scale to solve,reducing the memory consumption of calculation.Secondly,LU decomposition is used to decompose the coefficient matrix of Laguerre domain after solution scale reduction,which realizes the first step of avoiding large sparse matrix,and then the chase after method is used to solve the conventional electromagnetic field equation,so as to improve the computational efficiency and reduce memory consumption.Numerical examples show that comparing with conventional WLP-FDTD,the proposed scheme can reduce the memory consumption by 57%and increase the computing efficiency by 49%,without losing the accuracy,the proposed method has a good electromagnetic wave absorption effect,and the reflection error can reach-70 dB.

关 键 词：加权Laguerre-FDTD法 分解思想 LU分解 追赶法 降低求解规模 

分 类 号：O441[理学—电磁学]