[求助]调整网格步长以后发吸收边界有反射

用的是susan的3D UPML程序，调整网格步长为1cm（原来是2mm），总场扩大到100×100×100，然后在中间放置一个微分高斯脉冲，发现有反射，很郁闷～～ us8HXvvp{  
ks %arm&  
我把PML层从10扩大到20，出图以后感觉还是有点 .. :t;i2Ck  
/1D.Ud^  

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

跟网友讨论感觉像是dx变化倒是PML的设置发生了变化 ;j[:tt\k  
delbc=upml*delta; 9'e<{mlM  
sigmam=-log(rmax)*(orderbc+1.0)/(2.0*eta*delbc); =fdW H4  
sigfactor=sigmam/(delta*(delbc^orderbc)*(orderbc+1.0)); F
Fe{=H,=  
sigma=sigfactor*(x1^(orderbc+1)-x2^(orderbc+1)); \rg;xZa5  
P_H_\KsH*(  
通过以上程序发现，dx也就是delta变大后，sigfactor变小了，导致sigma变小了，现在有反射，说明电导率不够大，没有让进去的波完全衰减掉 [d-Y1  
请教大家该如何调整

你网格长度变了5倍，从2mm到1cm，那你高斯脉冲的中心频率还是没变吗？

频率改了啊，大约为0.15G,lamda=c/f，dx＝lamda,一个波长的取样点数为20。 9>vB,8  
我用的是susan写的（allen taflove书附带程序）3D UPML程序，他那个程序，我在想是不是频率变了，那个脉冲宽度tau也要改变，波源也要作出调整 ?P#\CW  
5lM 3In@  
J(t)=J0*(t-t0)*exp(-(t-t0)^2/tau^2) aD6!x3c/  
tau=50ps,原来的频率为15GHz m5HMtoU  
t0=3*tau

是啊，要改就都改了，除非你想加入那样的源，你可以仔细研究下，做一下频谱分析。 Fta=yH}  
发散的源，频谱上看能量怎么分布的。 uaCI2I  
不发散的源，频谱怎么分布的，然后再考虑一下，为什么会这样。

把仿真频率、网格步长、时间步长绑定好，避免改了一两个参数，影响仿真的稳定性

嗯，我回头分析一下，谢谢师兄

谢谢老师的回复，我认为频率改动后，根据c/f波长就会发生变动，然后dx就改变了，根据courant稳定条件，dt=dx/2c，dt也变动，我的程序里面这些都是随着频率的变化改动的，但是出来的图，观察发现还是反射，后来根据波源 <uo@k'   
J(t)=J0*(t-t0)*exp(-(t-t0)^2/tau^2)   3k/X;:,.  
tau=50ps,t0=3*tau @R`Ao9n9V  
我把原来的时间常数tau从50ps改动为160ps反射就几乎看不出来 tK6=F63e  
但是我不太明白，这个脉冲宽度的时间常数tau到底起了个什么作用影响到我的吸收边界效果

FDTD存在数值的色散问题，高斯脉冲的宽度只有满足一定条件才能减小这个问题

数tau是脉冲宽度，决定了仿真信号的带宽，fmax=1/(2*tau) ux6p2Sk;K  
;s-@m<  
tau=50ps,原来的频率为15GHz ？？15G是仿真最大频率？ }y vH)q  
QlS_{XV  
一般PML吸收边界，像你这种仿真，5层的效果就很好了，这么还有反射？是高斯波的尾波？还是反射？