1032阅读
2回复

[讨论]关于总散场的设置问题

天天向上

UID ：83276

注册：2011-09-26
登录：2012-09-01
发帖：86
等级：仿真二级

0楼发表于: 2012-03-16 19:18:57

最近在做一维薄层色散介质的反射，用到总散场的设置问题，希望大家可以一起讨论一下，
我的思路是先算出点源在场各个网格处的场值，然后再进行迭代，总散场边界的迭代公式如下，（其余场点均为一般的迭代公式）：

总散场左边界：e(q)=e(q)-dt/(epsilon0*dz)*(h(q)+H(q-1)-h(q-1));

h(q-1)=h(q-1)-dt/(mu0*dz)*(e(q)-E(q)-e(q-1));

总散场右边界：e(q)=e(q)-dt/(epsilon0*dz)*(h(q)+H(q)-h(q-1));

h(q)=h(q)-dt/(mu0*dz)*(e(q+1)-E(q)-e(q));

其中E(q)，H(q)为源场值
但现 .. <?h`
P4H%pm{-
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线m4a1ken

UID ：33675

注册：2009-05-26
登录：2014-01-26
发帖：178
等级：积极交流四级

1楼发表于: 2012-03-17 20:18:34

楼主，其实这个小问题，你自己仔细推导下公式应该可以解决的，论坛上的朋友都那么忙，一般看到这种帖子是不会静下心来看的。
感觉存在两处错误：
第一，我推导了一下，感觉电场和磁场的递推公式有问题，减号都应该是加号吧，用红色标注。
第二，在每一个公式里，处理散射场时，加减有问题，仔细想想应该知道，散射场左边界和右边界对边缘入射场的处理应该是不一样的，用绿色标注。
我推导的公式如下：
总散场左边界：
e(q)=e(q)+dt/(epsilon0*dz)*(h(q)-H(q-1)-h(q-1));
h(q-1)=h(q-1)+dt/(mu0*dz)*(e(q)-E(q)-e(q-1));

总散场右边界：
e(q)=e(q)+dt/(epsilon0*dz)*(h(q)+H(q)-h(q-1));
h(q)=h(q)+dt/(mu0*dz)*(e(q+1)+E(q)-e(q));

我在想，楼主应该不会把红色的四个符号都弄反吧，所以可能是咱们坐标建立不同导致的。如果是这样，我建议你检查总散边界处入射波的加入，左右两个边界的符号一定是不同的。

离线myp841025

天天向上

UID ：83276

注册：2011-09-26
登录：2012-09-01
发帖：86
等级：仿真二级

2楼发表于: 2012-03-18 10:37:48

谢谢你这么细心又耐心的解答！

天天向上

发帖回复