[求助]色散材料界面的处理？急急急

真空---界面---色散材料-----界面----真空，请问以上结构中的界面应该怎么进行FDTD迭代呢？我当时的处理方法是界面和色散材料内部都用drude模型下由ADE得出来的FDTD时域迭代公式，结果变量很快就变得非常大，程序跑飞了，请教各位大神应该怎 .. 8gJ"7,}-'  
TWl':}  

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

真空也用Drude模型吧，只不过omega_p=0；gamma=0；e_inf=1。这样计算可行，程序写起来也方便，就不用管界面处理了。

不过参数的突变还是会出现在介质界面上啊

这个参数突变不会造成计算的发散，如果你想更精确的，可以推导下界面的参数应该与两种参数有怎样的关系，但是这个对结果的影响不大。 2h@/Q)z  
总之，如果你这么做了，计算还是发散，那就应该是其它的问题了。

我按一楼说的进行编程了，omega_p=0；gamma=0；e_inf=1的时候确实是真空传输的形式，但是我将模拟区域中间100~120的20个网格的参数omega_p、gamma设置为银的参数，而其他地方仍为0后，当又出现发散了。模拟的是一维情形

你在现有程序的基础上，把空间步长取为原来的十分之一，看还发散不发散