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

引用第18楼cem-uestc于2009-06-30 14:42发表的 : ,!~U5~  
损耗正切是频率的函数，也可以考虑整个仿真频率范围内当成一个常数处理 F7p`zf@O]  
损耗正切一般情况下转换成电导率sigma来处理，加入到FDTD中 KE)^S [Da  

y9:o];/  
这种方法我试了，好像不行 )P/~{Ci:T&  
@S012} xH  
因为推导出来的损耗正切转换成sigma为： lKirc2  
+w?1<Z  
sigma=tand*2*pi*frequency*er0*er. Mf:x9#  
因此其付里叶逆变换是一个复数，并且导致计算能量不收敛（强谐振）。

是一个平行板结构

说明一下 In4T`c?kQ  
楼一中的图2我是在这样的条件下得到的 20uR?/|@  
将tand=0.019在频率为0.9GHz时得到的传导率代入到FDTD中迭代得到的 i `8Y/$aT  
qoD M!~  
并用CST在设置传导率为该值时得到的结果与我用FDTD得到的一样（红色）。 j@s,5:;[  
!FnH;  
M#>f:_`<  
绿色是用CST的const tand fit模型得到的。

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

图片:未命名.jpg

z
UxF"g-W  
有哪里错误了么？

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

)Cyrs~  
这种情况下是不会发散。 jr@<-.  
考虑Debye模型的发散，可能是你代码或公式的书写错误造成的。

感谢大家特别是gwzhao和cem-uestc 的帮忙，发散的问题已经解决了， K}! VY`  
正在验证计算结果的正确性，发现使用debye的一阶模型来求解这类问题的收敛速度好慢 \nPa>2r  
@;d(>_n  
本来用const conductivity模型只需要两万步就可以收敛的，用了debye模型后两百万步长过后都还在不停的振荡，还不知道收敛没有（激励用的是一个低频分量很小的高斯脉冲） ?OvtR:hC  
.#J'+LxFr  
请问这是不是正常现象？（但是在CST中好像也没有我的这样慢呀？）不知道为什么。

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

是自己写的代码， |Xlc2?e  
const conductivity模型是用了几千步就收敛了 @w[WG:-+  
P'KaWu9z  
nUz2~z  
const tangent delta 用了debye模型 i_)j K  
是不是因为这个原因？

debye模型后两百万步长？ %qcBM~efT  
两百万步仿真，数值误差累积起来也很恐怖啊 y/K%F,WMf  
是二维还是三维FDTD？

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