第二章时域有限差分法IV-三维FDTD.ppt

* 3. 三维FDTD 重写Maxswell方程组 4.1 自由空间FDTD公式 （4-1-a） （4-1-b） （4-1-c） 二维Maxswell方程组标量方程 三维Maxswell方程组标量方程 (4.2a) (4.2b) (4.2c) (4.2d) (4.2e) (4.2f) （4-1-a） （4-1-c） 对时间和空间差分后，以(4.2c)和(4.2f)为例，迭代公式 dx(i,j,k) = dx(i,j,k) + 0.5 *( hz(i,j,k) - hz(i,j-1,k) ... - hy(i,j,k) + hy(i,j,k-1)); dy(i,j,k) = dy(i,j,k) + 0.5 *( hx(i,j,k) - hx(i,j,k-1) ... - hz(i,j,k) + hz(i-1,j,k)); dz(i,j,k) = dz(i,j,k) + 0.5 *( hy(i,j,k) - hy(i-1,j,k) ... - hx(i,j,k) + hx(i,j-1,k)); hx(i,j,k) = hx(i,j,k) + 0.5 *( ey(i,j,k+1) - ey(i,j,k) ... - ez(i,j+1,k) + ez(i,j,k)); hy(i,j,k) = hy(i,j,k) + 0.5 *( ez(i+1,j,k) - ez(i,j,k) ... - ex(i,j,k+1) + ex(i,j,k)); hz(i,j,k) = hz(i,j,k) + 0.5 *( ex(i,j+1,k) - ex(i,j,k) ... - ey(i+1,j,k) + ey(i,j,k)); 迭代公式 在三维基本FDTD程序和二维FDTD程序中，都是在问题空间的中心设置简单点源。但在三维仿真中，由点源产生的电场向外传播的过程中，与距离的平方成反比衰减。相反，我们可用电偶极子天线作为源。如下图所示。 （1）在两臂空隙设置一电 场值。 （2）依照安培环路定律 设置围绕空隙的磁场值。 设置激励的方法 T=50 T=20 T=40 T=70 4.2 三维PML 将二维问题的PML扩展至三维 其中： 由第三章知 %%%%%%%%%%% Loop for D_z %%%%%%%%%%% for i=2:IE for j=2:JE for k=1:ka curl_h = hy(i,j,k) - hy(i-1,j,k)... -hx(i,j,k) + hx(i,j-1,k); idzl(i,j,k) = idzl(i,j,k) + curl_h; dz(i,j,k) = gi3(i)*gj3(j)*dz(i,j,k)... +gi2(i)*gj2(j)* 0.5 *... (curl_h+ gk1(k)*idzl(i,j,k)); end end end 关于Dz的迭代程序，分三段计算以节约存储资源 for i=2:IE for j=2:JE for k=ka+1:kb+1 curl_h = hy(i,j,k) - hy(i-1,j,k)... -hx(i,j,k) + hx(i,j-1,k); dz(i,j,k) = gi3(i)*gj3(j)*dz(i,j,k)... +gi2(i)*gj2(j)* 0.5 *curl_h; end end end for i=2:IE for j=2:J