关于 计算精度的问题？

各位，我是用FDTD计算弹性波散射的 Z3S\@_/;  
散射体是圆柱形的空气柱，基体是钢 1z8.wdWJ}  
wv1?v_4  
dt,dx,时间步长，空间步长均满足收敛条件 )>=`[$D1t  
zis-}K<  
但是在计算过程中，场值会越来越大，最后出现 SBYRN##n_  
NaN Vv zd>yII  
费了很长时间一直找不出原因， <q_H 3|  
请教各位大侠指点迷津 .. SK-|O9Ki  
^6Xio6W  

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

FDTD就是这样，编程稍有不慎，就容易发散。

计算发散啊，程序本身有问题 Z;lE-`Z*(F  
不是精度问题

我觉得可能是计算模型的问题吧，比如二维正常FDTD公式中CA=(1-sigma*dt/2/epsilon)/(1+sigma*dt/2/epsilon),,,当遇到金属时，介电常数为负，CA的绝对值要大于1，就会出现发散现象，我也遇到这样的问题，我的介质是银，用Debye色散卷积模型计算还是发散，希望得到大家的指点！

但是，因为空气的密度是1.29很小， tWL9>7]G  
如果将空气密度数值增大到300之后，就没有发散问题了 3p$ZHH.UP  
并且如果将空气换成其他固体介质，程序就不会有发散问题了并且结果也很合理， >TwOL  
这是程序问题吗？ >aN@)=h}  
>8 VfijK  

重新检查介质分界面处场量的边界条件试试。
`hI1  
jAN(r>zVL  
或者看看葛老师书的分界面处介质参数的处理的那章试试。

恩，好的，谢谢