[讨论]用FDTD仿真平面波传播遇到问题

离线dcmzf

天行健君子以自强不息

UID ：58042

注册：2010-04-25
登录：2018-12-01
发帖：38
等级：仿真新人

0楼发表于: 2010-07-09 09:49:47

书上描述平面波（TEM）传播方向不存在场分量，电场只在x方向，磁场只在y方向，电场磁场同相，即同时出现最大值同时出现最小值，那么请问各位平面波的能量是怎么传送的？我用FDTD仿真结果，磁场和电场并不 ..

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

离线caocheng82

UID ：10116

注册：2008-03-28
登录：2025-05-26
发帖：697
等级：积极交流六级

1楼发表于: 2010-07-09 10:08:42

你的模型是什么样的，能够画个示意图么？
我觉得你的结果并不是TEM波

离线funnyhaha

UID ：16429

注册：2008-08-05
登录：2015-03-17
发帖：173
等级：积极交流四级

2楼发表于: 2010-07-09 12:14:42

平面波，电场和磁场0相位，应该可以计算出波印廷矢量不为０的。能量可以传输。

”FDTD仿真结果，磁场和电场并不同相“　－　你的平面波的加入有问题。

离线dcmzf

天行健君子以自强不息

UID ：58042

注册：2010-04-25
登录：2018-12-01
发帖：38
等级：仿真新人

3楼发表于: 2010-07-09 17:49:22

Ein[10]=sin(2.0*pi*f*ts* float (n));   //入射波源
    for(i=0;i<incidentend-1;i++)
        Hin=Hin-rh1*((Ein[i+1]-Ein)/dz);//磁场递推 iD])E/
    for(i=1;i<incidentend-1;i++) aot2F60J,
        Ein=Ein-rh2*((Hin-Hin[i-1])/dz);//电场递推 Fo;:GX,b
:&:P4Y1 E
    Ein[0]=Einal[1]+uz1*(Ein[1]-Einal[0]);//吸收边界 H8dS]N~[Y
    Ein[incidentend-1]=Einbl[0]+uz1*(Ein[incidentend-2]-Einbl[1]); "%a<+D
Iu`S0#+
    Einal[1]=Ein[1];//更新变量 Zby3.=.e
    Einal[0]=Ein[0]; y&O?`"Uv/M
    Einbl[0]=Ein[incidentend-2]; birc&<
    Einbl[1]=Ein[incidentend-1]; r%ebC
xVX:kDX
这是我写的实现过程，不知道是否有问题，谢谢关注

离线funnyhaha

UID ：16429

注册：2008-08-05
登录：2015-03-17
发帖：173
等级：积极交流四级

4楼发表于: 2010-07-11 06:37:31

同意一楼的。你的代码好象不是总场／散场的一维波。边界条件像是从Sullivan那里得到的。

发帖回复