[转载]时域有限差分法（FDTD）系列之二

进行一个具体计算，假设网格的大小是L*M*N, 方向分别是x, y, z。当计算y=M这个面上的电场时，是需要y=M+1上的磁场的数值。也就是说计算里面一层的数值时需要用到外面一层的数值，但总不能有无限的层，因为内存和计算时间是有限的，我们总不能在计算机上模拟出无限的空间行无限的计算。因此，在计算y=M的电场时就要用一种特殊的方法来计算，称之为“吸收边界条件(ABC)”,简称为“边界条件”（在电磁场这门学科种有多种边界条件^-^）。 wOQ-sp0q0  
g{wIdV  
       实际情况中，电磁场是要向无限的空间传播的，举个最简单的例子：空间中一个点源辐射电场的情况。在计算中，当电场传播到人为设定的边界上时，如果我们强制电场为零，就相当于在边界上放置了导体，就会虚拟出了反射，这可以称为“数值反射”。但如果我们设计一种算法使得电场到了边界时被吸收了，这就虚拟出了透射的情况，模拟出了在计算中无限空间的情况（这就好像微波实验室的吸波材料的作用：在实验室模拟出无限空间）。 ;V]EF  
^/\OS@CT\  
       1981年G. Mur在[1]中给出了一种较好的ABC。它的原理是：根据数学推导，得到：当微波经过某个面时不被反射的数学条件，即边界上的数值满足这个条件时理论上就是没反射的。由于这个数学条件没法直接离散，需要对它取近似。这种方法在计算边界上的某一点的值时，要用到边界里面一层的对应点和周围的点的数值。根据近似的程度不同，所用周围的点的情况也不同。具体来说，利用二阶近似计算边界y=M .. CL%?K<um  
*m:
h0[[J  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册