径向基函数及其耦合方法在电磁场数值计算中的应用.pdfVIP

笛．：靼苗．相 计丌 算畀 切物 理理 第2警磊2期年月3 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Mar．．2(伪 [文章编号】1001．246x(2009)02．0299．12 径向基函数及其耦合方法在电磁场数值计算中的应用 张淮清， 俞集辉， 郑亚利 (重庆大学输配电装备及系统安全与新技术国家重点实验室，重庆400044) [摘 要】 针对带电导体附近急剧变化的位函数和场函数这一难于处理的边界条件，将小波函数的紧支撑特性 和全域径向基函数(RBF)的高精度逼近能力相结合，提出电磁场边值问题求解的耦合方法并应用于接地金属槽 ，箱的数值计算中；将径向基函数无网格方法引入波导本征值的计算中，给出其求解本征问题的思路，建立相应 的离散方程，分析矩形、圆形和脊形波导的本征值并与有限元方法进行比较．数值仿真实验表明，径向基函数 及其耦合方法在分析电磁场边值和本征值问题时是有效的且具有实现简单、节点少和精度高的优势． [关键词】 径向基函数；无网格法；拟小波；波导；本征值 [中图分类号】TMl51 [文献标识码】 A 0 引言 近年来，无网格法作为求解偏微分方程的新方法在物理和力学领域掀起了研究热潮，在电磁场计算中， 也有大量的文献予以报道．由于无网格法在求解中仅需节点信息，部分或彻底摆脱网格的限制，避免了传统 方法在网格剖分中的困难，而对近似方案的改进，解决了网格法因局部线性近似造成的场函数连续性变差的 问题．目前，无网格法按计算模型的不同主要分为配点型和Galerkin型…，因后者在求解中需依赖背景网格 进行积分近似计算，而配点式却仅与节点相关，所以是完全的无网格方法．以径向基函数作为插值函数的径 Basis 向基函数法是一类配点型无网格法，径向基函数(Radial 散布数据的插值，随着Kansa成功将其应用于抛物、双曲和椭圆方程求解，RBF方法逐渐成为求解偏微分方 程的重要方法． 边界条件的处理一直是无网格法研究中的关键问题，在Galerkin形式中，由于采用拟合而非插值近似， 边界条件主要通过拉格朗13乘子法、罚函数法、变换法及修正变分等方法来处理，极大地增加了计算的复杂 程度．在配点型方法中，虽然近似函数具有插值特性，但对一些强制边界条件，通过增加配点数的方式仍有很 大局限性，而且配点数越多就越可能导致离散方程奇异和近似解的振荡．通常，全域型RBF在逼近光滑、连 续和缓变函数时，具有很高的精度，但对强制边界下的突变函数，在边界附近存在较大误差．而在电气工程领 域，导体和介质问电性差异悬殊，由导体表面形成的边界附近，位函数和场函数变化剧烈，如在长接地金属槽 和输电线路模型中，带电导体附近函数变化率远大于其它区域，增加配点的改进能力十分有限．Kansa在进行 接地金属槽模型的研究时，提出改变RBF形状参数的方法，取得了一定效果，但提供的结论仅适用于特定模 型，不具普适性．然而，类似的强制边界条件却广泛存在于工程电磁场中，因此，寻求更好的解决方法具有重 要的意义．考虑到小波函数的紧支撑特性以及它对突变函数良好的逼近能力，将其引入到边界的处理中，用 RBF逼近域内缓变的函数，用小波逼近边界处突变函数，结合各自特点，发挥耦合方法的优势．在配点型的小 后者无解析表达式． 在电磁工程和微波技术领域，求解波导本征值问题，对于场型分析、传输特性研究和工程设计具有重要 意义．一般，对截面为矩形、圆形、椭圆等结构波导，可求得本征值的准确解或半解析解；对脊形、三角形等较 复杂截面结构，目前主要借助差分法、有限元、边界元等基于网格的方法求得数值解．网格法基本思想是网格 剖分和局部线性近似，它通过对波导截面的剖分形成单元，后在单元上采用线性近似，将微分方程化为代数 [收稿13期】2007—09—27；[修回日期】2008—02—24 [作者简介】张淮清(1979一)，男，安微霍山，讲师，从事电磁场数值计算、数字信号处理工作，重庆大学电气工程学院400044． 万方数据 计 算 物 理 第26卷 方程的求解，其精度取决于剖分质量和近似方案的选择．由于受限于剖分和局