首页|柱坐标系"分解"思想WLP-FDTD的PML实现

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

扫码查看
从柱坐标系下传统时域有限差分(FDTD)的基本方程出发,解决传统加权 Laguerre 多项式(WLP)-FDTD算法在计算上存在内存消耗大和计算效率低的问题.本文的"分解"思想分为两部分,首先,在频域上对电磁场方程进行第一次分解,并代入 PML参数,将分解后的时域方程转换至 Laguerre域,使原三维双向求解问题转换为二维单向规模来求解,降低计算的内存消耗.其次,通过 LU分解对求解规模降低后的 Laguerre域系数矩阵实施第二次分解,实现规避大型稀疏矩阵求解的第一步,进而通过追赶法求解三对角电磁场方程,以提高计算效率,带动内存消耗的降低.算例证明:对比传统 WLP-FDTD法,本文算法可在内存消耗上降低 57%,计算效率提高 49%左右,且在不丢失精度的情况下,具有较好的电磁波吸收效果,误差反射系数可达-70 dB.
PML Implementation of WLP-FDTD Algorithm Based on"Decomposition"Ideology in Cylindrical Coordinate System
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.

weighted-Laguerre-polynomial-FDTDdecomposition ideologyLU decompositionchase after methodreduce the scale of solution

朱大伟、陈海林、徐佳琛、周晓晓、徐博奥

展开 >

南京工业职业技术大学电气工程学院,江苏 南京 210023

陆军工程大学电磁环境效应与光电工程国家级重点实验室,江苏 南京 210007

南京理工大学电子工程与光电技术学院,江苏 南京 210094

中国人民解放军 96901 部队,北京 100000

展开 >

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

国家自然科学基金南京工业职业技术大学引进人才科研启动基金

61701221YK21-02-01

2024

计算物理
中国核学会

计算物理

CSTPCD北大核心
影响因子:0.366
ISSN:1001-246X
年,卷(期):2024.41(4)