登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
时域有限差分法 FDTD
>
TE模,Mur边界问题是不是需要处理8个角点?
发帖
回复
1283
阅读
8
回复
TE模,Mur边界问题是不是需要处理8个角点?
离线
lwfz
UID :7619
注册:
2008-01-15
登录:
2022-08-21
发帖:
1081
等级:
积极交流六级
0楼
发表于: 2009-11-11 06:32:13
最近在编2维TE模的Mur吸收边界的程序
$C[z]}iOi
请问,是不是要对Ex,Ey各自的4个角点都需要处理?
9Dx~!(
那本书上有准确的 ..
P#O"{+`
P6gkbtg
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
离线
cem-uestc
UID :9061
注册:
2008-03-07
登录:
2019-01-05
发帖:
2575
等级:
荣誉管理员
1楼
发表于: 2009-11-11 11:21:53
是对Hz的处理啊,怎么是Ex,Ey
共
条评分
欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535
离线
lwfz
UID :7619
注册:
2008-01-15
登录:
2022-08-21
发帖:
1081
等级:
积极交流六级
2楼
发表于: 2009-11-11 11:24:40
回 1楼(cem-uestc) 的帖子
请问下版主,TE模式的话需要处理的是Hz的边界条件?
X>dQK4!R
不是说要处理边界上切向的分量吗?
8\il~IFyi
8?~>FLWTXZ
我要做的是辐射边界条件RBC,用Mur 来做
共
条评分
离线
vincentwl
UID :15476
注册:
2008-07-16
登录:
2015-06-24
发帖:
189
等级:
八级仿真大师
3楼
发表于: 2009-11-11 12:04:17
二维TE情况,E,H在网格上相差半个网格距离,角点都是Hz,为什么还要处理Ex Ey呢?
共
条评分
http://www.meta-materials.com/
离线
lwfz
UID :7619
注册:
2008-01-15
登录:
2022-08-21
发帖:
1081
等级:
积极交流六级
4楼
发表于: 2009-11-11 12:19:09
回 3楼(vincentwl) 的帖子
边界截断的话是处理Ex Ey 的边界,还是处理Hz的边界呢?
h#rziZ(
这里有点闹糊涂了........
共
条评分
离线
cem-uestc
UID :9061
注册:
2008-03-07
登录:
2019-01-05
发帖:
2575
等级:
荣誉管理员
5楼
发表于: 2009-11-11 20:58:41
是切向场分量
共
条评分
欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535
离线
lwfz
UID :7619
注册:
2008-01-15
登录:
2022-08-21
发帖:
1081
等级:
积极交流六级
6楼
发表于: 2009-11-12 10:00:24
hz(2:ie,1)=-hz2(2:ie,2)-(1/3)*(hz(2:ie,2)+hz2(2:ie,1))+...
O@-|_N*;K
(4/3)*(hz1(2:ie,2)+hz1(2:ie,1))+...
fx3oA}
(1/12)*(hz1(1:ie-1,1)-2*hz1(2:ie,1)+hz1(3:ie+1,1)+hz1(1:ie-1,2)-2*hz1(2:ie,2)+hz1(3:ie+1,2)); % Mur x=0
uoi~JF
0n-S%e5
hz(2:ie,je+1)=-hz2(2:ie,je)-(1/3)*(hz(2:ie,je)+hz2(2:je,je+1))+...
6O0CF}B*
(4/3)*(hz1(2:ie,je+1)+hz1(2:ie,je))+...
DT#F?@LG(
(1/12)*(hz1(3:ie+1,je+1)-2*hz1(2:ie,je+1)+hz1(1:ie-1,je+1)+hz1(3:ie+1,je)-2*hz1(2:ie,je)+hz1(1:ie-1,je)); %Mur x=h
e` {F7rd:
G'IRqO*]
hz(ie+1,2:je)=-hz2(ie+1,2:je)-(1/3)*(hz(ie,2:je)+hz2(ie+1,2:je))+...
5|_El/G
(4/3)*(hz1(ie+1,1:je-1)-2*hz1(ie+1,2:je)+hz1(ie+1,3:je+1))+...
Nm4 h
(1/12)*(hz1(ie+1,1:je-1)-2*hz1(ie+1,2:je)+hz1(ie+1,3:je+1)+hz1(ie,1:je-1)-2*hz1(ie,2:je)+hz1(ie,3:je+1)); %Mre y=0
q)rxv7Iu\
G A+#'R
hz(1,2:je)=-hz2(1,2:je)-(1/3)*(hz(2,2:je)+hz2(1,2:je))+...
ki=]#]rg
(4/3)*(hz1(1,2:je)+hz1(2,2:je))+...
'S#D+oF(1~
(1/12)*(hz1(1,3:je+1)-2*hz1(1,2:je)+hz1(1,1:je-1)+hz1(2,3:je+1)-2*hz1(2,2:je)+hz1(2,1:je-1)); %Mur y=h
`m+o^!SGe
%
gO%oA} !i
9L7jYy=A#
这是处理Mur 边界处用到的程序,
w)7 s]Ld
不知道什么地方出问题了,老是不行。。。。
oj%(@6L
hz2是n-2时刻的hz分量, hz1是n-1时刻的hz分量。
共
条评分
离线
lwfz
UID :7619
注册:
2008-01-15
登录:
2022-08-21
发帖:
1081
等级:
积极交流六级
7楼
发表于: 2009-11-12 10:02:11
回 5楼(cem-uestc) 的帖子
你好 请问下,处理四个边界处的hz分量,为什么运行结果总是不正确呢?
\$W\[s4I
请问下你那里有没有四个边界处hz分量的近似公式?
共
条评分
离线
zks800
UID :16797
注册:
2008-08-15
登录:
2017-02-21
发帖:
188
等级:
积极交流六级
8楼
发表于: 2009-11-12 12:25:38
角点处按照对角线的传播的方式处理的,你是否是这样做的?
共
条评分
发帖
回复