本蒟蒻初学FDTD的课程大作业,查阅了很多资料, 得到了TM模式从左侧入射实心圆柱的仿真。现将核心思想写出来,完整源代码较简易,基本思想与参考的的C语言版本较为相似。
对于物体的处理:
入射圆柱,在xy平面上处理为一个近似的圆形,用格点和中心格点的距离判断,但可能会出现其实是一个锯齿形的东西。这时可以考虑使用介电平均来处理最外一层的格点。目前考虑的是使用算出来和圆有关的格子总数,然后用方格子总面积减去圆的面积,而后将剩余面积平均到最外层格点上,即实现一个类似外面是包层里面是纤芯的圆柱。
【待填坑】