应用傅立叶伪谱时域算法与高阶时域有限差分方法的混合算法到25维问题中研究.pdfVIP

2005。全国微波毫米波会议论文集 应用傅立叶伪谱时域算法与高阶时域有限 差分方法的混合算法到2．5维问题中 史琰 梁昌洪 (西安电子科技大学电子工程学院西安710071) 摘要一本文将傅立叶伪谱时域算法(FPsTD)与高阶时域有限差分算法(FDTI))相混合来 求解2,5维的电磁问题。对于均匀的或者均匀变化的媒质，FPSTD算法仅需要每最小波 长2个采样点就可达到所需要的精度。然而对任意变化的不均匀媒质而言，该算法的精 度就会大大地降低。在这种情况下，FDTD是一个很好的选择。混台后的算法不仅维持 了各自算法的简单性而且精度也得到了提高。文中所给的数值结果充分验汪了混台算法 的性能。 I．引言 近年来，高阶和谱域算法作为一种精确有效的计算技术越来越引起人们的注意。傅立叶 匀变化的媒质，该算法仅仅需要每最小波长2个采样点就可以达到所需要的精度。对任 意变化的不均匀媒质而言。该算法的精度会大大地降低。尽管时域有限差分算法 地处理各种复杂的电磁媒质。因此将FPSTD和FDTD混合来处理在某一维或二维方向 上是均匀的而其它方向上是不均匀的问题是十分有益的。Leung和Chan[3]给出了二阶 FDTD算法和FPSTD算法的混合算法。在本文中，我们将这种混和算法推广到高阶FDTD 算法和FPSTD算法相混合，并将所得的算法应用到2．5维的电磁问题中。 II．傅立叶伪谱时域算法与高阶时域有限差分方法的混合算法 ’、2．5维Maxwell方程 在一些应用中，例如接地雷达探测地面下的物体，所研究的媒质一般仅仅在两个方I自J(x 方向和p方向)上是不均匀的，在第三个方向(z方向)上是均匀的，但是电磁场量仍然是 (五Y，21的函数。对于这类问题我们称为2,5维电磁问题。在这种情况下，我们通过正旋 变换和余旋变换可以将Maxwell方程进行简化。将电磁场量沿z方向分解成奇数部分和 偶数部分，即雷：营。+雷。，豆=詹。+膏。。例如 E。x，Y，女：，t)=【层8b，Y，鼻t)cos(k：z弦 (1) 豆。(x,y…k (2) f)=r豆。x，Y，z，f)sin伍：z迎 其中，k是个傅立叶积分参数。通过傅立叶变换后，标准的Maxwell可以分成两类： 416 2005’全国微波毫米波会议论文集 Groupl：第一类是由电流源0：，E，彤)所产生的场分量幢，E，茸，疗；，茸，膏；)。 a “I却 ‘…J 警一去【等～：髟] 2‘ 譬：一土k茸一譬] a ∥I 缸J 旃；一1f嘭啦] Ot ∥l缸 咖J 钟 s’ 警={[等批蚓一吉z sl钞 。’J 鲁=舭弘刳一!。y，o 重：三f堡一譬1一土z ㈣ Ot 占【Ox 砂J 5。 Group2：第二类是由电流源(鬈，，；，f)所产生的场分量瞎：，髟，孽，疗：，膏；，膏；)。其场方 程可以由第一类的场方程中k．换成一≈，，奇数部分与偶数部分互换得到。下面我们仅仅 将傅立叶伪谱时域算法和高阶时域有