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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 对于色散介质FDTD程序的CPML吸收层问题
121/1返回列表
查看: 2193  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
关闭
请填处理意见
关闭