[求助]关于葛老师书后的程序的问题

葛老师书后(第二版)的第二个程序中遇到一点问题请教大家 m^rrbU+HM?  
1.real EAC(4,2,0:1),括号里面的数字是什么意思? 4是不是代表有4个边连接边界? >s+TD4OfY  
o$FYCz n  
2,   real EBin(4) !入射波吸收边界条件临时变量 +rA:/!b)Y  
    Ein(IncidentStart)=EBin(4); :VEy\ R>W  
    EBin(4)=EBin(3); =$t  
    EBin(3)=Ein(IncidentStart+1); f-b#F2I  
    Ein(IncidentEnd)=EBin(2); loPBHoE3@H  
    EBin(2)=EBin(1); tQ > IJ  
    EBin(1)=Ein(IncidentEnd-1); ;YK{[$F  
谁能给解释一下这段程序,另外括号里面的数字代表什么意思?? ehCZhi~  
要是在三维里面的话,是不是入射波吸收边界也是这么设定(看了sullivan的程序发现他二维三维都是一样的程序设定入射波吸收边界)  Skn2-8;10  
3. T1=float(I)*cp+(float(Jtmin)-0.5)*sp (底部连接边界处HX .. lp=8RbQYC  
+W-,74A  

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

我自己顶一下,期待高人解答

hoping,waiting高人的解答

xfw)0S  
这个我觉得他程序里面都已经注明了啊。 TfK$tTkM  
比如第二个，是存储临时变量的。 rkW*C'2fz  
括号里面的数字？ zRsG$)B  
A(2)就表示数组A里面的第二个数据，这个，这个。。。很难么，：（

引用第3楼gwzhao于2009-04-06 11:17发表的 : UqJ}5{rt  
[表情] ,T:Uk*Bj  
这个我觉得他程序里面都已经注明了啊。 Y|%s =0M  
比如第二个，是存储临时变量的。 %#5yC|o9Pn  
括号里面的数字？ [=f(u wY>g  
A(2)就表示数组A里面的第二个数据，这个，这个。。。很难么，：（

4KH8dau.fF  
\v|nRn,`-  
 ,这第二问题问的确实......也不知道自己但是脑子怎么想的.......我检讨! 1>5l(zK!9  
第二段程序,问题我想明白了,一开始自己想偏了,因为一维入射波的dt=dx/c,而二维TM波的dt=dx/2c,所以传播同样的距离一维入射波需要2个时间步,传播同样的距离一维与二维相比有一个时间延迟,为了解决这个问题所以要这么设置,不知道这么理解对不对

如果你知道为什么一维吸收边界 &/\Q6$a  
E(n) = E(n-1) | dt = dz/(c)， ;z Qrree#  
E(n) = E(n-2) | dt = dx/(2*c) j0X Jf<  
就更好了。  K}OY!|  

引用第4楼wq_463于2009-04-06 12:09发表的  :
[表情] ,这第二问题问的确实......也不知道自己但是脑子怎么想的.......我检讨! UG~/   
第二段程序,问题我想明白了,一开始自己想偏了,因为一维入射波的dt=dx/c,而二维TM波的dt=dx/2c,所以传播同样的距离一维入射波需要2个时间步,传播同样的距离一维与二维相比有一个时间延迟,为了解决这个问题所以要这么设置,不知道这么理解对不对

引用第5楼gwzhao于2009-04-06 16:25发表的 : [AA*B  
如果你知道为什么一维吸收边界 g.]S5(  
E(n) = E(n-1) | dt = dz/(c)， VJoobu1h  
E(n) = E(n-2) | dt = dx/(2*c) {Y6;/".DM  
就更好了。 i2yE-sgF  

GU,ztO.w3  
根据courant稳定性要求dT<=dx/sqrt(n)/c0,二维和三维我们取sqrt(n)=2,一维取sqrt＝1。 vFx0B?  
R[bI4|t  
E(n) = E(n-1) | dt = dz/(c)  这个我明白 [+2iwfD  
E(n) = E(n-2) | dt = dx/(2*c)  是不是可以理解为当一维入射波引入到二维的TM（TE）波时是这么设置呢？ 6qF9+r&e?  
请师兄指正

有点不对。 5FQtlB9F  
就是一维时，我也可以设置dt = dx / (3 * c)或者 dt = dx / ( n * c)的啊。 V0/O T~gS8  

引用第6楼wq_463于2009-04-06 20:54发表的  :
根据courant稳定性要求dT<=dx/sqrt(n)/c0,二维和三维我们取sqrt(n)=2,一维取sqrt＝1。
E(n) = E(n-1) | dt = dz/(c)  这个我明白 Lcow2 SbH  
E(n) = E(n-2) | dt = dx/(2*c)  是不是可以理解为当一维入射波引入到二维的TM（TE）波时是这么设置呢？ EN/e`S$)  
.......

是这样的, 一维dt = dx/c0 R\^XF8n6/  
                 二维dt=dx/sqrt(2)/c0 N}5'Hk4
+  
                 三维dt=dx/sqrt(3)/c0 c$AwJhl^]  
这是走过一个网格（无论几维）走对角线（最长距离）时的所需要的时间。因此，收敛条件就应该是 Z`{ZV5  
dt<=dx/sqrt(n)/c0，n是维数。那么即使当n=3时，取dt=dx/2/c0也是可以符合收敛条件的。可以看Talove的论文。 X|wg7>kh*`  
我也不知道说得清楚不。

引用第8楼editwolf于2009-04-07 17:12发表的 : 6tKCY(#oO+  
是这样的, 一维dt = dx/c0 `Eq~W@';Q0  
                 二维dt=dx/sqrt(2)/c0 ~f5g\n;  
                 三维dt=dx/sqrt(3)/c0 '#Pg:v_  
这是走过一个网格（无论几维）走对角线（最长距离）时的所需要的时间。因此，收敛条件就应该是 U2TR>0l  
dt<=dx/sqrt(n)/c0，n是维数。那么即使当n=3时，取dt=dx/2/c0也是可以符合收敛条件的。可以看Talove的论文。 2(5<Wj"  
.......

;(I')[R"  
;.EW7`)Z  
你跟我一样还没搞明白gwzhao师兄的问题,师兄的问题是为什么一维的时候也可以设置dt=dx/2c或者dt=dx/3c!! 2n|]&D3V"'  
我还没想明白,正在思考中