[求助]激励为软源时波传播发生变化的问题

我在susan3D－UPML的基础上改动了一下，把他的高斯脉冲源用双指数脉冲来代替，我的这个模型是基于实际工程的，所以我用双指数函数来模拟电磁脉冲 $f<Rj/`&  
我固定Y然后令 3&JsYQu  
ez(ih_bc:ih_tot-upml,js,kh_bc:kh_tot-upml)= ez(ih_bc:ih_tot-upml,js,kh_bc:kh_tot-upml)+J0*(exp(-a*n*dt)-exp(-b*n*dt)); rru `%~'O  
式子中ih:bc:ih_tot-upml为X方向的总场网格数，kh_bc:kh_tot-upml为Z方向的总场网格数 X<"W@  
J0*(exp(-a*n*dt)-exp(-b*n*dt))为双指数脉冲函数J0，a,b都为常数 .W%{j()op  
EZ=EZ+EIN来构成一个简单的延Y轴传播的平面波，可是传播过程中波发生了变化 PfVjfrI[  
运行40步时 )Ikx0vDFQ  

图片:40.jpg

+Gvf5+ 5VR  
运行70步时  kQm\;[R  
  

图片:70.jpg

r&ex<(I{  
运行100步时 xpVYNS{c+|  
  

图片:100.jpg

/ZK
O\q  
从图中可以看到两边的电场值和中心相比有延迟 QFf lx  
我在调试过程中一开始认为是稳定性的问题，后来我按照葛德彪书中所写调整了dt和dx，效果依然不好 
i'9  
再后来我调整了PML层数设置成PML为20层，但结果依然不对。 #fYz367>  
是不是我EZ=EZ+EIN这个方法设置形成简单的平面波是有问题的呢？ PDssEb7  
或者dt不够小，不能充分反映双指数脉冲的变化，就像gwzhao师兄说 .. H\<C@OkJS}  
nZM|8  

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

Ez = Ez + Ein ? {cKKTDN  
5+'1 :Sa(i  
是因为简化才这么写的吗？ pl V]hu27K  
do%6P^qA  
平面波，线源，点源 的区别 你应该比较了解了吧？

是啊，这是为了发帖简单才这么写的，原式是 c Q|nL  
ez(ih_bc:ih_tot-upml,js,kh_bc:kh_tot-upml)= ez(ih_bc:ih_tot-upml,js,kh_bc:kh_tot-upml)+J0*(exp(-a*n*dt)-exp(-b*n*dt)); SI)u@3hl&w  
                                                                         EZ=EZ +EIN CbA!  
我固定Y然令把这个XZ面上的所有EZ=EZ+EIN，从而形成一个简单的平面波

你只是在一个平面上加了源，而不是加了平面波。

这种方法不能形成平面波吗，我一开始是按葛老师书中的方法来做，可后来我导师跟我说这样加源，可以形成简单的平面波，于是我就照做了，不过出来图就这了样  

这不是连续性的问题，也不是稳定性，步长的问题，应该是你平面波加了不对。 @xAfD{}f!  
>c@1UEwkm  
想一下，在一个 X * X的平面上，加了一个幅度相同的源，随着时间，还是会慢慢扩散出去，形成类似点源的那种图形。 Y.8mgy>   
6H!"oC&  
而平面波的传播不是这样的，你要加个连接边界，确保传播是平面波才行。

因为你加了吸收边界，所以可以近似为无限真空。 8z h{?0  
T#e ;$\  
你想一下 /a!M6:,pX  
d%:   
无限真空中，你在一个 小的区域里 加了一个源 ， 随着时间传播， 电磁场分布会怎么样，肯定不会是平面波啊。

当初认为模型的确实有限的真空区域，加连接边界是不是可以理解为有了连接边界才能形成一个有限的真空区域，保证EZ=EZ+EIN按平面波的方式传播

按照师兄这个说法，就是按时葛德彪书上的方法加入平面波，只不过我加入源的时候不要那么复杂，只要在连接边界一侧加入源EZ=EZ+EIN,然后再另一侧把源EZ=EZ+EIN减掉就可以是吧

你不是3D的吗，平面波在这个区域里面传播，总共涉及到6个面你都要处理一下吧。 Vy:I[@6@+  
你要保证在连接边界范围内的波是以平面波的方式传播出去，特别是边界上的点，这也很容易调试的，因为对一个平面上的物理量应该是相同的。