登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
透射系数大于1的老问题,应该有前辈解决过
发帖
回复
8122
阅读
6
回复
[
已解决
]
透射系数大于1的老问题,应该有前辈解决过
离线
wei522538
UID :77280
注册:
2011-05-12
登录:
2013-01-14
发帖:
62
等级:
仿真一级
0楼
发表于: 2012-04-06 09:31:34
关键词:
仿真
频率
激励
FDTD
文献
本人采用FDTD在算声子晶体的透射系数时,一些频率的透射系数大于1,令人很费解。具体模型是这样的,X方向采用PML吸收法,Y方向采用周期边界法。X方向左边一端给的是线激励源,右边一端在中间处(一个点)收集透射过来的信号。声子晶体X方向是五个点阵。但即使收集信号已经完全收敛了,即最后浮动振幅趋近于零,透射系数中仍有一些频率大于1。和文献中的图形对照(也是文献作者仿真的图形),大致通带和阻带都能对应的上,就是不知别人是怎样让频率都小于1 的。请前辈或已经解决这个 ..
hq.z:D
Bv$;yR
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
离线
wei522538
UID :77280
注册:
2011-05-12
登录:
2013-01-14
发帖:
62
等级:
仿真一级
1楼
发表于: 2012-04-06 09:48:15
这是求得的透射系数图形
.AH#D}m
图片:未命名.jpg
共
条评分
离线
m4a1ken
UID :33675
注册:
2009-05-26
登录:
2014-01-26
发帖:
178
等级:
积极交流四级
2楼
发表于: 2012-04-06 14:08:22
我也做过类似的问题,确实有时候出现透射大于1的情况,但是我觉得这种情况可能由很多原因造成,比如PML设置的不够好(对消逝波吸收不好),或PML厚度不够半个波长,感觉这类问题主要是吸收边界造成的。
;E l"dqH
但是,看了楼主的情况,我觉得问题应该这样解决:既然你是线激励,那么右端检测也应该在一条平行于y方向的线做检测(而不是一个点),即把线上的所有点都检测。一般透射率用能流表征,你应该求解沿x方向的能流Sx(电磁波用能流,不知道声波用什么),然后把这条线上的能流做积分。也就是出射端沿平行于y的一条线把能流积分,作为透射。
v}Gq.(b
如果我没猜错,你计算这种情况要分两步走,第一步是没有结构,在你所谓的出射点监测(作为入射波),第二次加结构,同样在你的出射点监测(作为透射波),然后fft后作比值得到透射谱。我建议你不用一个点,要用整条线,这样才能保证能流守恒,也就是透射小于1。如果仅仅有一个点的话,出现透射大于1是正常现象,因为很有可能你选的点透射大了,其它点小了,能流整体(在沿y方向的一个周期内)是守恒的。也正是这样的原因,你可以算出通带和阻带,只是透射出现了问题。
共
2
条评分
hefang
技术分
+4
积极参与讨论+技术分 论坛感谢您的参与
2012-04-06
hefang
rf币
+40
积极参与讨论+技术分 论坛感谢您的参与
2012-04-06
离线
wei522538
UID :77280
注册:
2011-05-12
登录:
2013-01-14
发帖:
62
等级:
仿真一级
3楼
发表于: 2012-04-06 14:33:03
谢谢楼上的建议,计算透射率确实是按照你猜想的方法。非常感谢你的建议!
共
条评分
离线
yangcong
UID :14921
注册:
2008-07-07
登录:
2016-03-30
发帖:
45
等级:
仿真新人
4楼
发表于: 2012-04-08 18:30:22
线激励怎么理解,是平面波么?若不是周期边界难实现。
共
条评分
离线
m4a1ken
UID :33675
注册:
2009-05-26
登录:
2014-01-26
发帖:
178
等级:
积极交流四级
5楼
发表于: 2012-04-09 08:24:38
回 4楼(yangcong) 的帖子
线激励就是平面波,不用周期边界就只能靠总场/散射场引入了。
共
条评分
离线
zxhuang
UID :15286
注册:
2008-07-14
登录:
2020-02-08
发帖:
75
等级:
仿真二级
6楼
发表于: 2012-06-29 09:46:02
是你的迭代时间不够长
共
条评分
发帖
回复