1180阅读
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); XN Uw
} lL:!d.{
else if(i>=IMAX+1-cpml_layers) E7iAN\vo
     { ]-X6Cl
    kedx=1.0/(ke_x1[ii]*dx); 9&f+I@K
    ii=ii-1; A_S7z*T
} U&BCd$
else Jk(V ]
{ GxD`M2
kedx=1.0/dx; {guOAT-w
} x6R M)rr
} AH&RabH2
[/td][td][/td][/tr][tr][td]诸位好，我想知道为什么别人编写CPML吸收边界的时候，在计算迭代系数的时候，为什么用ii变量来计算kedx? mA$y$73=T
直接计算不就可以了么？比如： y@~ VE5N
ii=cpml_layers; u:|^L]{
for(i=1;i<=IMAX-1;++i) qA!p7"m|
{ Fnak:R0
if(i<=cpml_layers) stajTN*J
{ <A Hzs
       kedx=1.0/(ke_x1*dx); *w4#D:g
} ?V!5VHa
else if(i>=IMAX+1-cpml_layers) fEL 9J{
     { u~M$<|;
    kedx=1.0/(ke_x1*dx); C @<T(`o
} vO2WZ7E!
else uOzoE_i
{ <L4$f(2
kedx= .. G$TO'Ciu:
;&6 {c
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线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与所在空气边界距离变化规律一样。

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

发帖回复