1179阅读
1回复

[求助]编写CPML吸收边界过程中遇到的问题

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

0楼发表于: 2010-08-29 09:29:33

[table=100%][tr][td]
ii=cpml_layers;
for(i=1;i<=IMAX-1;++i)
{
if(i<=cpml_layers)
{
       kedx=1.0/(ke_x1*dx); ~rVKQ-+4&
} Z;|0"K
else if(i>=IMAX+1-cpml_layers) x.0k%H
     { c9Y2eetO
    kedx=1.0/(ke_x1[ii]*dx); _A@fP[C
    ii=ii-1; [u`17hyX
} aF;]7i@
else FYx `o\
{ : ;E7+m
kedx=1.0/dx; w>`h3;,2
} 1qB!RIau
} <3i4NXnL2
[/td][td][/td][/tr][tr][td]诸位好，我想知道为什么别人编写CPML吸收边界的时候，在计算迭代系数的时候，为什么用ii变量来计算kedx? lpM>}0v
直接计算不就可以了么？比如： W+F<P@[u<$
ii=cpml_layers; ]Ssw32yn
for(i=1;i<=IMAX-1;++i) rW=k%# p
{ iBGSBSeL&
if(i<=cpml_layers) *`@XKK
{ U}A|]vi@
       kedx=1.0/(ke_x1*dx); s=\LewF1<
} @%I_&!d
else if(i>=IMAX+1-cpml_layers) Q:-%3)g<<
     { *G2)@0 {
    kedx=1.0/(ke_x1*dx); `:-@E2
} reqfgNg
else RTgQ#<W8
{ N$=(1`zM=
kedx= .. NB!'u) lFD
s2(w#n)
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线yanghchuan

学习

UID ：42864

注册：2009-09-30
登录：2012-03-11
发帖：43
等级：仿真一级

1楼发表于: 2010-09-13 19:47:53

我想楼主这段代码有问题哈。kedx应该是数组哈。。。
它这里用ii是因为在x方向最右边设置fdtd叠代系数时方便而已，保证k参数在x方向上两端的pml参数的变化规率是一样的哈。
比如说设 K参数在x左边方向设为k1,k2,k3，那么在右边应该设为k3,k2,k1变化，中间为要截断的计算区域，这样就是保证
同一参数随pml与所在空气边界距离变化规律一样。

不知道有没有说清楚。。。

发帖回复