论文部分内容阅读
无网格法是在对一个问题域建立离散的系统方程时不用事先定义网格的一种数值方法,它是利用一组散布在问题域中以及域边界上的节点表示(而非离散)该问题域和其边界,并采用完全基于点的近似,彻底或部分地消除了网格,完全避免了网格的再生成。径向点插值法(RPIM)采用耦合径向基函数和多项式基函数构造形函数,不仅有效的解决了点插值法(PIM)系数矩阵的奇异性问题,而且插值函数具有δ函数性质,从而克服了以往无网格法难以施加位移边界条件的难点。 近几年,无网格径向点插值法(RPIM)作为一种新型无网格法,被越来越多的用于求解电磁场模型问题所带来的类似于网格依赖性强、精度低、计算时间长、边界条件难于施加等问题,并且取得了较好的结果。 首先,本文在深入研究时域有限差分法理论的基础上,研究了电磁斗篷材料的理论方法,通过引入Drude方程推导出RPIM方法对电磁斗篷材料的实现方法; 其次,以电磁斗篷为物理模型,详细介绍了无网格(RPIM)方法原理,给出了形函数构造方法,并结合麦克斯韦方程组结合Drude离散系统方程理论推导出电磁斗篷材料的场分量迭代方程,以及完美匹配层(PML)的RPIM实现方法; 第三,论述了程序的设计思想,绘制出程序结构流程图,在此基础上使用Matlab编写了基于径向点插值法程序,通过计算,对电磁斗篷算例中的不同参数(包括形参、节点分布等)对精度的影响进行了实验验证; 最后,通过设置不同时步长和斗篷半径大小,验证了电磁波在斗篷中的传播情况,以及该方法的可行性;通过与时域有限差分法数值比较发现,该方法具有更高的精度和更快的计算速度。