波源的设置.pptVIP

第4讲 波源的设置 时域有限差分法 4.1 引言 为了用FDTD法模拟电磁场工程问题，必须在FDTD网格中引入电磁波激励源。 常用的波源种类： 平面波源—用于电磁散射问题； 导波源—用于微波网络参数计算； 电流源或电压源—用于微波电路或天线的激励。 常用的设置方式： 初始条件； 硬源； 电流源； 总场/散射场公式 4.2 入射波源条件 在Maxwell方程的FDTD求解中，入射波源的设置必须满足下列主要要求： 对于平面波源，应保证计算域内部的数值模拟入射波就好象是从计算域外部进入的； 在任何时间步保证所产生的数值平面波在波前保持不变； 允许数值模拟入射波具有任意场极化、任意波形和任意周期； 允许任何散射波无任何干扰和相互作用地通过源所在位置。 4.2 波源作为初始条件插入（1） 在空间网格中每一个电场和磁场分量所在位置上作为初始条件插入入射波。要求所有网格上入射波的所有电场和磁场分量的初值预先给定。选取每一个初始场分量的符号及振幅以给出所希望的波极化和传播方向。 由于Yee算法中电场分量与磁场分量的初值之间相差半个时间步的迟延。所以，当设置入射波初值时，必须考虑由于这个延迟对波传播所造成的波位移。位移取决于在模拟中所使用的时间步长。如果网格填充非色散介质，波速低于c，这些位移要做适当地调整。但如果填充的介质是色散的，则无法直接计算波位移。这时使用后面介绍的其他波源可能更合适。 初始条件波源基本满足入射波源条件。 4.2 波源作为初始条件插入（2） 存在两个基本问题，大大限制了这种波源的用途。 对于长周期脉冲或连续正弦波，必须增加许多自由空间网格以加长计算域，使得初始条件可以被包含，这大大增加了计算机储存和代数运算。 在二、三维空间网格中由这种方法产生的斜入射数值平面波，当波前拖到网格外部边界时由于绕射效应不可避免地产生波前失真。 由于这些限制，初始条件波源目前一般只有一些特殊用途，主要在一维FDTD模拟中。 4.3 “硬”波源 (1) 硬波源可以通过规定FDTD空间网格中电场或磁场分量满足所希望的时间函数进行简单地设置。如：在网格点（is,js,ks)建立Ez的硬源为: 硬波源相当于点源，模拟从源点向外辐射的具有与源函数性质相同的数值波。如果材料结构放在离波源一定距离的地方，辐射的数值波最终传播到这个结构，一部分传输过去，一部分反射。在原理上，时间步进可以连续进行，直到所有的瞬态波消失。 4.3 “硬”波源 (2) 常用源波形有以下三种 正弦波源为时谐源，提供了频率为f0的正弦波。 FDTD模拟时，一般采用Gaussian脉冲信号激励，因为它可以提供宽频带特性。 Gaussian脉冲信号的频谱也是Gaussian函数。 4.3 “硬”波源 (3) 为了减少数值色散误差、满足初值条件并尽可能地利用信号能量，Gaussian脉冲的脉宽和频宽必须受到限制。 定义Gaussian脉冲信号的最高频率为其频谱幅度是最大频谱幅度的a倍时所对应的频率，即 其中0a1，一般a较小，使得大于最高频率的信号能量很小，因此可以忽略大于最高频率的信号。Gaussian脉冲信号的最高频率为 4.3 “硬”波源 (4) 令 ， ， 则 定义常数b为Gaussian脉冲信号初始时刻的幅度与最大幅度之比，即 由于FDTD法把电磁问题作为初值问题处理，因此要求Gaussian脉冲信号初始时刻的幅度很小，即要求b很小。 由于 因此 4.3 “硬”波源 (5) 但是，在连续时间步进时，当反射波到达波源位置，由于波源处总场已被规定，并没有考虑网格中可能的反射波（所以称为硬源），硬源对这些反射波就会造成寄生的非物理的再次反射。当规定了一个表面上的总场而没有考虑反射场值时，寄生反射就一定会发生。 解决这一问题的一个简单方法是当反射波到达波源位置时去掉波源。这就要求源离散射体足够远，使散射波到达源位置时源入射波已基本为零，所以这种方法不适合正弦入射波源。 不过硬源在某些结构中也是非常合适的，如探针激励就是典型的硬源，而且在二、三维模拟中，硬源引入的反射比一维中要小得多。 4.4 电流源 在Maxwell方程中引入电流密度J： （4-1） 则FDTD公式为 其中， 对于电流源I， FDTD公式中源项S也可以看成是电场源。 电流源不会造成寄生反射。 4.5 总场/散射场公式(1) 总场/散射场方法是基于Maxwell方程的线性特性和下列电磁场的分解 (