[已解决]三维空间中在边界面处加入二阶mur吸收边界条件时，时间步怎么处理？谢谢了

      计算具体一个点时，要涉及到他周边的9个点，此点和不是边界面上的6个点容易处理，但是在边界面上的其它四个点怎么知道此时的时间步是多少？和具体循环过程有关，每个点按情况考虑，和计算没计算有关？这样的话，紧贴棱边部分（也是紧靠一阶mur吸收的4条边界上）的点又该怎么处理？ {E3;r7  
        本来都快编完程序了，突然想到这个问题了！钻牛角了！请教各位了！谢了！
UnJi& ~O  
$&bU2]  

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

简单介绍一下 je，ke 为计算空间在y和z轴上的大小，P，Q位置关系：P0为计算点，P1234 为在边界面上P0周围的点，Q01234 为紧邻边界面上的点！ kW*W4{Fth  
pZNlcB[Qn-  
这样编程计算，是不是部分P点计算的是上上时间步的值？请教了！ 9&?tQ"@x  
p1[WGeV  
%     p0_ez_x=zeros(2,je,ke);              % 前后面P0处 ez  t-1步的值 x`&P}4v0  
%     q00_ez_x=zeros(2,je,ke);             % 前后面Q0处 ez  t-1步的值 w+:+r/!g  
%     q0_ez_x=zeros(2,je,ke);              % 前后面Q0处 ez  t步的值 :Map,]]B_  
%     p1_ez_x=zeros(2,je,ke);              % 前后面P1处 ez  t步的值 dA3`b*nC  
%     p2_ez_x=zeros(2,je,ke);              % 前后面P2处 ez  t步的值 >B|ofwm*  
%     p3_ez_x=zeros(2,je,ke);              % 前后面P3处 ez  t步的值 >rQ)|W=i  
%     p4_ez_x=zeros(2,je,ke);              % 前后面P4处 ez  t步的值 fh/)di  
%     q1_ez_x=zeros(2,je,ke);              % 前后面Q1处 ez  t步的值 @pyA;>U  
%     q2_ez_x=zeros(2,je,ke);              % 前后面Q2处 ez  t步的值 oW` *FD  
%     q3_ez_x=zeros(2,je,ke);              % 前后面Q3处 ez  t步的值 JC1BUheeb  
%     q4_ez_x=zeros(2,je,ke);              % 前后面Q4处 ez  t步的值 ~8S4Kj)%  
    for j=3:je-2                                                              %前面 EZ 吸收边界 kr &:;  
        for k=2:kb-1 >-o?S O(M,  
            ez(1,j,k)=-q00_ez_x(1,j,k)+(ez(2,j,k)+p0_ez_x(1,j,k))*(cz*dt-dx)/(cz*dt+dx)-(q0_ez_x(1,j,k)+ez(1,j,k))*2*(cz*dt-dx)/dx+... @5uyUSt]  
                (p1_ez_x(1,j,k)+p2_ez_x(1,j,k)+p3_ez_x(1,j,k)+p4_ez_x(1,j,k)+q1_ez_x(1,j,k)+q2_ez_x(1,j,k)+q3_ez_x(1,j,k)+q4_ez_x(1,j,k))*(cz*dt)^2/(2*dx*(cz*dt+dx));   t/4/G']W  
            p0_ez_x(1,j,k)=ez(1,j,k); r)*_,Fo|  
            q00_ez_x(1,j,k)=p0_ez_x(1,j,k); |wl")|b%  
            q0_ez_x(1,j,k)=ez(2,j,k); lBYS>4~  
            p1_ez_x(1,j,k)=ez(1,j+1,k); 2l;ge>DJ  
            p2_ez_x(1,j,k)=ez(1,j,k-1); ^+YGSg7  
            p3_ez_x(1,j,k)=ez(1,j-1,k); (5CX*)R  
            p4_ez_x(1,j,k)=ez(1,j,k+1); (]GY.(F{  
            q1_ez_x(1,j,k)=ez(2,j+1,k); oQE_?">w  
            q2_ez_x(1,j,k)=ez(2,j,k-1); IrVM|8vT3  
            q3_ez_x(1,j,k)=ez(2,j-1,k); 2g`uC}  
            q4_ez_x(1,j,k)=ez(2,j,k+1);     5wVi{P5+  
        end aTsy)=N  
    end

引用楼主lizi于2008-09-10 11:08发表的 三维空间中在边界面处加入二阶mur吸收边界条件时，时间步怎么处理？谢谢了 : cUX]tiC0  
      计算具体一个点时，要涉及到他周边的9个点，此点和不是边界面上的6个点容易处理，但是在边界面上的其它四个点怎么知道此时的时间步是多少？和具体循环过程有关，每个点按情况考虑，和计算没计算有关？这样的话，紧贴棱边部分（也是紧靠一阶mur吸收的4条边界上）的点又该怎么处理？ =&<$I  
        本来都快编完程序了，突然想到这个问题了！钻牛角了！请教各位了！谢了！ VA*79I#_q  
       不胜感激！

[~&:`I1  
'FvhzGn9Q
 
因为FDTD的迭代式是显式迭代方法，所以“边界面上的其它四个点”的时间肯定是前一个时刻或者几个时刻（由迭代式决定）。

引用第2楼zks800于2008-09-10 21:28发表的  :
因为FDTD的迭代式是显式迭代方法，所以“边界面上的其它四个点”的时间肯定是前一个时刻或者几个时刻（由迭代式决定）。

_=M'KCL*)  
xbC-ueEj  
       以前你做过这样的程序吧，再给小弟讲解详细一点，抓住一个牛人就不放了！呵呵！ 4Xho0lO&  
       计算点周边4个点，其中一侧的点和计算点时间同步，另一侧就是上一时间步的，是不是应该这样理解？假设这样的话，紧邻用一阶mur计算点时就都是上一步了？这样讨论有点麻烦呀？ |~vQ0D  
       我用上面的程序时，发现运行时间长了吸收边界就出问题了，而刚开始没有问题！并且好像还是从一条棱边开始的发散的！明天再来仔细看看！谢了！

     自己解决了，呵呵，谢谢各位了！ 4^^rOi0  
     程序中P14取消，然后把棱边部分一部分放到前面另一部分放到后面就好了！呵呵