1178阅读
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); ;PCnEs
} \T`InBbf
else if(i>=IMAX+1-cpml_layers) P^&%T?Y6z
     { AyKvh
    kedx=1.0/(ke_x1[ii]*dx); {F6>XuS=u
    ii=ii-1; m2F2
} V(I!HT5.W
else Ebw1 %W KC
{ 29("gB
kedx=1.0/dx; tJ>%Xop
} G0oY`WXOB
} HY0q!.qog
[/td][td][/td][/tr][tr][td]诸位好，我想知道为什么别人编写CPML吸收边界的时候，在计算迭代系数的时候，为什么用ii变量来计算kedx? K+ZJSfO6
直接计算不就可以了么？比如： D99g}
ii=cpml_layers; '3TwrY?-
for(i=1;i<=IMAX-1;++i) BgRfy2:
{ `Fnl<C<
if(i<=cpml_layers) 7|%|w
{ m?'H7cFR
       kedx=1.0/(ke_x1*dx); U_i%@{
} ZKy)F-yX
else if(i>=IMAX+1-cpml_layers) q,`"Z)97
     { "O{j}QwY
    kedx=1.0/(ke_x1*dx); ^0)Mc"&{
} =r"-Pm{
else w^OV;gp
{ 3<}\{jT
kedx= .. O<7Q>m
!~Vo'ykwx'
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

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

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

发帖回复