[求助]1维FDTD TFSF 程序源码以及求教

oned_tfsf_fdtd.rar (1 K) 下载次数:16     根据Allen Taflove 的TFSF一维FDTD 编的程序，所用入射场激励为正弦函数，幅度为1；总场区的Ez分量和入射场看上去比较吻合； 3>7{Q_5  
可是没有加散射源，理论上散射场部分的常量应该为0，在运行我的程序得到的结果中， 在散射场部分的值幅度在10^-4， 请问这种情况正常么？ k;cX,*DIn  
hu0z 36  
结果以及程序见附件。 hR|xU
p  
9yTdbpY  
.. AYY(<b  
QObVJg,GD  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

没人回复么？ I&i6-xp  
自己顶啊顶 'W9[Vm  
_\IA[-C+O  
Kp7)my  

引用楼主lwfz于2010-10-21 10:01发表的 1维FDTD TFSF 程序源码以及求教 : hunlKIg  
[图片]    根据Allen Taflove 的TFSF一维FDTD 编的程序，所用入射场激励为正弦函数，幅度为1；总场区的Ez分量和入射场看上去比较吻合； 6v scu2  
可是没有加散射源，理论上散射场部分的常量应该为0，在运行我的程序得到的结果中， 在散射场部分的值幅度在10^-4， 请问这种情况正常么？
结果以及程序见附件。
....... !.#g  

rkWy3X{%2<  

散射场有泄露是正常的，但一维的情况，感觉这个在散射场有点大。。。

描述:你的附件程序运行结果

图片:untitled.jpg

z`dnS]q9  
uE#i3( J  
你程序运行结果与描述的结果对不上。。 Bq
,Pk5b  
3[kl` *`  
程序中有迭代错误 jlD3SF~2  
for i=2:99 dr"@2=Z  
    Hy(i-1)=Hy(i-1)+dt/(mu*dx)*(Ez(i)-Ez(i-1)); 7Z81+I|&8  
    Ez(i)=Ez(i)+dt./(eps*dx)*(Hy(i)-Hy(i-1)); MLDAr dvK  
end 7B)@ aUj$  
Hy(99)=Hy(99)+dt/(mu*dx)*(Ez(100)-Ez(99))-dt/(mu*dx)*Ez_inc100;      % incident field Ez_inc at 100 needed [Fd[(  
Ez(100)=Ez(100)+dt/(eps*dx)*(Hy(100)-Hy(99))-dt/(eps*dx)*Hy_inc99;  % incident field Hy_inc at 100 needed 0|rd
I,z  
M?ElD1#Z  
在上面的for 循环中， Ez(99) 的成为当前时刻的值， 在求 Hy(99)的值时用到的应该是 Ez(99)的前一次迭代求之后的值，而不是当前的Ez(99)的值。 R~dWblv  
同样Ez(100)的求解也有同样的问题， Ez的求解应该用当前迭代得到的Hy的值求，而程序中求解Ez(100)用的Hy(100)是前一次迭代之后的值。 Z= pvoTY  
(b.Mtd  
在写迭代时，最好把电场和磁场分开写，因为电场和磁场有1/2的时差。 BJZGQrsz  
.MxMBrM  
:|ytw=3>  
K}LF ${bS  

谢谢指点 4--[.j*W  
我按照你所说的对H和E分别进行了更新，程序运行完之后， 输入plot（Ez_save）是总场区的Ez分量随时间的变化， plot（Ez_save2）是散射场区的Ez分量随时间的变化。不过貌似在散射场区还是有一定的泄漏

oned_tfsf_fdtd_2.rar (1 K) 下载次数:26 =-|,v*  
"AUHe6Yv  
入射场也使用FDTD的形式计算，在所求的总场边界处引入入射场。你看看上面的程序。 xnD"LK  
试一下这样的方式引入入射波到总场区，误差会变小了。。 $fl+l5?9  

我也是学FDTD编程的，学习一下！

请问下，因为这个FDTD里面没有加入吸收边界，是不是计算时间比较多的时候，E H 来回的反射叠加，造成场量不断的增加？ ;fe~PPT  
k$3Iv"gbx  
在有PML的情况下，  我在计算Ez_INC 和Hy_INC的时候是不是也得加个吸收边界？  这样的话，一个程序里面就得有两个计算循环，一个 dwJnPJ=z  
Ym|%ka  
计算入射场，一个算E H ？ 师兄能不能帮我解答下？

1）请问下，因为这个FDTD里面没有加入吸收边界，是不是计算时间比较多的时候，E H 来回的反射叠加，造成场量不断的增加？ S$n?  
对这种情况，对边界什么都不处理，其实就相当于在边界全反射，类似理想金属边界条件。 0+T:};]  
从能量角度来看，如果你对这个密封腔体不断的持续加入能量，又没有能量损失，那波幅肯定会越来越大。 \+j:d9?  
NA=m<n#  

你可以去看看sullivan 的书上第四章写的一个二维的程序例子，程序中也是对入射场单独计算的。 |K" nSXzk  
希望对你有帮助