登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
求助UPML
发帖
回复
1156
阅读
5
回复
[
求助
]
求助UPML
离线
geniuszn
UID :65271
注册:
2010-08-24
登录:
2011-06-21
发帖:
36
等级:
仿真新人
0楼
发表于: 2011-02-20 11:36:39
如果计算区域中有两种介电常数不同的介质。那是不是应该将这两种介质的介电常数延伸到UPML中,针对Taflov的那本书上的说法,需要让sigma与介电常数成正比,以保证他们的比是个常数。另外,在与不同介质对应的UPML区域内,所有的计算时用到的介电常数都用所对应的工作区域的介质的介电常数?
*p.P/w@1
G!54 e
是这样吗?
PT|W{RlNl
SZCFdb
请问有没有高手写过类似的程序?
?hS n)
) tGC&l+?/
..
b "4W` A
X}yEMe{T
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
离线
m4a1ken
UID :33675
注册:
2009-05-26
登录:
2014-01-26
发帖:
178
等级:
积极交流四级
1楼
发表于: 2011-02-24 20:56:52
Taflove的那本书上第七章讲到这个问题了吧,我是参考(7.57)式的,那里讲的就是你的问题,所以Taflove用(7.57)式把前面的(7.55)修正了,这就保证了不同介质伸进UPML都能实现。你的理解和我的一样,我觉得应该是对的。
jGd{*4{3+
也就是说,按照Taflove后面的公式写程序,已经保证了,你应该不需要担心了。
Rw*l#cr=.
至于你说的‘在与不同介质对应的UPML区域内,所有的计算时用到的介电常数都用所对应的工作区域的介质的介电常数’,我想应该是这样的。你可以试试嘛,把空间设成不同的介质,然后设个点源之类的,检测下UPML的吸收效果。记得有了结果告诉我下啊。
共
条评分
离线
m4a1ken
UID :33675
注册:
2009-05-26
登录:
2014-01-26
发帖:
178
等级:
积极交流四级
2楼
发表于: 2011-02-26 10:11:43
楼主,我按照LS的说法做了程序,计算光从一侧入射到半无限介质的问题,虽然问题很简单,但是由于在这种情况下TF/SF方法引入入射波要用到三波技术,并且空间UPML的设置要以相邻的介质为准,即介质要伸到UPML里。计算结果非常稳定,UPML吸收效果很好,总场外散射场为0。证明了三波技术和介质伸进UPML都是正确的。所以,LS的说法是正确的,楼主放心用Taflove的公式吧。
G?{BVWtl}
gXZ.je)NM
其实可以用一个简单的一维的UPML验证,比如下面的图,入射光在第20个网格,横坐标>100是n=1.5,横坐标小于100是n=1.两侧的介质都伸进了UPML(UPML为10个网格),看电场振幅说明UPML吸收效果很好。很简单的一维程序,附在下面。
bBc<yaN
|rka/_
_y#t[|}w
共
条评分
离线
lw300880
UID :39822
注册:
2009-08-20
登录:
2022-03-20
发帖:
176
等级:
仿真二级
3楼
发表于: 2011-11-07 21:28:23
顶一下,观摩楼主
共
条评分
离线
xiaoyuan
UID :53338
注册:
2010-02-28
登录:
2015-12-05
发帖:
156
等级:
仿真二级
4楼
发表于: 2011-11-10 02:28:18
我觉得根本就不需要延伸,pml层里的相对介电常数就用真空的,你看pml推pml最根本的公式,理论上说只要符合那个无反射条件就可以,任何相对介电常数的值都行,只要符合那个条件。跟你主区域没必要一样。
>KmOTM<{
不过我都试了,用epsr=1效果比较好。所以就用epsr=1,无需延伸。
w$WN` =
共
条评分
离线
etooooo
UID :62978
注册:
2010-07-06
登录:
2013-02-19
发帖:
109
等级:
仿真二级
5楼
发表于: 2011-11-17 14:33:43
我怎么记得是 在与不同介质对应的UPML区域内所用到的介电常数要取一样的
共
条评分
发帖
回复