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

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

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

简单介绍一下 je，ke 为计算空间在y和z轴上的大小，P，Q位置关系：P0为计算点，P1234 为在边界面上P0周围的点，Q01234 为紧邻边界面上的点！ Ehy(;n)\  
<n_?$ TJ  
这样编程计算，是不是部分P点计算的是上上时间步的值？请教了！ VnuG^)S  
6>?qBWW  
%     p0_ez_x=zeros(2,je,ke);              % 前后面P0处 ez  t-1步的值 lU 9o"2  
%     q00_ez_x=zeros(2,je,ke);             % 前后面Q0处 ez  t-1步的值 1^IMoC7$#  
%     q0_ez_x=zeros(2,je,ke);              % 前后面Q0处 ez  t步的值 v9T3=  
%     p1_ez_x=zeros(2,je,ke);              % 前后面P1处 ez  t步的值 8;M,l2pmR{  
%     p2_ez_x=zeros(2,je,ke);              % 前后面P2处 ez  t步的值 1dD%a91  
%     p3_ez_x=zeros(2,je,ke);              % 前后面P3处 ez  t步的值 ]W3u~T*  
%     p4_ez_x=zeros(2,je,ke);              % 前后面P4处 ez  t步的值 cg )(L;  
%     q1_ez_x=zeros(2,je,ke);              % 前后面Q1处 ez  t步的值 n3e,vP? R  
%     q2_ez_x=zeros(2,je,ke);              % 前后面Q2处 ez  t步的值 O!PGZuF  
%     q3_ez_x=zeros(2,je,ke);              % 前后面Q3处 ez  t步的值 DzC`yWstP  
%     q4_ez_x=zeros(2,je,ke);              % 前后面Q4处 ez  t步的值 lB}?ey  
    for j=3:je-2                                                              %前面 EZ 吸收边界 D(z
#)oDr  
        for k=2:kb-1 5?3v;B6  
            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+... gd[muR ~  
                (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));   8hV]t'/;  
            p0_ez_x(1,j,k)=ez(1,j,k); 4n#u?)  
            q00_ez_x(1,j,k)=p0_ez_x(1,j,k); erI&XI  
            q0_ez_x(1,j,k)=ez(2,j,k); X8}r= K~  
            p1_ez_x(1,j,k)=ez(1,j+1,k); /}u:N:HA%  
            p2_ez_x(1,j,k)=ez(1,j,k-1); X&Oo[Z  
            p3_ez_x(1,j,k)=ez(1,j-1,k); $y;w
@^  
            p4_ez_x(1,j,k)=ez(1,j,k+1); Tp;W  
            q1_ez_x(1,j,k)=ez(2,j+1,k);  H_g
]q  
            q2_ez_x(1,j,k)=ez(2,j,k-1); uNewWtUb(  
            q3_ez_x(1,j,k)=ez(2,j-1,k); eg~ Dm>Es  
            q4_ez_x(1,j,k)=ez(2,j,k+1);     mR.j8pi  
        end 5REH`-  
    end

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

nUY)LnI  
!B[Y?b:  
因为FDTD的迭代式是显式迭代方法，所以“边界面上的其它四个点”的时间肯定是前一个时刻或者几个时刻（由迭代式决定）。

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

Y5>'(A>  
,!~U5~  
       以前你做过这样的程序吧，再给小弟讲解详细一点，抓住一个牛人就不放了！呵呵！ aBF<it>  
       计算点周边4个点，其中一侧的点和计算点时间同步，另一侧就是上一时间步的，是不是应该这样理解？假设这样的话，紧邻用一阶mur计算点时就都是上一步了？这样讨论有点麻烦呀？ F7p`zf@O]  
       我用上面的程序时，发现运行时间长了吸收边界就出问题了，而刚开始没有问题！并且好像还是从一条棱边开始的发散的！明天再来仔细看看！谢了！

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