[求助]新手请教Traflove书中附带的2D程序

to cem-uestc  师兄 0k0c   
&9w%n  
按你说的Ey是切向场量，应该归到PML边界处理，那么为什么左右两测的不一样啊 <`j[;>O  
eybcl=zeros(iebc,je);%fields in left PML region _JJKbi  
eybcr=zeros(ibbc,je); %fields in right PML region dIa(</ }  
？？？？？ Jo aDX ,  
我当时画了个网格图对照看了看，front和back的PML层中总场跟PML交界上的Ex算到总场上的话才跟这个程序对应如下所写 pUMB)(<k  
exbcf=zeros(iefbc,jebc); %fields in front PML region =#2qX>?  
exbcb=zeros(iefbc,jbbc);   %fields in back PML region R|J>8AL}BY  
m2q;^o:J  
说实话我还是没看懂你的写的啊～～

引用第6楼cem-uestc于2008-12-21 14:03发表的 :
从假设没有PML的一般FDTD仿真，边界上的切向场量是放在边界处理，如加上PML层，实际边界上的切向场量当成PML边界处理，就顺利成章。

<% 7P  
5}MlZp  
师兄我说得不是很严格，这个程序他没有严格分出散射场，总场跟吸收边界。我查到葛老师的书78页下面是这么说的：由于在PML内界面处sigx=0,所以对于界面上的Ex和Ey节点可以用散射场区的迭代公式，无需对其作单独处理。所以我认为总场与PML交界面上的Ey应该是放到总场里面算的，而不是放到PML里算，那样的话eybcr=zeros(ibbc,je);这个式子我认为是不对的应该是eybcr=zeros(iebc,je);可我改了以后，程序运行不对啊，希望师兄指点一下，十分感谢

引用第9楼cem-uestc于2008-12-21 17:04发表的 : nped  
很简单，为了提高计算效率

>ysriPnQ  
师兄你好人做到底，就给我解释一下这段程序吧，我盯了老半天了 U^B"|lc:[  
%     FRONT region p9/bzT34.  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 jVX._bEGX  
cbexbcf(1:iefbc,1)=0.0; W=~id"XtJ  
for j=2:jebc uP
FRh~ (b  
  y1=(jebc-j+1.5)*dx; %1。5和0。5是怎么来的啊，这y1,y2又什么参数 w|&,I4["  
  y2=(jebc-j+0.5)*dx; f8^"E $"  
  sigmay=bcfactor*(y1^(orderbc+1)-y2^(orderbc+1)); vAi"$e  
  ca1=exp(-sigmay*dt/(epsz*eps(1))); La\Q'0  
  cb1=(1.0-ca1)/(sigmay*dx); 1T!cc%ah  
  caexbcf(1:iefbc,j)=ca1; {VBR/M(q  
  cbexbcf(1:iefbc,j)=cb1; kXigX
-  
end lx
m*;?j`W  
sigmay = bcfactor*(0.5*dx)^(orderbc+1); 63.( j P1;  
ca1=exp(-sigmay*dt/(epsz*eps(1))); ; Zh9^0  
cb1=(1-ca1)/(sigmay*dx); *]Fg
fttES  
caex(1:ie,1)=ca1; HpW"lYW4  
cbex(1:ie,1)=cb1; A)OdQFet(  
caexbcl(1:iebc,1)=ca1; ?x:m;z/  
cbexbcl(1:iebc,1)=cb1; D._{E*vg  
caexbcr(1:iebc,1)=ca1; xy2\'kS`G  
cbexbcr(1:iebc,1)=cb1; !K!)S^^Po?  
for j=1:jebc l<$rqz3D  
  y1=(jebc-j+1)*dx; \^yXc*C  
  y2=(jebc-j)*dx; d "2wO[  
  sigmay=bcfactor*(y1^(orderbc+1)-y2^(orderbc+1)); IS-}:~Pi  
  sigmays=sigmay*(muz/(epsz*eps(1))); !{$qMhT  
  da1=exp(-sigmays*dt/muz); \.5F](:  
  db1=(1-da1)/(sigmays*dx); XxhsPFv  
  dahzybcf(1:iefbc,j)=da1; s5Pq$<  
  dbhzybcf(1:iefbc,j)=db1; ]t*33  
  caeybcf(1:ibfbc,j)=ca(1); dI!8S  
  cbeybcf(1:ibfbc,j)=cb(1); T0g0jr{  
  dahzxbcf(1:iefbc,j)=da(1); /d[Mss  
  dbhzxbcf(1:iefbc,j)=db(1); R'Sa?6xS4  
end OSU{8.  
哎～～师兄啊，麻烦你给我发个邮件详细解释一吧，我实在没有勇气把疑问一个个写下去了

引用第10楼wq_463于2008-12-21 18:41发表的 : W+K.r?G<j  
to cem-uestc  师兄
按你说的Ey是切向场量，应该归到PML边界处理，那么为什么左右两测的不一样啊 07FT)Q
TE  
eybcl=zeros(iebc,je);%fields in left PML region eSn$k:\W  
eybcr=zeros(ibbc,je); %fields in right PML region drI\iae{^  
.......

w763zi{  
各边界的PML层数可以不一样，但为了编程方便，一般设置一样

引用第11楼wq_463于2008-12-21 19:28发表的 :
师兄我说得不是很严格，这个程序他没有严格分出散射场，总场跟吸收边界。我查到葛老师的书78页下面是这么说的：由于在PML内界面处sigx=0,所以对于界面上的Ex和Ey节点可以用散射场区的迭代公式，无需对其作单独处理。所以我认为总场与PML交界面上的Ey应该是放到总场里面算的，而不是放到PML里算，那样的话eybcr=zeros(ibbc,je);这个式子我认为是不对的应该是eybcr=zeros(iebc,je);可我改了以后，程序运行不对啊，希望师兄指点一下，十分感谢

12OlrU  
mNDuwDd$S  
首先有个问题要搞清楚，这个程序只是仿真包含有金属体PML层的电磁传输，不是总场/散射场的FDTD仿真

引用第12楼wq_463于2008-12-21 19:46发表的 :
师兄你好人做到底，就给我解释一下这段程序吧，我盯了老半天了 (9
'be\  
%     FRONT region uO"8aD`W  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 3#mE( `|P  
cbexbcf(1:iefbc,1)=0.0; \(bj(an
y  
.......

eJaUmK:  
dEET}s\  
就在论坛讨论，这样的问题，许多初学者都会遇到，以后让更多人看到讨论的过程 R@$+t:}  
------------------------- k=|K|  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 JV%nH!Fs  
因各PML层中的各计算系数的PML材料不都有值，大多为0，这样先初始化为真空中的材料值，这是一种编程上的方便

y1=(jebc-j+1.5)*dx; %1。5和0。5是怎么来的啊，这y1,y2又什么参数 OelU D/[$  
'H!V54 \j  
这是为了提高积分精度（二阶）引入的，这在Taflove的书上有

source(n)=sin(omega*(n-delay)*dt)*exp(-((n-delay)^2/tau^2));  %?????????? GahaZ F  
s;anP0-O  
激励源为调制-高斯源 #,t2*tM  
O(∩_∩)O哈哈~，今天挣了不少技术分，不好意思

引用第17楼cem-uestc于2008-12-21 20:20发表的 : +G<}JJ'V  
source(n)=sin(omega*(n-delay)*dt)*exp(-((n-delay)^2/tau^2));  %??????????
激励源为调制-高斯源 8w03{H 0  
O(∩_∩)O哈哈~，今天挣了不少技术分，不好意思

u wH)$Pl  
X'YfjbGo  
mYntU^4f  
谢谢哈，看来Traflove得书我得好好看看了

引用第15楼cem-uestc于2008-12-21 20:17发表的 :
就在论坛讨论，这样的问题，许多初学者都会遇到，以后让更多人看到讨论的过程 PTuCN  
------------------------- ({d,oU$>y  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 U<ku_(2"#  
.......

L337/8fh  
括号中的（1:iefbc,1）怎么理解啊。。。。。。。。 wf~5lpI[  
在网格得什么位置呢。。。