12

3463阅读
10回复

[讨论]用FDTD做光子晶体计算的进来看看

离线echo_8029

UID ：26115

注册：2009-02-25
登录：2012-05-24
发帖：6
等级：旁观者

0楼发表于: 2009-06-01 09:34:07

— 本帖被 gwzhao 执行加亮操作(2009-06-01) —

我需FDTD已经有一年多了。期间自己用matlab编写了计算二维，三维光子晶体的程序。入射波的引入，吸收边界条件都没有问题。现在遇到一个近场-远场外推的问题跟大家讨论讨论。
大家知道，光子晶体是周期结构，在计算光子晶体反射率时，葛德彪老师书中介绍了两种方法。一种是，用周期边界条件（但是这方法好像不适合斜入射的情况），一种是，利用高斯窗函数，对入射波进行修饰(在计算周期光子斜入射的时候，转动光子晶体本身)。
我在编写程序的时候使用的是第二种方法来计算斜入射的情况，但是在近场-远场外推的出现困难。按照葛老师书中的介绍，用投盒子方法，将近场投射到远场，在某些方向得到的远场电场或者磁场强度，比原始入射波 ..

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

离线echo_8029

UID ：26115

注册：2009-02-25
登录：2012-05-24
发帖：6
等级：旁观者

1楼发表于: 2009-06-01 09:36:42

本人做光子晶体计算也有部分心得，愿意于大家一起分享。不知道大家在做计算时是否遇到过同样的问题。大家在计算光子晶体反射率时，是用近场还是远场？

离线vincentwl

UID ：15476

注册：2008-07-16
登录：2015-06-24
发帖：189
等级：八级仿真大师

2楼发表于: 2009-06-01 23:07:15

我做metamaterial的FDTD 提取透射反射率只用到近场也在研究中多交流 QQ 100488514

http://www.meta-materials.com/

离线echo_8029

UID ：26115

注册：2009-02-25
登录：2012-05-24
发帖：6
等级：旁观者

3楼发表于: 2009-06-05 14:56:51

呵呵通过一段时间摸索，问题基本解决！望今后大家遇到问题多交流！

离线lijing

UID ：34345

注册：2009-06-03
登录：2010-01-17
发帖：13
等级：仿真新人

4楼发表于: 2009-06-06 17:22:18

二维，三维光子晶体的MATLAB程序可不可以提供给小弟学习一下？不胜感激

离线yangcong

UID ：14921

注册：2008-07-07
登录：2016-03-30
发帖：45
等级：仿真新人

5楼发表于: 2009-06-11 20:20:07

透射反射用近场

离线echo_8029

UID ：26115

注册：2009-02-25
登录：2012-05-24
发帖：6
等级：旁观者

6楼发表于: 2009-06-16 14:03:30

回 5楼(yangcong) 的帖子

透射反射用近场。但是如果我要求在某个方向上的反射率还能用近场？近场提取不出某个方面的量啊！只能用远推啊！

离线yangcong

UID ：14921

注册：2008-07-07
登录：2016-03-30
发帖：45
等级：仿真新人

7楼发表于: 2009-07-01 23:18:14

是的
要求某个方向上的得近推远
不知楼主做了没,
我只求个近场的
但TM模结果不太理想

图片:TM_TRANS.GIF
图片:TE8layer介质柱透反2.GIF

离线yangcong

UID ：14921

注册：2008-07-07
登录：2016-03-30
发帖：45
等级：仿真新人

8楼发表于: 2009-07-02 00:24:05

回 7楼(yangcong) 的帖子

回 7楼(yangcong) 的帖子

图片:平面波展开透反谱.jpg

离线sannxi

UID ：13088

注册：2008-05-28
登录：2012-10-13
发帖：70
等级：仿真一级

9楼发表于: 2010-07-30 22:29:49

回 8楼(yangcong) 的帖子

结果真好，能不能指导我一下啊？能不能告诉我你的邮箱或者QQ啊？我的邮箱为sannxi@163.com

发帖回复