一种提高时域有限差分法计算精度的新算法.pdfVIP

一种提高时域有限差分法计算精度的新算法 一种提高时域有限差分法计算精度的新算法 李永明俞集辉谭邦定 (重庆大学电气工程学院， 重庆市沙坪坝区400044) 摘 要 介绍了一种利用外推方法来提高时域有限差分法计算精度的算法，给出了这种外摊法应用到时域有限差分法 中的具体计算格式。通过数值实验可知，这种应用到时域有限差分法中的外推法只需进行简单的外推计算就能比 ‘ 较好的提高计算精度。 ． 关键词：外推法，时域有限差分法 l引言 1966由feell】提出的时域有限差分法(Finite．Difference 和电磁兼容问题等方面起了巨大的作用。由于这种方法的时域特点，电磁场作用的细节过程可以 被详细的描述。但这个方法也受到一些制约。比如计算时间受稳定性条件的限制。并且如果提高 计算精度，就必须大大加密计算区域的剖分等。为了解决这些矛盾，先后提出了若干种方法。如 2阶时间、4阶空间近似鸭局部网格加密法【3l’但它们的计算格式都比较复杂。例如，用局部网 格加密法，那么在不同的网格交界处的场量的衔接关系就十分复杂且有较大的近似。另外，对计 Namikil41等人利用了ADI．FDTD取消了一 算的时间步长和空间步长的稳定性条件，虽然Takefumi 定的约束，但计算格式过于繁琐。 ’ 本文基于外推法提出了一种提高时域有限差分法计算精度的新算法，在这种算法中，只需通 过利用原计算区域加密剖分一次前后两次的计算解，再通过简单的外推计算，便可得到一个新的 解一一外推解【外。数值实验结果表明，这种外推法得到的解能有效的提高时域有限差分法的计算 + 1 精度。 、 2外推法的概念 外推法是一种通过组合同一问题的相邻两次计算，得到更精确解的一种方法。我们通过例子 来介绍一下外推法。对一圆内接正多边形，如果以口II表示内接于一个半径为l的正r／边形的边长， ’ 那么圆周率的近似值为： 1 万一。虿阳n 只要我们取足够多的边数，就可以很精确地估计圆周率。但是，如果我们利用一个内接多边 形及在此基础上加倍得到的多边形的近似值，如磊和况n，作一个线性组合： 7／-。=妻(4月。：py／2。一万。)n一万n夕 j 可大大提高计算精度【6J。 第十属全国电工数学学术年会论文集 ■●■■