[讨论]大空间步长FDTD衍生算法

本人用FDTD算法模拟激光的传播。 >M^ 1m(  
5@Q4[+5&_  
激光的波长特别短，一般为532nm或者1064nm。由于标准FDTD算法中要求空间步长要比电磁波波长小很多，所以现在可以计算的物理空间仅仅能够达到微米等级。 MOG[cp  
K0\a+6kh  
涉及到的物理问题要求模拟的物理空间要在毫米或者里面等级（使用标准FDTD计算机资源要求很高）。 L/GM~*Xp(O  
D.<CkDB  
故在此发帖，向大家询问有没FDTD的衍生算法，使得计 .. O`dob&C  
S~F:%@,*  

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

FDTD 的衍生类算法比如High-order FDTD, MRTD，PSTD都可以使空间步长变大。 但是对时间步长的要求会加强（时间步长要变小）。这个是个矛盾。是FDTD的算法研究的一个努力方向。 kmc9P&  
gv[7h'}<  
对于你的问题，如果空间网格尺寸要求即使能提高三到五倍。可能内存的需求还是很大,  波长在微米级， 结构在毫米级----- -千倍的比率呀!  ---即使空间网格可以按照Nyquist取样定律来极限取样，就是一个波长内只取两个点。 粗略估计一毫米内至少是２０００点。三维FDTD问题如果网格数是２０００X２０００X２０００。变量使用Float(4bytes.  .那么内存需２０００＾３＊６＊４=１９２GB。 有点恐怖啊!   这就是很多激光器的模拟是２维的原因， 三维的都非常消耗内存和CPU time 。CLustering Simulation是个好的方向, 但对个人编程比较难,, 软件又比较贵。 如果结构不是非常复杂。或者有波导引导，那么建议尝试BPM或者别的算法。 l(]\[}.5  
"j a0,%3  
另外问一下， 您的类似的模拟激光传播的文献中， 多数是用什么算法?