24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Yaniel

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 89.7
散金: 509
红花: 2
帖子: 445
在线: 73.6小时
虫号: 740078
注册: 2009-04-04
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】【电磁场数值计算】对称边界的差分

运用有限差分法对如下图所示的区域进行数值求解，利用该结构的对称性，只需求

解右上角1/4区域的场即可。请问在对称边界上，如果差分格式中涉及1/4区域以外

的场，怎么处理？实际计算中，分TE/TM两种情况，最终方程的未知量只有Hz/Ez.

[ Last edited by woshilsh on 2009-6-13 at 22:21 ]

回复此楼

» 猜你喜欢

stronger than yesterday!

1楼 2009-06-04 15:47:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Yaniel

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 89.7
散金: 509
红花: 2
帖子: 445
在线: 73.6小时
虫号: 740078
注册: 2009-04-04
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by redcrag at 2009-6-5 09:06:
Yee网格差分格式？看着这么熟呢，呵呵。
对称边界外面的点给跟内部对应的点场值相同，反对称则值相反。用到的边界外的点的值要折算到边界内对应点上。因此在处理涉及到边界外的点那些矩阵元素时，要对相应的边界 ...

眼神果然犀利啊，反对称指什么？对于TE（计算Hz）/TM（计算Ez）都是这样处理吗？

考虑点（i，j）上的一个差分格式（介电常数用 & 表示，网格步长det表示）

jw&Ey（i，j+1/2）=-【Hz（i+1/2，j+1/2）- Hz（i-1/2，j+1/2）】/det（x）

假如（i，j+1/2）恰好处于对称边界，则（i+1/2，j+1/2）或者Hz（i-1/2，

j+1/2）处于计算区域外，如果折射到计算区域内的话，+1还是 -1？

[ Last edited by Yaniel on 2009-6-5 at 09:23 ]