动用全波离散镜象理论计算二维空域格林函数(一)%3a理论.pdfVIP

运用全波离散镜象理论计算 二维空域格林函数(一)：理论 周桂蓉冯宁宁方大纲 南京理I大学l毡光学院蓬米波技术研究室，南京，210094 【摘要】：本文先同顾丁全波离散镜象(FWDI)理论的发展过程．对于虽早的 一层近似方法与后来发腮的两层近似方法作了比较研究，并且从理论上解释了在 计算二维格林函数中必须抽取表砸波极点的机理。 l。引言 在用积分方程方法分析分层介质结构问题时要遇到高振荡慢收敛的 Sommerfeld积分，采用传统的数值积分方法将花费较长的CPU时间。全波离散 镜象(FWDI)理论(也称复镜蒙理论)的提出划解决了耩确快速计算三维空域 格林豳数的问题f11，，{：上|．经过后来的发展已在二维问题中得到，广泛的应用f2卜 f4】。而在二维结构的传播问题中有关FWDI的应用在国际上目前还非常少见。 最近，FWDI被用于分析多导体的传播闷题【5]．但所分析的结构仅限于单层介质 并且没有对二维问题必须抽取表面波极点的机理作出解释。本文将先对最早的一 层近似方法【1]同后采发展的两屡近似方法【4】作理论和数值实验上舶比较。再对 一：维计算中事先抽取表面波极点贡献的必要性作出了理论解释。 2．两种近似方法的比较 最早的一层近似方法将空域格林函数分戊三部分：准静项、表面波项和复镜 析积分得到；表面波项则由谱函数极点的留数积分得到；剩余的谱经过Prony方 法沿变形路径拟台，再运用sT得到闭式空域解，印复镜象项。而后来发展的两 层近似方法有两点改进：(1：、以第一层近似的数值抽取代替解析准静抽取，由于 准静项的解析空域表达式仪征～些简单情况诸如单层介质的情形下才存在，从而 给这一方法的广泛应用带来了困难。而在两层近似方法中第一层近似是采用与复 镜象项类似的数值方法得到，因此比较任意，便于推广到多层介质结构。图l显 OHz．)的谱函数分别用解析抽取准静项和第一层近似数值抽取情况的比较，结 果表明两种方法都使谱函数很快由原来的o．85附近衰减至o．02附近，从而加快 了后续计算的速度，从图中还可以看出，虽然在物理本质上这两种抽取是对应的， 方法拟合谱函数，对于由表面波引起的近谱的剧烈变化也能准确拟台，从而不必 玉计算极点『句留数．这一做法的止确性已住三维问题中得到验证【4】。但是这种 处理力洼在汁算一二维问题时是父效的。阳2显示的是对一个两层介质结构(s．=3，9． d=0．5u OHz)的标量位二维空域格林 m，s。=43．0，s。=28．u．h=O．5ram，，=50 函数的计算情况，结果表明直接HJ两层近似而不抽表面波极点将使结果在远区偏 离实际结果(即数值积分结泉)较远；相反如粱事先抽取了极点贡献就能使结果 --66-- 准确。关丁-在三维问题中1I常有效阿层近似的．方法在二维问题火效的机理将在F 一部分讨论。 屯，岛 LogIoko{x—x’l 图1对标量位谱函数的准静项进行解析 舀2表面波拔点羹献对标量位空域格林 抽取和数值抽取的比较。 函数的计算的影响。 3。 表面波极点的影响： 点再沿实轴到k。，然后沿原路径返回。根据复变函数理论，沿图3(a)中原始积 分路径C。和变形拟台路径C，积分是完全等价的，因为在这两条路径之间没有 附加极点。 所以即使不抽极点，远场的表面波信息还是能够卷入格林函数的计 算中。而对丁二维情形则不同，此时．|}二=昧一k：一∥2并且积分是对k，沿 r——?———一 绕回到一00，此时极点的主值积分贡献就未能卷入到原始积分路径C。中。为了 更清楚地揭示极点抽取的必要性，我们在图4(a)和(b)中分别画出了在二维情形 tin(k：o)J lm(k：o)· Co^o 五o