[求助]FDTD谐振点位置计算误差较大

引用第18楼cem-uestc于2009-06-30 14:42发表的 : ;-84
cpfu  
损耗正切是频率的函数，也可以考虑整个仿真频率范围内当成一个常数处理 %4QCUc*lr  
损耗正切一般情况下转换成电导率sigma来处理，加入到FDTD中 dLOUL9hf  

0_JbE  
这种方法我试了，好像不行 x" 2
1 Jh  
SPKen}g  
因为推导出来的损耗正切转换成sigma为： bm1+|gssn  
Ht^MY  
sigma=tand*2*pi*frequency*er0*er. W}B4^l
 
因此其付里叶逆变换是一个复数，并且导致计算能量不收敛（强谐振）。

是一个平行板结构

说明一下 'R6D+Vk/  
楼一中的图2我是在这样的条件下得到的 Bv"Fx*{W  
将tand=0.019在频率为0.9GHz时得到的传导率代入到FDTD中迭代得到的 }i_[wq{E&  
GEP YSp  
并用CST在设置传导率为该值时得到的结果与我用FDTD得到的一样（红色）。 }8WpX2U  
D[{p~x^  
h`&@>uEiq  
绿色是用CST的const tand fit模型得到的。

为什么采用以下的1st debye order 模型结果会发散？（葛老书上的） xDADJ>u2K  

图片:未命名.jpg

MiRH i<g0  
有哪里错误了么？

引用第20楼caocheng82于2009-06-30 14:46发表的  :
这种方法我试了，好像不行
因为推导出来的损耗正切转换成sigma为：
....... N1$P6ZF  

Lr~K3nb  
这种情况下是不会发散。 Eg:p_F*lr  
考虑Debye模型的发散，可能是你代码或公式的书写错误造成的。

感谢大家特别是gwzhao和cem-uestc 的帮忙，发散的问题已经解决了， z~_\onC  
正在验证计算结果的正确性，发现使用debye的一阶模型来求解这类问题的收敛速度好慢 -jy"?]ve.  
%xruPWT:k  
本来用const conductivity模型只需要两万步就可以收敛的，用了debye模型后两百万步长过后都还在不停的振荡，还不知道收敛没有（激励用的是一个低频分量很小的高斯脉冲） S ] &->5"  
;0*^98K  
请问这是不是正常现象？（但是在CST中好像也没有我的这样慢呀？）不知道为什么。

你是自己的code吗？ 要几万步才收敛？ kVnyX@  
难道你是特殊的模型吗，一般自己写的code，最多也就几千步，想看的东西都可以看到了。

是自己写的代码， j#P4Le[t  
const conductivity模型是用了几千步就收敛了 X-#mv|3  
i%PHYSJ.  
lMez!qx,=  
const tangent delta 用了debye模型 5,BkwAr+6[  
是不是因为这个原因？

debye模型后两百万步长？ KY
mWfM3^  
两百万步仿真，数值误差累积起来也很恐怖啊 C`~4q<W'  
是二维还是三维FDTD？

二维的，后面的数据的处理我都快受不了了 k1Mxsd  
特别是计算FFT变换