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

引用第18楼cem-uestc于2009-06-30 14:42发表的 : SXn1v.6  
损耗正切是频率的函数，也可以考虑整个仿真频率范围内当成一个常数处理 _ a|zvH  
损耗正切一般情况下转换成电导率sigma来处理，加入到FDTD中 t/\J  

#N?EPV$  
这种方法我试了，好像不行 xZ}1dq8  
-aPvls   
因为推导出来的损耗正切转换成sigma为：
~r6qnC2  
Tp&03  
sigma=tand*2*pi*frequency*er0*er. o })k@-oL  
因此其付里叶逆变换是一个复数，并且导致计算能量不收敛（强谐振）。

是一个平行板结构

说明一下 ^W|B Xxo  
楼一中的图2我是在这样的条件下得到的 4viP lO  
将tand=0.019在频率为0.9GHz时得到的传导率代入到FDTD中迭代得到的 [;3` Aw  
aV?dy4o$  
并用CST在设置传导率为该值时得到的结果与我用FDTD得到的一样（红色）。 zyyt`  
dYdZt<6W<(  
eHUr!zH:  
绿色是用CST的const tand fit模型得到的。

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

图片:未命名.jpg

6 isz  
有哪里错误了么？

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

iCN@G&rVw  
这种情况下是不会发散。 LF%1)x  
考虑Debye模型的发散，可能是你代码或公式的书写错误造成的。

感谢大家特别是gwzhao和cem-uestc 的帮忙，发散的问题已经解决了， <zZAVGb4I  
正在验证计算结果的正确性，发现使用debye的一阶模型来求解这类问题的收敛速度好慢 LEhku4U.  
e+y< a~N  
本来用const conductivity模型只需要两万步就可以收敛的，用了debye模型后两百万步长过后都还在不停的振荡，还不知道收敛没有（激励用的是一个低频分量很小的高斯脉冲） *{4cc  
eiEZtu  
请问这是不是正常现象？（但是在CST中好像也没有我的这样慢呀？）不知道为什么。

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

是自己写的代码， xw-q)u  
const conductivity模型是用了几千步就收敛了 >WDpBn:  
uc<@ Fh(  
l3afuD:  
const tangent delta 用了debye模型 ZF t^q/pw  
是不是因为这个原因？

debye模型后两百万步长？ p*0Ve21i,  
两百万步仿真，数值误差累积起来也很恐怖啊 o x^lI  
是二维还是三维FDTD？

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