搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

光子晶体传输特性的时域精细积分法分析

杨红卫 孟珊珊 高冉冉 彭硕

引用本文:
Citation:

光子晶体传输特性的时域精细积分法分析

杨红卫, 孟珊珊, 高冉冉, 彭硕

Analysis of photonic crystal transmission properties by the precise integration time domain

Yang Hong-Wei, Meng Shan-Shan, Gao Ran-Ran, Peng Shuo
PDF
导出引用
  • 将精细积分法应用于时域有限差分法中,提出了一种求解光子晶体传输特性的时域精细积分法,并对其计算精度及稳定性进行了分析. 从一阶麦克斯韦方程出发,在空间上采用Yee元胞进行差分离散,结合吸收边界条件及激励源表达式将方程整理为标准的一阶常微分方程组形式. 通过时间步长的精细划分和指数矩阵的加法定理,在时间上利用精细积分法对齐次微分方程进行积分求解,并结合激励向量的特解得到空间离散的场分量,最终通过傅里叶变换求得方程的解. 利用时域精细积分法对光子晶体进行了实例计算,并将其结果分别与时域有限差分法和四阶龙格库塔法在精度、稳定性等方面进行了比较,结果表明时域精细积分法具有更高的计算精度,并且克服了时域有限差分法以及四阶龙格库塔法在计算稳定性上对时间步长的限制. 提出的方法具有精确、稳定的特点,为光子晶体传输特性的研究提供了一种新的有效的分析方法.
    Photonic crystals are materials patterned with a periodicity in the dielectric constant, which can create a range of forbidden frequencies called as a photonic band gap. The photonic band gap of the photonic crystal indicates its primary property, which is the basis of its application. In recent years, photonic crystals have been widely used to design optical waveguides, filters, microwave circuits and other functional devices. Therefore, the study on the transmission properties in photonic crystal is significantly important for constructing the optical devices. The finite difference time domain (FDTD) is a very useful numerical simulation technique for solving the transmission properties of the photonic crystals. However, as the FDTD method is based on the second order central difference algorithm, its accuracy is relatively low and the Courant stability condition must be satisfied when this method is used, which may restrict its application. To increase the accuracy and the stability, considerable scientific interest has been attracted to explore the schemes to improve the performance of the FDTD. The fourth order Ronge-Kutta (RK4) method has been applied to the FDTD method, which improves the accuracy and eliminates the influence of accumulation errors of the results, but the stability remains very poor if the time step is large. An effective time domain algorithm based on the high precision integration is proposed to solve the transmission properties of photonic crystals. The Yee cell differential technique is used to discretize the first order Maxwell equations in the spatial domain. Then the discretized Maxwell equations with the absorption boundary conditions and the expression of excitation source are rewritten in the standard form of the first order ordinary differential equation. According to the precise division of the time step and the additional theorem of exponential matrix, the high precision integration is used to obtain the homogeneous solution. To obtain the discretized electric and magnetic fields, the particular solution must be solved based on the excitation and then be added to the homogeneous solution. The transmission properties of photonic crystals are obtained by the Fourier transform. Practical calculation of photonic crystals is carried out by the precise integration time domain, and the accuracy and the stability are compared with those from the FDTD and the RK4 methods. The numerical results show that the precise integration time domain has a higher calculation precision and overcomes the restriction of stability conditions on the time step, which provides an effective analytical method of studying the transmission properties of photonic crystals.
      通信作者: 杨红卫, yanghongwei@bjut.edu.cn
    • 基金项目: 国家自然科学基金(批准号:11172008,11272020)资助的课题.
      Corresponding author: Yang Hong-Wei, yanghongwei@bjut.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11172008, 11272020).
    [1]

    Sun H T, Song Z X, Weng Z K, Wang D P, Jiang Y S, Yu Y 2011 Acta Photon. Sin. 40 1

    [2]

    Su A 2009 M. S. Dissertation (Nanning: Guangxi University) (in Chinese) [苏安 2009 硕士学位论文 (南宁: 广西大学)]

    [3]

    Lin M 2010 Ph. D. Dissertation (Chengdu: University of Electronic Science and Technology of China) (in Chinese) [林密 2010 博士学位论文 (成都: 电子科技大学)]

    [4]

    Wang Q 2009 Ph. D. Dissertation (Nanjing: Southeast University) (in Chinese) [王琼 2009 博士学位论文 (南京: 东南大学)]

    [5]

    Liu P 2012 M. S. Dissertation (Tianjin: Tianjin University of Technology) (in Chinese) [刘佩 2012 硕士学位论文 (天津: 天津理工大学)]

    [6]

    Ding T 2011 Ph. D. Dissertation (Beijing: University of Chinese Academy of Sciences) (in Chinese) [丁涛 2011 博士学位论文 (北京: 中国科学院大学)]

    [7]

    Zhao N S, Guan J M 2014 Laser Optoelectron. Prog. 51 042302 (in Chinese) [赵年顺, 官骏鸣 2014 激光与光电子学进展 51 042302]

    [8]

    Chen H M, Meng Q 2011 Acta Phys. Sin. 60 014202 (in Chinese) [陈鹤鸣, 孟晴 2011 物理学报 60 014202]

    [9]

    Pan W, Yu H J, Zhang X G, Xi L X 2012 Acta Phys. Sin. 61 034209 (in Chinese) [潘伟, 余和军, 张晓光, 席丽霞 2012 物理学报 61 034209]

    [10]

    Cheng J, Wang W Y, Xiong Y B, Tan W J 2012 Electro-Opt. Technol. Appl. 27 34 (in Chinese) [承军, 王玮钰, 熊耀兵, 谭文疆 2012 光电技术应用 27 34]

    [11]

    Zhang H J 2007 J. Ankang Univ. 19 74 (in Chinese) [张洪江 2007 安康学院学报 19 74]

    [12]

    Song Q, Gao J S, Wang X Y, Wang T T, Chen H, Zheng X M, Shen Z F, Ling W 2006 Opt. Instrum. 28 37 (in Chinese) [宋琦, 高劲松, 王笑夷, 王彤彤, 陈红, 郑宣鸣, 申振峰, 凌伟 2006 光学仪器 28 37]

    [13]

    Namiki T 1999 IEEE Trans. Microw. Theory Tech. 47 2003

    [14]

    Young J L, Gaitonde D, Shang J S 1997 IEEE Trans. Antennas Propag. 45 1573

    [15]

    Zhong W X 2002 Dual System in Applied Mechanics (1st Edition) (Beijing: Science Press) pp4-10 (in Chinese) [钟万勰 2002 应用力学对偶体系 (1 版) (北京: 科学出版社) 第4-10页]

    [16]

    Liu P J, Gu L C, Ren J H 2008 J. Anhui Univ. Nat. Sci. Ed. 32 61 (in Chinese) [刘沛津, 谷立臣, 任继红 2008 安徽大学学报(自然科学版) 32 61]

    [17]

    Ge D B, Yan Y B 2011 Finite Difference Time Domain Method for Electromagnetic Waves (3rd Edition) (Xi'an: Xidian University Press) pp9-21 (in Chinese) [葛德彪, 闫玉波 2011 电磁波时域有限差分方法 (3版) (西安: 西安电子科技大学出版社) 第9-21页]

    [18]

    Chen J, Lee J H, Liu Q H 2009 IEEE Trans. Antennas Propag. 57 3223

    [19]

    Liu Q H 1997 Microw. Opt. Technol. Lett. 14 134

    [20]

    Lin J, Shen W, Williams F W 1995 Comput. Struct. 56 113

    [21]

    Liu W S, Li J L 1981 Sampling Technology Principle and Application (1st Ed.) (Beijing: Science Press) pp10-360 (in Chinese) [刘文生, 李锦林 1981 取样技术原理与应用 (第1版) (北京: 科学出版社) 第10-360页]

  • [1]

    Sun H T, Song Z X, Weng Z K, Wang D P, Jiang Y S, Yu Y 2011 Acta Photon. Sin. 40 1

    [2]

    Su A 2009 M. S. Dissertation (Nanning: Guangxi University) (in Chinese) [苏安 2009 硕士学位论文 (南宁: 广西大学)]

    [3]

    Lin M 2010 Ph. D. Dissertation (Chengdu: University of Electronic Science and Technology of China) (in Chinese) [林密 2010 博士学位论文 (成都: 电子科技大学)]

    [4]

    Wang Q 2009 Ph. D. Dissertation (Nanjing: Southeast University) (in Chinese) [王琼 2009 博士学位论文 (南京: 东南大学)]

    [5]

    Liu P 2012 M. S. Dissertation (Tianjin: Tianjin University of Technology) (in Chinese) [刘佩 2012 硕士学位论文 (天津: 天津理工大学)]

    [6]

    Ding T 2011 Ph. D. Dissertation (Beijing: University of Chinese Academy of Sciences) (in Chinese) [丁涛 2011 博士学位论文 (北京: 中国科学院大学)]

    [7]

    Zhao N S, Guan J M 2014 Laser Optoelectron. Prog. 51 042302 (in Chinese) [赵年顺, 官骏鸣 2014 激光与光电子学进展 51 042302]

    [8]

    Chen H M, Meng Q 2011 Acta Phys. Sin. 60 014202 (in Chinese) [陈鹤鸣, 孟晴 2011 物理学报 60 014202]

    [9]

    Pan W, Yu H J, Zhang X G, Xi L X 2012 Acta Phys. Sin. 61 034209 (in Chinese) [潘伟, 余和军, 张晓光, 席丽霞 2012 物理学报 61 034209]

    [10]

    Cheng J, Wang W Y, Xiong Y B, Tan W J 2012 Electro-Opt. Technol. Appl. 27 34 (in Chinese) [承军, 王玮钰, 熊耀兵, 谭文疆 2012 光电技术应用 27 34]

    [11]

    Zhang H J 2007 J. Ankang Univ. 19 74 (in Chinese) [张洪江 2007 安康学院学报 19 74]

    [12]

    Song Q, Gao J S, Wang X Y, Wang T T, Chen H, Zheng X M, Shen Z F, Ling W 2006 Opt. Instrum. 28 37 (in Chinese) [宋琦, 高劲松, 王笑夷, 王彤彤, 陈红, 郑宣鸣, 申振峰, 凌伟 2006 光学仪器 28 37]

    [13]

    Namiki T 1999 IEEE Trans. Microw. Theory Tech. 47 2003

    [14]

    Young J L, Gaitonde D, Shang J S 1997 IEEE Trans. Antennas Propag. 45 1573

    [15]

    Zhong W X 2002 Dual System in Applied Mechanics (1st Edition) (Beijing: Science Press) pp4-10 (in Chinese) [钟万勰 2002 应用力学对偶体系 (1 版) (北京: 科学出版社) 第4-10页]

    [16]

    Liu P J, Gu L C, Ren J H 2008 J. Anhui Univ. Nat. Sci. Ed. 32 61 (in Chinese) [刘沛津, 谷立臣, 任继红 2008 安徽大学学报(自然科学版) 32 61]

    [17]

    Ge D B, Yan Y B 2011 Finite Difference Time Domain Method for Electromagnetic Waves (3rd Edition) (Xi'an: Xidian University Press) pp9-21 (in Chinese) [葛德彪, 闫玉波 2011 电磁波时域有限差分方法 (3版) (西安: 西安电子科技大学出版社) 第9-21页]

    [18]

    Chen J, Lee J H, Liu Q H 2009 IEEE Trans. Antennas Propag. 57 3223

    [19]

    Liu Q H 1997 Microw. Opt. Technol. Lett. 14 134

    [20]

    Lin J, Shen W, Williams F W 1995 Comput. Struct. 56 113

    [21]

    Liu W S, Li J L 1981 Sampling Technology Principle and Application (1st Ed.) (Beijing: Science Press) pp10-360 (in Chinese) [刘文生, 李锦林 1981 取样技术原理与应用 (第1版) (北京: 科学出版社) 第10-360页]

  • [1] 柯航, 李培丽, 施伟华. 基于下山单纯形算法逆向设计二维光子晶体波导型1×5分束器. 物理学报, 2022, 71(14): 144204. doi: 10.7498/aps.71.20220328
    [2] 兰伟霞, 顾嘉陆, 高晓辉, 廖英杰, 钟宋义, 张卫东, 彭艳, 孙钰, 魏斌. 基于光子晶体的有机太阳能电池研究进展. 物理学报, 2021, 70(12): 128804. doi: 10.7498/aps.70.20201805
    [3] 左依凡, 李培丽, 栾开智, 王磊. 基于自准直效应的光子晶体异质结偏振分束器. 物理学报, 2018, 67(3): 034204. doi: 10.7498/aps.67.20171815
    [4] 赵绚, 刘晨, 马会丽, 冯帅. 基于波导间能量耦合效应的光子晶体频段选择与能量分束器. 物理学报, 2017, 66(11): 114208. doi: 10.7498/aps.66.114208
    [5] 周雯, 陈鹤鸣. 基于磁光效应的二维三角晶格光子晶体模分复用器. 物理学报, 2015, 64(6): 064210. doi: 10.7498/aps.64.064210
    [6] 赵浩, 沈义峰, 张中杰. 光子晶体中基于有效折射率接近零的光束准直出射. 物理学报, 2014, 63(17): 174204. doi: 10.7498/aps.63.174204
    [7] 刘广东, 张开银. 二维电磁逆散射问题的时域高斯-牛顿反演算法. 物理学报, 2014, 63(3): 034102. doi: 10.7498/aps.63.034102
    [8] 高晖, 孔凡敏, 李康, 陈新莲, 丁庆安, 孙静. 双层光子晶体氮化镓蓝光发光二极管结构优化的研究. 物理学报, 2012, 61(12): 127807. doi: 10.7498/aps.61.127807
    [9] 刘发, 徐晨, 赵振波, 周康, 解意洋, 毛明明, 魏思民, 曹田, 沈光地. 氧化孔形状对光子晶体垂直腔面发射激光器模式的影响. 物理学报, 2012, 61(5): 054203. doi: 10.7498/aps.61.054203
    [10] 岳庆炀, 孔凡敏, 李康, 赵佳. 基于缺陷光子晶体结构的GaN基发光二极管光提取效率的有关研究. 物理学报, 2012, 61(20): 208502. doi: 10.7498/aps.61.208502
    [11] 颛孙旭, 马西奎. 一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法. 物理学报, 2012, 61(11): 110206. doi: 10.7498/aps.61.110206
    [12] 刘广东, 张业荣. 二维有耗色散介质的时域逆散射方法. 物理学报, 2010, 59(10): 6969-6979. doi: 10.7498/aps.59.6969
    [13] 魏兵, 董宇航, 王飞, 李存志. 基于移位算子时域有限差分的色散薄层节点修正算法. 物理学报, 2010, 59(4): 2443-2450. doi: 10.7498/aps.59.2443
    [14] 刘漾, 巩华荣, 魏彦玉, 宫玉彬, 王文祥, 廖复疆. 有效抑制光子晶体加载矩形谐振腔中模式竞争的方法. 物理学报, 2009, 58(11): 7845-7851. doi: 10.7498/aps.58.7845
    [15] 张玉强, 葛德彪. 一种基于数字信号处理技术的改进通用色散介质移位算子时域有限差分方法. 物理学报, 2009, 58(12): 8243-8248. doi: 10.7498/aps.58.8243
    [16] 殷海荣, 宫玉彬, 魏彦玉, 岳玲娜, 路志刚, 巩华荣, 黄民智, 王文祥. 有限开敞介质光子晶体的模式及其带结构分析. 物理学报, 2008, 57(6): 3562-3570. doi: 10.7498/aps.57.3562
    [17] 姜彦南, 葛德彪. 层状介质时域有限差分方法斜入射平面波引入新方式. 物理学报, 2008, 57(10): 6283-6289. doi: 10.7498/aps.57.6283
    [18] 刘大刚, 周 俊, 刘盛纲. 用时域有限差分法实现金属支撑杆的计算机模拟. 物理学报, 2007, 56(12): 6924-6930. doi: 10.7498/aps.56.6924
    [19] 于天宝, 刘念华. 光脉冲通过含有色散与增益型缺陷的一维光子晶体的传播. 物理学报, 2004, 53(9): 3049-3053. doi: 10.7498/aps.53.3049
    [20] 王 刚, 温激鸿, 韩小云, 赵宏刚. 二维声子晶体带隙计算中的时域有限差分方法. 物理学报, 2003, 52(8): 1943-1947. doi: 10.7498/aps.52.1943
计量
  • 文章访问数:  5000
  • PDF下载量:  205
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-11-10
  • 修回日期:  2017-01-16
  • 刊出日期:  2017-04-05

/

返回文章
返回