dgtd算法在电磁问题中的应用与研究word格式论文.docxVIP

摘要现在电磁工程的计算大都采用传统的有限元法(Finite Element Method, FEM) 或时域有限差分法(FiniteDifferenceTimeDomain,FDTD)等方法，然而这些方法在求解大规模多尺度的电磁问题上还面临着很多挑战。由于越来越多的电磁问题趋于大规模复杂化，因此选择一种高效的算法具有很重要的意义。本文研究的时域间断伽辽金方法(Discontinuous GalerkinTimeDomain, DGTD) ，它是结合FEM 和有限体积法(FiniteVolume Method, FVM)优点的算法，具有非结构化网格离散优势，具有高阶基函数特性的高计算精度，具备内在的并行性，更适用于复杂多尺度的电磁问题计算。本文研究DGTD算法在电磁问题中的应用与分析，主要是针对一维、二维时域Maxwell方程的DGTD法，讨论分析不同阶数基函数下的稳定性，以及不同网格密度和不同阶数对结果精度的影响。本文首先介绍了传统的FEM以及时域有限体积法(FiniteVolumeTime Domain, FVTD)的基本原理以及关键技术，包括Galerkin法和数值通量的概念，以此为基础引出DGTD法的原理及其关键技术。然后讨论了DGTD算法的基本原理及关键技术，包括半离散公式的导出，质量矩阵和刚度矩阵的计算，基函数的选择，单元内离散点的选择，数值通量的概念以及时间偏导数的处理。接着讨论一维、二维时域Maxwell 方程的DGTD 法及其关键技术。DGTD算法采用Legendre多项式作为单元高阶基函数，并采用Legendre-Gauss-Lobatto 离散点作为单元内积分点，时间偏导的处理采用了4 阶低存储Runge-Kutta 时间方案。最后对一维、二维时域电磁问题进行DGTD 计算和分析。先通过Matlab 语言编程实现一维算例，验证了DGTD法的可行性以及程序的正确性。然后又利用Gmsh 软件绘图绘制出二维矩形波导和圆形波导并剖分网格，通过程序读取网格数据到DGTD算法模块；计算出波导截面的电磁场，绘制出电磁场分布图、电磁场随时间变化的曲线图以及谐振频率。通过改变基函数的阶数，并绘制不同阶数下的电磁场随时间变化曲线图；研究了在相同网格密度条件下，选择不同的基函数阶数对电磁场随时间变化曲线产生的影响；研究了不同网格密度、不同基函数阶数等参数下的运行时间、计算结果的误差以及占用的存储空间，得到不同条件下的计算效率。关键词：DGTD，伽辽金法，勒让德多项式，4 阶低存储龙格库塔法，时域Maxwell方程ABSTRACTThetraditionalmethodsappliedtoelectromagneticproblems,suchasfinite elementmethod(FEM)andfinitedifferencetimedomain(FDTD),arefacingmany challengesforsolvinglarge-scaleelectromagneticproblems.Asmoreandmore electromagneticproblemstendtobelarge-scaleandcomplicated,itsofgreat significancetoseekoutahigherefficientalgorithm.Inthispaper,thediscontinuous Galerkintimedomainmethod (DGTD),whichiscombining theadvantagesof FEM and finitevolumemethod(FVM),issuitableforcalculatingcomplexmulti-scale electromagneticproblems.Ithastheadvantagesofunstructuredgrids,higher-order basisfunctionswithhighprecision andinherent parallelism.Inthisstudy,thecoreis the applicationandanalysisof DGTD algorithminelectromagnetic problemsfor one-dimensionalandtwo-dimensionalMaxwellequationsintime-domain,discussing andanalyzingthestabilitywithdifferentorderofthebasisfunctions,aswellasthe impact on the results accuracy by different mesh densit