在对电磁场模型数值求解中,除了需要在足够的网格空间中模拟煤质,设置合理有效的激励波源也十分重要,需要尽可能地使源的特性与实际物理模型性质一致。
从源随时间的变化看有两类激励源:一是随时间周期变化的时谐场源,另一类是对时间呈脉冲函数形式的波源。从空间分布来看,有面源、线源、电源等。
几种常见的时谐源和脉冲源:
-
时谐场源
在早起的FDTD法计算中,由于缺乏合适波导中色散严重的宽带吸收边界条件,波导中的不连续问题往往采用这种单频激励。由于任何波模都可以展开为平面波谱,一个算法若对平面波不稳定,则对其他波模也不稳定,因此以平面波为对象来进行研究是一种行之有效的方法。网格空间按真空波长 λ 0 \lambda_0 λ0(煤质中波长为 λ n \frac{\lambda}{n} nλ,其中 n n n为介质的折射率)设定,再次空间模拟真空波长为 λ 1 \lambda_1 λ1的光波。
假设入射场为:
E i ( t ) = 0 ( t < 0 ) E i ( t ) = E 0 sin ( ω t ) ( ( t ≥ 0 ) E_i(t)=0(t<0)\\ E_i(t)=E_0\sin(\omega t)((t\geq0) Ei(t)=0(t<0)Ei(t)=E0sin(ωt)((t≥0)
上式所示激励源情况下达到时谐场的稳态,通常需要3-5个周期。对于散射问题所需的周期数还与散射体大小及形状有关。例如,对于具有凹腔结构的物体,Taflove等指出,达到稳定状态所经过的周期大约等于所模拟散射体结构的Q值。为缩短稳态建立时间,减小冲击效应,可以引入开关函数,如采用升余弦函数。 -
高斯脉冲
高斯脉冲函数的时域形式为:
E i ( t ) = e − 4 π ( t − t 0 ) 2 τ 2 E_i(t)=e^{-\frac{4\pi(t-t_0)^2}{\tau^2}} Ei(t)=e−τ24π(t−t0)2
其中 τ \tau τ为常数,决定了高斯脉冲的宽度。脉冲的峰值出现在 t = t 0 t=t_0 t=t0时刻,如图(a)所示。
上式的傅立叶变换为:
E i ( f ) = τ 2 e − j 2 π f t 0 − π f 2 τ 2 4 E_i(f)=\frac{\tau}{2}e^{-j2\pi ft_0-\frac{\pi f^2 \tau^2}{4}} Ei(f)=2τe−j2πft0−4πf2τ2
其频谱如图(b)所示,其中负频率部分已经去掉。通常可取 f = 2 τ f=\frac{2}{\tau} f=τ2为高斯脉冲的频宽,这时频谱为最大值的4.3%;在 f = 1 τ f=\frac{1}{\tau} f=τ1时为最大值的45.6%;大约在 f = 1.7 τ f=\frac{1.7}{\tau} f=τ1.7时为最大值的10%。
高斯脉冲中包含了丰富的频谱,其波形在时间按分布上是光滑的,他比较适合作为FDTD法研究波导问题的初始激励信号。但高斯脉冲并非理想的激励源。首先他的激励效率不高,有高斯脉冲的频谱函数可知,他的主要能量分布在零频率(直流)和低频段,而实际中人们感兴趣的频段往往是微波毫米波或跟高的频段,但在这些频段上高斯脉冲的能量分布比零频率的能量分布往往小几个数量级,这对有效的激励所研究的电磁系统非常不利。其次,由于极宽的频谱分布在导波系统中往往造成主模式截止模式和高阶模式并存的局面,特别是在波导中能量分布较大的截止模式的存在会严重干扰正常的FDTD计算,为例摆脱这些截止模式的干扰,常常需要加大计算区域衰减到他们足够小的程度。但实际计算表明,当计算区域太小时,这些截止频率分量将造成计算的不稳定性甚至计算失败。 -
升余弦脉冲
升余弦脉冲的时域形式为:
E i ( t ) = 0.5 [ 1 − cos ( 2 π t / τ ) ] ( 0 < t ≤ τ ) E i ( t ) = 0 ( o t h e r ) E_i(t)=0.5[1-\cos(2\pi t/\tau)](0<t\leq \tau)\\ E_i(t)=0(other) Ei(t)=0.5[1−cos(2πt/τ)](0<t≤τ)Ei(t)=0(other)
τ \tau τ为脉冲底座宽度。时域波形如图(a)所示。
频域形式为:
E i ( f ) = τ e − j π f τ 1 − f 2 τ 2 sin ( π f τ ) π f t E_i(f)=\frac{\tau e^{-j\pi f \tau}}{1-f^2\tau^2}\frac{\sin(\pi f\tau)}{\pi ft} Ei(f)=1−f2τ2τe−jπfτπftsin(πfτ)
-
微分高斯脉冲
将高斯脉冲求导后得到微分高斯脉冲函数:
E i ( t ) = t − t 0 τ e − 4 π ( t − t 0 ) 2 τ 2 E_i(t)=\frac{t-t_0}{\tau}e^{-\frac{4\pi(t-t_0)^2}{\tau^2}} Ei(t)=τt−t0e−τ24π(t−t0)2
其优点是不含零频率分量。波行和频谱图如下。
频域形式为:
E i ( f ) = − j τ 2 f 8 e − j 2 π f t 0 − π f 2 τ 0 4 E_i(f)=-\frac{j\tau^2f}{8}e^{-j2\pi ft_0-\frac{\pi f^2\tau_0}{4}} Ei(f)=−8jτ2fe−j2πft0−4πf2τ0
- 调制高斯函数
调制高斯函数的时域表达式为
E i ( t ) = − cos ( ω t ) e − 4 π ( t − t 0 ) 2 τ 2 E_i(t)=-\cos(\omega t)e^{-\frac{4\pi(t-t_0)^2}{\tau^2}} Ei(t)=−cos(ωt)e−τ24π(t−t0)2
上式右边第一项为基波表达式,中心频率为 f 0 = ω 2 π f_0=\frac{\omega}{2\pi} f0=2πω;第二项为高斯函数形式, t 0 t_0 t0通常取基波的 2 1 4 2\frac{1}{4} 241个周期,即 t 0 = 9 π 2 ω t_0=\frac{9\pi}{2\omega} t0=2ω9π。下图为调制高斯脉冲的时域波行和频域波形。其中高斯脉冲的频谱为:
E i ( f ) = τ 4 e − π ( f − f 0 ) 2 τ 2 4 e − j 2 π ( f − f 0 ) t 0 + τ 4 e − π ( f + f 0 ) 2 τ 2 4 e − j 2 π ( f + f 0 ) t 0 E_i(f)=\frac{\tau}{4}e^{-\frac{\pi(f-f_0)^2\tau^2}{4}}e^{-j2\pi(f-f_0)t_0}+\frac{\tau}{4}e^{-\frac{\pi(f+f_0)^2\tau^2}{4}}e^{-j2\pi(f+f_0)t_0} Ei(f)=4τe−4π(f−f0)2τ2e−j2π(f−f0)t0+4τe−4π(f+f0)2τ2e−j2π(f+f0)t0
由上式可知,调制高斯脉冲的频谱与高斯脉冲的频谱相比向零频率点两侧移动来 f 0 f_0 f0,如果要使直流分量为零,有效频谱中心位于 f 0 f_0 f0,则可采用下列调制高斯波作为激励波:
E i ( t ) = E 0 e − 4 π ( n − n 0 ) 2 n d e c a y sin ( 2 π f 0 ( n − n 0 ) Δ t ) E_i(t)=E_0e^{-\frac{4\pi(n-n_0)^2}{n_{decay}}}\sin(2\pi f_0(n-n_0)\Delta t) Ei(t)=E0e−ndecay4π(n−n0)2sin(2πf0(n−n0)Δt)
http://www.docin.com/p-721135108.html