[已解决]二阶mur边界吸收条件为什么发散了

我是fdtd新手 ，编了一个电磁波在矩形波导中传播的fdtd代码。开始的时候截断边界也是PEC，虽然有点发散，但在1000个时间步里还可以忍受。现在我学习二阶mur边界吸收条件，于是在原程序的基础上在截断处加了二阶mur，但是当电磁波一到截断边界就发散相当厉害，值变得很大。如图： [_wenlkm  
第165步，电场最大值为1.365 i2(lqhaP  
15tT%TC  
第264步，电场最大值为3.650 T9=55tpG9  
5>*~1}0T  
第363步，电场最大值为285.3 :ZP`Y%dt'
 
{|ChwM\x  
,CA3Q.y>|  
以下是部分二阶mur边界条件代码，没有包含棱边 P gK> Z,  
do i=x2,x2,1 _vgFcE~E@  
     do j=1,(ny-1),1 Y
;OqdO  
        do k=1,(nz-1),1 j9]H~:g$d  
           ey(i,j,k) = -ey_scratch(i-1,j,k) + (LIGHT_SPEED*dt-dx)*(ey(i-1,j,k)+ey_scratch(i,j,k))/(LIGHT_SPEED*dt+dx)+2*dx*(ey_current(i,j,k)+ey_current(i-1,j,k))/(LIGHT_SPEED*dt+dx)+dx*LIGHT_SPEED*dt*LIGHT_SPEED*dt*(ey_current(i,j+1,k)-2*ey_current(i,j,k)+ey_current(i,j-1,k)+ey_current(i-1,j+1,k)-2*ey_current(i-1,j,k)+ey_current(i-1,j-1,k-1))/(2*dy*dy*(LIGHT_SPEED*dt+dx))+dx*LIGHT_SPEED*dt*LIGHT_SPEED*dt*(ey_current(i,j,k+1)-2*ey_current(i,j,k)+ey_current(i,j,k-1)-ey_current(i-1,j,k+1)-2*ey_current(i-1,j,k)+ey_current(i-1,j,k-1))/(2*dz*dz*(LIGHT_SPEED*dt+dx)) i*-L_!cc:  
               ez(i,j,k) = -ez_scratch(i-1,j,k) + (LIGHT_SPEED*dt-dx)*(ez(i-1,j,k)+ez_scratch(i,j,k))/(LIGHT_SPEED*dt+dx)+2*dx*(ez_current(i,j,k)+ez_current(i-1,j,k))/(LIGHT_SPEED*dt+dx)+dx*LIGHT_SPEED*dt*LIGHT_SPEED*dt*(ez_current(i,j+1,k)-2*ez_current(i,j,k)+ez_current(i,j-1,k)+ez_current(i-1,j+1,k)-2*ez_current(i-1,j,k)+ez_current(i-1,j-1,k-1))/(2*dy*dy*(LIGHT_SPEED*dt+dx))+dx*LIGHT_SPEED*dt*LIGHT_SPEED*dt*(ez_current(i,j,k+1)-2*ez_current(i,j,k)+ez_current(i,j,k-1)-ez_current(i-1,j,k+1)-2*ez_current(i-1,j,k)+ez_current(i-1,j,k-1))/(2*dz*dz .. x;p7n2_  
52>,JHq  

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

你用啥编的,Matlab吗

引用第1楼chenxl0426于2008-10-08 12:07发表的  : \@}G'7{  
你用啥编的,Matlab吗

;tZ;C(;
<  
我用fortran编的

Ey中的第二个系数怎么会有+dx，应该是 -dy 吧？

顺便问一句，你是用VisIt绘制的图吗？电场值在什么范围算正常？你的激励源是怎么加入的？

引用第3楼air81666于2008-10-11 12:31发表的  : PNs*+/-S  
Ey中的第二个系数怎么会有+dx，应该是 -dy 吧？

jAcrXB*  
T!ww3d  
哪儿？葛德彪老师的书上就是这样的啊。

引用第4楼air81666于2008-10-11 19:49发表的  : LS$zA>:  
顺便问一句，你是用VisIt绘制的图吗？电场值在什么范围算正常？你的激励源是怎么加入的？

"~mY4WVG  
2?{'(iay  
不是哈。 ~F^tLi!5  
电场值的正常范围在几千个时间步里面，应该在激励源的电场幅值左右。因为mur边界吸收条件也不能完全吸收波，所以会有小的反射。 ^IKT!"J&?  
我的激励源很简单，就是在x=0的界面令ez等于一个正弦函数，ey和ex都为0。

问题解决了！ '$c9S[  
我总结了两点原因： v<t?t<|J  
1. 我所参考的书上误把一个加号写成了减号。我开始编写的时候的没有注意到这个问题，所以mur边界吸收条件出现错误。 _!E&%=f  
2. 之前编写程序的时候，边界棱边没有处理好，进一步造成部分发散。 l2kGFgc  
谢谢大家的帮助！

有QQ吗？在QQ上聊吧，我的：86693550

引用第7楼shawny于2008-10-12 12:12发表的  : Pmo<t6  
问题解决了！ #G.eiqh$a  
我总结了两点原因： p<{P#?4 g  
1. 我所参考的书上误把一个加号写成了减号。我开始编写的时候的没有注意到这个问题，所以mur边界吸收条件出现错误。 ) $_1U!z  
2. 之前编写程序的时候，边界棱边没有处理好，进一步造成部分发散。 I8YUq  
谢谢大家的帮助！

5 ty2e`~K  
4qz+cB_  
是哪里的符号错了？我参考的是葛德彪老师的那本书