[求助]FDTD程序的计算精度问题

我用UPML编了一个三维的程序，刚开始用正弦波实验，在连接边界的一个面上加正弦波源，用葛老师书上（第二版）145页的加源方法，波德传播方向是z方向，电场在x方向，近场Hy的误差有7%，程序中的公式已检查没有错误，请教各位高人 .. `@Tl7I\  
.7M:AS>  

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

如果减小步长，误差会减少多少？

版主说的减小步长是不是计算的时间步少一些，还是网格尺寸再小一些？ rIJv(&l  
ywj'S7~A  
我后来又用lamda/60计算了下，误差小了些，输出边界泄露的值分别为，H最大为Hy的误差为0.0287（源磁场峰值是1），E最大为Ex的误差21.3（源电场峰值为120*pi），这些数据比原来的lamda/40时小了许多。 czu?]9;^ Z  
WO]9\"|y  
时间步数我设置的是1000步，整个计算区域为65*65*65，这样可以保证时谐场稳定，但不知道时间步数大了会不会对结果的影响也越大？

计算区域不变，又再用lamda/80计算了下，误差又小了些，输出边界泄露的值分别为，H最大为Hx的误差为0.0142（源磁场峰值是1），E最大为Ex的误差4.9740（源电场峰值为120*pi），这些数据比原来的lamda/60时小了许多。 Iy8fN"I9D  

odsLFU(  
初步结论，网格尺寸小，精度好，泄露不多，但是网格又不能太小，因为计算时间和目标大小有限制的，请教各位其他什么办法可以减小误差和泄漏？

纠正下，2楼和3楼说的计算区域是130*130*130，不是65*65*65，有正负标号，所以2倍才是

还在等待中，请大家给些建议