FDTD的DQM改进方法.pdfVIP

FDTD的DQM改进方法 钱晓宁 毛军发 (上海交通大学电子丁．程系，f’海200030) 而减少时域有限差分法(FDTD)算法中的网络节点数目，节省计算机的内存，提高FDTD的 分析效率。 随着信号处理速度的不断提高及电路规模的不断扩大，高速集成电路的互连及封装结 构已越来越有必要用全波方法进行建模和分析。甚至已出现用全波方法对整个微电子系统 进行电磁场模拟的趋势。但是，传统的全波方法(如FDTD方法【1】)计算消耗大，提高全 波方法的计算效率至关重要。有两条途径可以提高全波方法的效率：一是改进吸收边界条 件以减小计算区域【2】【3】；二是对全波方法本身加以改进，如文献[41151中提出的基于小波 多分辨分析的MRTD方法和类似算法。 Quadrature 高其计算效率。DQM是一种宜接的数值计算方法，与差分和有限元法相比，它的计算量 可以大大降低。而且具有较高的精度。它的主要思想是将某一坐标方向上的偏微分算予用 在该方向一系列的离散点函数值加权逼近，也即数值积分。与传统的FDTD相比，DQM可 以用较少的剖分点，从而提高FDTD的计算效率。 由于传统FDTD本身吸收边界条件带来的数值反射，在此基础上若应用原始的DQM， 误差的积累会造成计算结果的不稳定。另一方面，电磁场在某一坐标方向上的分布可能有 很多零点，它们对偏微分算子的贡献也为零。因此，本文采用跟踪加窗的方法对DQM加 以改进．即规定在某～节点处场的偏微分值只使用其相邻的一定数目节点处的场值来加权 逼近，以保证计算的稳定性。并提高计算的效率。 一．原理 1．DQM 如果函数W(x，y,zt)在某个坐标方向上是足够光滑的，在这一方向上的偏微分就可以 用积分求导法来逼近[51，如在z方向： 型坐0型：兰咿b币，)z l=-J’，…。J7 其中z．为Z疗向分段点的坐标，系数a”的if簿可以由多种疗法确定，～艘忠譬选定刖錾 函数，使其对式(1)精确成立，从而得到相应的组全系数。常Hja々堪函数订覃，I：j 股的幂函数，将一细幂函数代入式r1)得到相应线性方程绢： 一 唧 才 ll拓 k=O，l，1、，N，i=1,2，·，N (2) 。∑川 解之得到系数～。这里，如果采用正交函数集结果将得到改善。 2．FDTD与DQM的结台 ～133一 在FD3D计算中．只要其场值分布足够光滑，司阻在一个或多个方向用DQM计算其 偏微分的数值近似，为简单起见，我们仅在z方向使用DQM计算电磁场量对z的偏微 6 分数值近似．以Maxwell旋度方程E。分_垦为例，将式中的偏微分6H。jz用类似式(1)的 DQM近似代替．而其它编微分用传统FDTD算法中的中心差分来近似．得到： E+，(i+；，／，々)： E(％1 J≈)+南 蔓逊：!丛燮蔓塑 口 。■ + 。∑Ⅲ ，，，●●L ，一2 缈厂=■ (3) DQM提高分析效率的指导思想是提高偏微分的计算精度，从而可以增大空间步长的 数值，减少FDTD网格节点的数目。但在FDTD网格边界上必须对吸收边界条件进行差分， 因此大的空间步长将带来误差。另一方面，DQM在靠近边界处的节点上效果相对铰差。 为解决这些问题，本文将采用网格节点两端细分，中间粗分(均匀)的方法(如图l所示)， 对细分网格节点处的场用传统的FD7FD方法或边界条件进行计算，而中间粗分网格节点处 的场则用FDTD结合DQM计算。 DQ分段计算点