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

to cem-uestc  师兄 l@rwf$-  
E8_Le  
按你说的Ey是切向场量，应该归到PML边界处理，那么为什么左右两测的不一样啊 RB7AI!'a?  
eybcl=zeros(iebc,je);%fields in left PML region I9un  
eybcr=zeros(ibbc,je); %fields in right PML region ITIj=!F*  
？？？？？ 8+ F}`lLA  
我当时画了个网格图对照看了看，front和back的PML层中总场跟PML交界上的Ex算到总场上的话才跟这个程序对应如下所写 Qy5Os?9"  
exbcf=zeros(iefbc,jebc); %fields in front PML region z<%g #bo  
exbcb=zeros(iefbc,jbbc);   %fields in back PML region xX/s1(P  
RefRoCD1  
说实话我还是没看懂你的写的啊～～

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

?QMclzh*-  
HS\3)Ooj>  
师兄我说得不是很严格，这个程序他没有严格分出散射场，总场跟吸收边界。我查到葛老师的书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发表的 : F;u7A]H^  
很简单，为了提高计算效率

)u;JwFstX  
师兄你好人做到底，就给我解释一下这段程序吧，我盯了老半天了 |8H_-n  
%     FRONT region G$C}?"l  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 Z?17Pu'Dp  
cbexbcf(1:iefbc,1)=0.0; h3kaD  
for j=2:jebc ~Rs#|JWB2V  
  y1=(jebc-j+1.5)*dx; %1。5和0。5是怎么来的啊，这y1,y2又什么参数 -0$:|p?@^  
  y2=(jebc-j+0.5)*dx; Edw2W8  
  sigmay=bcfactor*(y1^(orderbc+1)-y2^(orderbc+1)); 3cqQL!Gm  
  ca1=exp(-sigmay*dt/(epsz*eps(1))); # :#M{1I  
  cb1=(1.0-ca1)/(sigmay*dx); <FUon  
  caexbcf(1:iefbc,j)=ca1; 1 tPVP  
  cbexbcf(1:iefbc,j)=cb1; vt.P*Z5  
end bDDqaO ,8  
sigmay = bcfactor*(0.5*dx)^(orderbc+1); /s:w^g~  
ca1=exp(-sigmay*dt/(epsz*eps(1))); ocF>LR%P  
cb1=(1-ca1)/(sigmay*dx); p* RC  
caex(1:ie,1)=ca1; RvyuGU  
cbex(1:ie,1)=cb1; bhD ~4Rz  
caexbcl(1:iebc,1)=ca1; G}VDEC  
cbexbcl(1:iebc,1)=cb1; {)xWD%  
caexbcr(1:iebc,1)=ca1; EUy(T1Cl&&  
cbexbcr(1:iebc,1)=cb1; +Q)ULnie e  
for j=1:jebc x? N.WABr;  
  y1=(jebc-j+1)*dx; $Jp~\_X  
  y2=(jebc-j)*dx; aG(hs J)  
  sigmay=bcfactor*(y1^(orderbc+1)-y2^(orderbc+1)); ZqQ*}l5  
  sigmays=sigmay*(muz/(epsz*eps(1))); O2.'-  
  da1=exp(-sigmays*dt/muz); KY$k`f6?P  
  db1=(1-da1)/(sigmays*dx); J!hFN]M<<  
  dahzybcf(1:iefbc,j)=da1; ([R")~`(l2  
  dbhzybcf(1:iefbc,j)=db1; ,[isib3  
  caeybcf(1:ibfbc,j)=ca(1); 6YmP[%  
  cbeybcf(1:ibfbc,j)=cb(1); M+poB+K.  
  dahzxbcf(1:iefbc,j)=da(1); l4vTU=  
  dbhzxbcf(1:iefbc,j)=db(1); -#.< 12M  
end SO;N~D1Z6  
哎～～师兄啊，麻烦你给我发个邮件详细解释一吧，我实在没有勇气把疑问一个个写下去了

引用第10楼wq_463于2008-12-21 18:41发表的 : O|kOI?f  
to cem-uestc  师兄
按你说的Ey是切向场量，应该归到PML边界处理，那么为什么左右两测的不一样啊 CGbwm
Px  
eybcl=zeros(iebc,je);%fields in left PML region @? c2)0
 
eybcr=zeros(ibbc,je); %fields in right PML region 3g~'5Ao  
.......

Cbm\h/PXl  
各边界的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);可我改了以后，程序运行不对啊，希望师兄指点一下，十分感谢

j@YU|-\qh  
yE}}c{hSn  
首先有个问题要搞清楚，这个程序只是仿真包含有金属体PML层的电磁传输，不是总场/散射场的FDTD仿真

引用第12楼wq_463于2008-12-21 19:46发表的 :
师兄你好人做到底，就给我解释一下这段程序吧，我盯了老半天了 Q+[n91ey**  
%     FRONT region YtmrRDQs  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 x(1:s|Uyp{  
cbexbcf(1:iefbc,1)=0.0; 8Vr%n2M  
.......

t'n pG}`tE  
6 (]Dh;gC  
就在论坛讨论，这样的问题，许多初学者都会遇到，以后让更多人看到讨论的过程 e"|efE  
------------------------- l3)}qu  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 `sn^ysp  
因各PML层中的各计算系数的PML材料不都有值，大多为0，这样先初始化为真空中的材料值，这是一种编程上的方便

y1=(jebc-j+1.5)*dx; %1。5和0。5是怎么来的啊，这y1,y2又什么参数 n@<YI  
B-*+r`@
Bd  
这是为了提高积分精度（二阶）引入的，这在Taflove的书上有

source(n)=sin(omega*(n-delay)*dt)*exp(-((n-delay)^2/tau^2));  %?????????? _qF+tm  
.#EFLXs  
激励源为调制-高斯源 '6DBs8>1  
O(∩_∩)O哈哈~，今天挣了不少技术分，不好意思

引用第17楼cem-uestc于2008-12-21 20:20发表的 : d z|or9&  
source(n)=sin(omega*(n-delay)*dt)*exp(-((n-delay)^2/tau^2));  %??????????
激励源为调制-高斯源 28-RC>,@}  
O(∩_∩)O哈哈~，今天挣了不少技术分，不好意思

[z:
!j$K  
IYv`IS"  
9gW|}&-  
谢谢哈，看来Traflove得书我得好好看看了

引用第15楼cem-uestc于2008-12-21 20:17发表的 :
就在论坛讨论，这样的问题，许多初学者都会遇到，以后让更多人看到讨论的过程 SM#]H-3  
------------------------- lv<*7BCp  
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 I*{nP)^9  
.......

g)[V(yWu  
括号中的（1:iefbc,1）怎么理解啊。。。。。。。。 rU:
`*b<  
在网格得什么位置呢。。。