基于等几何分析方法求解任意截面波导本征问题.pdfVIP

第29卷第4期 计 算 物 理 V01．29，No．4 OF PHYSICS 2012年7月 CHINESEJOURNALCOMPUTATIONAL July，2012 文章编号：1001-246X(2012)04-0534-09 基于等几何分析方法求解任意截面波导本征问题 张 勇， 林 皋， 胡志强， 刘 俊 (大连理工大学建设工程学部，辽宁大连116024) 摘要：基于等几何分析方法具有自由度花费少、高精度、高阶连续性等特点，通过加权余量法对椭圆波导本征问 题的亥姆霍兹方程等几何离散得出等几何分析方程．解决了传统数值方法的求解域与几何模型的非一致性问题， 实现了将问题的分析计算构架于精确几何模型基础之上．分析任意截面波导的本征问题，对不同偏心率的椭圆波 导以及三角形和五边形波导的截止波数的求解结果显示等几何分析方法求解波导本征问题的高效及高精度特性． 与传统方法相比，此方法以较少的自由度消耗便会达到较高的求解精度，并且数值解的收敛率较快． 关键词：波导本征问题；等几何分析；非均匀有理B样条；椭圆波导；三角形波导；五边形波导 中图分类号：TMt51；024；018文献标识码：A O 引言 等几何分析方法(Isogeometric 法解决了传统数值方法的求解域与几何模型非一致性问题，从而实现了将问题的分析计算构架于精确几何 模型基础之上，达到了较高的计算精度与效率．目前这种方法已被用于求解固体力学和流体力学问题旧。。． 其基本思想是直接利用实体的几何描述进行分析计算，避免传统有限元的二次建模．这样，可以消除几何模 型和计算模型之问存在的误差，同时也保留了几何模型中的几何光滑和连续的特点，提高了未知场的连续 阶，从有限元无法突破的co连续到口。1的高阶连续性．另外，得益于计算几何中丰富的细分策略，等几何分 析方法可以自动对计算模型进行保形细分，并且在细分过程中无需和原始模型进行通讯以及保持细分前后 模型的不变性． 在电磁工程和微波技术领域中广泛应用各类截面波导来传导电磁波，按照电磁波传播模式可分为：横向 中电磁波的传输模式和模态的选择与此问题有关． 波导本征值问题的求解主要有解析法、数值解法两大类，解析法包括以分离变量法、保角变换伸。为代表 的严格解析法，以多级理论、模态分析法为代表的半解析法和以微扰法、变分法和算子理论方法为代表的近 似解析法．其中，严格解析法只适用于简单的几何边界波导本征值问题．半解析方法也是只适用一些几何形 状比较规则的波导本征值问题．在近似解析法中，微扰法依赖于微扰量，当微扰量很大时，计算结果的误差较 大；变分法可以求解很多本征值问题，但求解效率和精度在很大程度上依赖于泛函的构造、试探函数的形式 和展开项数；算子理论在去掉伪模的影响方面具有一定的优势，但很难处理复杂形状边界条件的本征值 问题． 0|、有限元法¨“、边界元法。1“、无 求解电磁场本征值问题的数值方法很多，主要有经典的有限差分法一。1 网格法¨3‘H3等．在以上各种数值计算方法中，有限元方法是分析波导本征问题适应性最强、最通用的方法之 一．各类数值方法都需要对求解域或域边界进行离散，而离散之后的计算模型可能与原模型不具有一致性， 从而丧失精度． 本文应用等几何分析方法求解任意截面波导本征问题．以等参概念构造纵向电场强度(TM波)或者纵向磁 收稿日期：2011一09—02；修回日期：2012—02—10 万方数据 第4期 张 勇等：基于等几何分析方法求解任意截面波导本征问题 535 场强度(TE波)，并利用加权余量法对波导控制方程，即亥姆霍兹方程进行处理，推导出波导特征问题的等几何 分析方程．施加边界条件后，可求得相应TE波或TM波的截止波数．数值算例中，针对两种典型偏心率的椭圆 波导和以正三角形和正五边形波导为例的正多边形波导，将本文方法的截止波数数值解与解析解和其他数值 解相比较，结果表明本文方法求解椭圆波导是有效的，