应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 对于色散介质FDTD程序的CPML吸收层问题
121/1返回列表
查看: 2141  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kevin123581

金虫 (正式写手)

[交流] 对于色散介质FDTD程序的CPML吸收层问题

记得看书和问大家都是说CPML的通用性非常的好，可以在普通介质（电介质）和色散介质（比如金属）通用，不需要改变的。
我用把电介质的fdtd程序改成了色散介质drude模型的fdtd，但是cpml的部分完全没改动，大家都说是可以的，但是我做出来的结果却有点问题。就是在运行时间变长的时候，会从真空介质的分界面还有和cpml的交界处开始发散。
确切上来说不是发散，只是像产生了一个新的源一样，但是在这之前的结果倒都是对的。
我试着把吸收层改成了0，也就是没有cpml，这个时候程序可以正常运行，没有问题。（当然会有由于没有吸收层导致的全反射，但是不会出现上面说的那个问题）
有朋友遇到过这个情况么？  我看taflove的书上说是在色散介质的cpml里面有个partial update，但是他没例子，我也不能理解，也就没有改动，加上大家都说没问题，我也就没考虑了
现在出现的这个问题，所有又来问问，谢谢了
对于色散介质FDTD程序的CPML吸收层问题
064.png


对于色散介质FDTD程序的CPML吸收层问题-1
067.png


对于色散介质FDTD程序的CPML吸收层问题-2
069.png


对于色散介质FDTD程序的CPML吸收层问题-3
071.png


对于色散介质FDTD程序的CPML吸收层问题-4
073.png


对于色散介质FDTD程序的CPML吸收层问题-5
075.png
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2013-05-20 18:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
对了，上面的图都是包含了cpml的，可以看得到左右两边明显的cpml
赞一下 回复此楼
2楼2013-05-20 18:15:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinayang5

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我也想知道要为什么。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
3楼2013-05-21 07:34:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
引用回帖:
3楼: Originally posted by chinayang5 at 2013-05-21 07:34:35
我也想知道要为什么。

你也有这个问题？
赞一下 回复此楼
5楼2013-05-21 08:59:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
经过测试，有以下结果： 全部真空-不发散； 全部金属-发散；cpml设定为0层，不发散；    可见发散是从金属和cpml的分界面开始的

增加cpml程序里面的指数m从8到20； 发散后延但是cpml反射增加，不可取
增加cpml的层数，由10个网格到20个网格，发散后延，但是效果不高
赞一下 回复此楼
6楼2013-05-21 12:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinayang5

铁杆木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dbb627: 金币+1, 感谢参与讨论 2013-06-07 19:44:23
我以前遇到过类似情况，一直没解决。现在换换方向了。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
7楼2013-05-21 16:33:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
引用回帖:
7楼: Originally posted by chinayang5 at 2013-05-21 16:33:16
我以前遇到过类似情况，一直没解决。现在换换方向了。

能详细的说下你当时遇到的情况么   

现在的推断就是cpml和色散介质接触面的问题
赞一下 回复此楼
8楼2013-05-22 09:29:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

waocei

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
请教一下各位：，在fdtd计算中为什么要看稳态磁场Hz分布，电场分布及电场【E】^2分布？这些有什么区别吗？在这方面我是新手，麻烦各位解答，不甚感激！！
赞一下 回复此楼
多沟通多交流多坚持
9楼2013-05-23 21:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asdry

木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
祝楼主诸事顺利，心想事成！
赞一下 回复此楼
10楼2013-05-23 22:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
引用回帖:
9楼: Originally posted by waocei at 2013-05-23 21:39:56
请教一下各位：，在fdtd计算中为什么要看稳态磁场Hz分布，电场分布及电场【E】^2分布？这些有什么区别吗？在这方面我是新手，麻烦各位解答，不甚感激！！

hz的我没看过啊，电场分布这个不用说了吧，这个不就是你FDTD要算的东西么。 至于平方，那个是intensity分布，也就是能量
赞一下 回复此楼
11楼2013-05-24 10:18:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kevin123581

金虫 (正式写手)
还发现一个新的问题， 就是所有条件不变，只改变计算区域的大小，比如从10um增大到15um,  取一个面计算波印廷矢量的结果居然不一样，真空中是一样的结果，但是当加入介质后结果就不一样了
赞一下 回复此楼
12楼2013-06-06 21:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
wangyile4楼
2013-05-21 08:47   回复  
相关版块跳转 我要订阅楼主 kevin123581 的主题更新
121/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭