5675阅读
38回复

[已解决]关于高斯制下完全匹配层的差分格式问题

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

0楼发表于: 2008-10-20 16:57:09

葛德彪老师在他的时域有限差分法书中给出了国际制麦克斯韦电磁场公式下的具体差分格式，我用他的差分格式实现了较好的吸收效果，但是当我试图运用berenger的相关思想对高斯制下的二维TE模电磁场公式作类似的磁场分裂构造时，我推出来的差分格式编程实现时发现完全没有吸收效果。我的理论功底比较差，对于电磁场的公式变换推导一直搞得不是很清楚。那位大虾曾做过关于高斯制下完全匹配 ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

1楼发表于: 2008-10-20 17:00:56

最主要的是我急于知道在高斯制下，完全匹配的阻抗匹配条件，我推导的结果是电导率和磁到率相等，但是编程运用发现效果又不好，恳请赐教！！

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

2楼发表于: 2008-10-20 18:10:42

高斯制和国际单位制下面的推导关系应该是完全一样的啊，唯一不同的就是乘（除）一个常系数而已。
你肯定是把前面的系数没弄好，把某个量放大（缩小）一个系数了。

逆流而上

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

3楼发表于: 2008-10-20 18:13:32

是啊，Maxwell方程的国际与高斯制的PML推导是一样的，可能是程序代码有问题

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

4楼发表于: 2008-10-20 19:31:06

多谢楼上两位高人的解答，只是我反复推导了好几次，都得不出正确的差分格式，功底太浅，深感惭愧！！不知道能不能发个具体的差分格式（高斯制下）的给我？？我的邮箱fz923@163.com.多谢！！

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

5楼发表于: 2008-10-21 09:33:50

我昨天晚上又仔细的重复了我的推导过程，觉得好像没有搞错，但是编程效果却不好，可能是我的定势思维影响了我的判断力，恳请高手帮我推导一下

离线wudawolf

智者不锐，慧者不傲。

UID ：7362

注册：2008-01-09
登录：2025-09-27
发帖：3709
等级：荣誉管理员

6楼发表于: 2008-10-21 10:12:04

既然已经有了国际制下的差分格式，高斯制下的的差分格式无需重新推导，只是差些常量系数而已。

在国际制下的差分格式中，将所有场量单位表示成高斯制即可。在这些单位转换过程中，你也可以同时得到那些常量系数了。

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

7楼发表于: 2008-10-21 10:23:49

本来想把具体的高斯制二维方程发上来，恳请高手帮我仿berenger思想推导一下，结果公式编辑器打出来的东西贴不上来，郁闷！！

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

8楼发表于: 2008-10-21 14:22:17

全屏拷贝，弄成图片上传上来。

逆流而上

离线fz923

UID ：18775

注册：2008-10-08
登录：2010-03-05
发帖：59
等级：仿真一级

9楼发表于: 2008-10-21 15:39:08

请帮忙看看，谢谢

发帖回复