[求助]UPML程序中不太理解 K 的设置，希望得到各位的点拨。

我的困惑出自以下这个三维UPML的程序，就是论坛里最常见的那个三维程序，其中由于程序比较长，我节选了我的困惑所在的那段程序。 jA<T p}$!  
%***********************************************************************   r.<
JDdj  
%     3-D FDTD code with UPML absorbing boundary conditions C *KJ7nRKx(w  
%*********************************************************************** -TS? fne)  
I7zn>^0}  
%     Program author: Keely J. Willis, Graduate Student L   si?HkJv5  
%                     UW Computational Electromagnetics Laboratory t%$@fjz
 
%                           Director: Susan C. Hagness az2Xch]  
%                     Department of Electrical and Computer Engineering N5!&~~  
%                     University of Wisconsin-Madison )$_,?*fq:  
%                     1415 Engineering Drive m9 ] D6|o5  
%                     Madison, WI 53706-1691 OB ;RHNRVP  
e "n
|jRh  
。。。。。。 DEtf(lW_  
。。。。。。             CC#C  
%   x-varying material properties `)tA YH  
delbc=upml*delta; HTR1)b  
sigmam=-log(rmax)*(orderbc+1.0)/(2.0*eta*delbc); H#Q;"r3  
sigfactor=sigmam/(delta*(delbc^orderbc)*(orderbc+1.0)); M BVOfEMj  
kmax=1; |7c`(.  
kfactor=(kmax-1.0)/delta/(orderbc+1.0)/delbc^orderbc; for i=1:upml @c]Xh:I  
*/_@a?  
% Coefficients for field components in the center of the grid cell Q7(eq0na  
CjKRP;5  
    x1=(upml-i+1)*delta; ?bI?GvSh  
    x2=(upml-i)*delta; W_}/O'l{  
    sigma=sigfactor*(x1^(orderbc+1)-x2^(orderbc+1));   '\t7jQ
 
    ki=1+kfactor*(x1^(orderbc+1)-x2^(orderbc+1)); O]ZC+]}/  
    facm=(2*epsr*epsz*ki-sigma*dt); q~O>a0f0  
    facp=(2*epsr*epsz*ki+sigma*dt); 75AslL?t  
    C5ex(i,:,:)=facp; w+=Q6]FxJ  
    C5ex(ie_tot-i+1,:,:)=facp; LS6ry,D"7  
    C6ex(i,:,:)=facm;
8t[t{"  
    C6ex(ie_tot-i+1,:,:)=facm; ox4W$YdMG  
    D1hz(i,:,:)=facm/facp; RUX8qT(Z  
    D1hz(ie_tot-i+1,:,:)=facm/facp; "@)lH  
    D2hz(i,:,:)=2.0*epsr*epsz*dt/facp; s5AgsMq  
    D2hz(ie_tot-i+1,:,:)=2.0*epsr*epsz*dt/facp; L$hc,  
    D3hy(i,:,:)=facm/facp;   {:)vwUe{  
    D3hy(ie_tot-i+1,:,:)=facm/facp; >qpqQ; bm  
    D4hy(i,:,:)=1.0/facp/mur/muz; _   lwJipIO  
    D4hy(ie_tot-i+1,:,:)=1.0/facp/mur/muz; \O? u*  
_z]v<,=3M  
。。。 .. wT/6aJoX  

N9`97;.X  

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

没有哪位朋友能给个建议或是解答么？ 顶一下，谢谢！！

如果你比较UPML和CPML。一个差别是K的值。好象K的值作用是伸展坐标（对应SIGMA的值是吸收波的能量。）K不为1的时候，可能回应想到对消逝波的吸收。但K直需要优化。 esteFLm`6  
O4N-_Kfp/  
这只是我的个人理解。一起学习哈。

谢谢你的回复，那UPML中，K一直设为1可以么，不按照talflove书中的公式，不影响结果？

taflove的那个公式？

就是第三版 talflove 书里第292页，7.6.2节中的 （7.60a,b)式的K迭代公式。如果Kmax设为1，那K不就一直为1了么？这样做可以么？

恩。你看他在书314页上给的UPML的例子，KMAX设为1。 fMgcK$  
W<Bxm|  
对于CPML，KMAX ，SIGMA，ALPHA是需要优化的参数。 0c%@e2(N  
Hf!9`R[  
恩恩。觉得taflove的书在这个问题上的讨论，单薄了些。

嗯，谢谢你的解释。我就先当它为1吧，肯定有他的道理。