Title

解電磁模態問題的譜方法新模型

Translated Titles

New Pseudospectral Penalty Models for Electromagnetic Mode Problems

DOI

10.6342/NTU.2008.01638

Authors

江舜凡

Key Words

譜方法 ; 波導 ; 光子晶體 ; 模態 ; Pseudospectral method ; Waveguide ; Photonic Crystal ; Mode

PublicationName

臺灣大學光電工程學研究所學位論文

Volume or Term/Year and Month of Publication

2008年

Academic Degree Category

碩士

Advisor

張宏鈞

Content Language

英文

Chinese Abstract

在本篇論文中,我們發展了新的電磁解析模型來分析一維波導、二維波導和光子晶體結構。同時,我們也推導了新的同步近似項的邊界條件來配合新的方程式。另外,由於樂勤得(Legendre)內差方程式準確的近似特性,我們也使用了譜方法中的樂勤得法來對空間分割。 不同於大部分是以赫姆霍茲(Helmholtz)方程式為基礎來解波導模態,我們提出了一種綜合了馬克斯威爾旋度以及散度方程式的新模型來產生特徵值(eigenvalue)問題的方法。對於在介面的邊界條件,我們則是使用了同步近似項的邊界條件,這種邊界條件可以在數學上證明具有穩定數值的效果。雖然在時域上,使用同步近似項邊界條件的譜方法是一個被證實具有不錯收斂性質的數值方法,但其使用的方程式跟本篇研究所發展的、在頻域使用的方程式有本質上的不同。在時域上我們所考慮的電磁場一律是實數,但在頻域上因為使用相位法的因素,所有的電磁場必須使用複數來表示。所以我們重新推導了分別給一維波導、二維波導和光子晶體使用的同步近似項邊界條件。 而為了確認新提出演算法的效率程度和數值收歛性,我們做了一些數值實例分析。這些光實例包括平板波導、部份填滿波導、圓柱形波導、具有尖角的埋入式波導、肋型波導以及四角晶格和三角晶格這兩種光子晶體結構。利用新的演算法來解這些結構時,我們可以得到高精度的傳播常數和特徵頻率以及觀察到指數收斂的特性。另外要特別說明的是,對於有尖角的介質波導,我們在沒有對邊界作任何特殊處理的情況下,得到了比以往更準確的數值解,這代表了在光波導分析領域的一項重大突破。

English Abstract

In this thesis, mode solvers for one-dimensional (1D) and 2D waveguides, and photonic crystals (PCs) with new electromagnetic formulations are developed. The new penalty-type boundary conditions are derived to work with new formulations, and the pseudospectral Legendre method is adopted to perform spatial discretization for its accurate approximation property. Unlike many waveguide mode solvers which are based on Helmholtz equations, we propose new formulations which combine Maxwell’s curl and divergence equations to derive the eigenvalue problem. For the interface boundary condition treatment, penalty-type boundary conditions are employed and mathematically proved that they can stablize the scheme. Athlough pseudospectral time-domain (PSTD) methods with penalty-type boundary conditions have been known to offer good convergence property, the related frequency-domain formulations developed in this work prossess intrinsic difference. In time-domain simulations, the electromagnetic fields considered are all real quantities, while in frequency-domain analysis, the fields are complex ones with the phasor technique. And new penalty-type boundary conditions in the frequency-domain mode analysis of 1D and 2D waveguides, and PCs, are respectively derived. Numerical examples are considered to examine the efficieney and numerical convergency property of the proposed algorithms. Optical structures in these examples include slab waveguides, partial-filled waveguide, fiber waveguides, channel waveguides with sharp corners, rib waveguides, PCs with square lattice, and PCs with triangular lattice. Spectral convergence property with very high-accuracy modal effective index and igenfrequency calculation is achieved. In particular, for the dielectric waveguide with corners, higher numerical accuracy than reported results is obtained without doing field singularity treatment at the corners as in the latter. This represents significant advancement in the numerical analysis of optical waveguide problems.

Topic Category 電機資訊學院 > 光電工程學研究所
工程學 > 電機工程
Reference
  1. [1] Baken, N. H. G, M. B. J. Diemeer, J. M. V. Splunter and H. Blok, “Computational modeling of diffused channel waveguides using a domain integral equation,” J. Lightwave Technol., vol. 8, pp. 576–586, 1990.
    連結:
  2. [2] Berenger, J.-P., “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys., vol. 114, pp. 185–200, 1994.
    連結:
  3. [3] Bierwirth, K., N. Schulz, and F. Arndt, “Finite-difference analysis of rectangular dielectric waveguides structures,” IEEE Trand. Microwave Theory Tech., vol. 34, pp. 1104–1113, 1986.
    連結:
  4. [5] Cendes, Z.J., and P. Silvester, “Numerical solution of dielectric loaded waveguides: Finite-Element Analysis,” IEEE Trans. Microwave Theory Tech., vol. MTT-18, pp. 1124–1131, 1970.
    連結:
  5. [7] Collin, R., Field Theory of Guided Waves. New York: McGraw-Hill, 1960.
    連結:
  6. [8] Doncker, Ph. De, “The use of transfinite elements in the methods of moments applied to electromagnetic scattering by dielectric cylinders,” Progress In Electromagnetics Research, PIER 25, pp. 77–94, 2000.
    連結:
  7. [9] Funaro, D., and D. Gottlieb, “A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations,” Math. Comp., vol. 51, pp. 599–613, 1988.
    連結:
  8. [10] Hadley, G. R., “Low-truncation-error finite difference equations for photonics simulation I: Beam propagation,” J. Lightwave Technol., vol. 16, pp. 134–141, 1998.
    連結:
  9. [11] Hadley, G. Ronald, “High-accuracy finite-difference equations for dielectric waveguide analysis: II. Dielectric corners,” J. Lightwave Technol., vol. 20, pp. 1219–1231, 2002.
    連結:
  10. [12] Hadley, G. R., and R. E. Smith, “Full-vector waveguide modeling using an iterative finite-difference method with transparent boundary conditons,” J. Lightwave Technol., vol. 13, pp. 465–469, 1999.
    連結:
  11. [13] Hesthaven, J. S., and D. Gottlieb, “A stable penalty method for the compressible Navier-Stokes equations. I. Open boundary conditions,” SIAM J Sci. Comput., vol. 17, pp. 579–612, 1996.
    連結:
  12. [14] Hesthaven, J. S., P. G. Dinesen, and J. P. Lynov, “Spectral collocation time-domain modeling of diffractive optical elements,” J. Comput. Phys., vol. 155, pp. 287–306, 1999.
    連結:
  13. [15] Hesthaven, J., S. Gottlieb, and D. Gottlieb, Spectral Methods for Time-dependent Problems. Cambridge University Press, 2007.
    連結:
  14. [16] Johnson, S. G., and J. D. Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Opt. Express, vol. 8, pp. 173–190, 2001.
    連結:
  15. [17] Lee, J. F., D. K. Sun, and Z. J. Cendes, “Full-wave analysis of dielectric waveguides using tangential vector finite elements,” IEEE Trans. Microwave Theory Tech., vol. 39, pp. 1262–1271, 1991.
    連結:
  16. [18] L ̈sse, P., P. Stuwe, J. Sch ̈le, and H.-G. Unger, “Analysis of vectorial mode fields in optical waveguides by an new finite difference method,” J. Lightwave Technol., vol. 12, pp. 487–493, 1994.
    連結:
  17. [19] Orszag, S. A., and G. S. Patterson, “Numerical simulation of three-dimensional homogeneous isotropic turbulence,” Phys. Rev. Lett., vol. 28, pp. 76–79, 1972.
    連結:
  18. [20] Rahman, B. M. A., and J. B. Davis, “Finite-element analysis of optical and microwave waveguides problems,” IEEE Trans. Microwave Theory Tech., vol. MTT-32, pp. 20–28, 1984.
    連結:
  19. [21] Sphicopoulos, T., V. Teodoridis, and F. E. Gardiol, “Dyadic Green function for the eletromagnetic field in multilayered isotropic media: an operator approach,” Inst. Elec. Eng. Proc., vol. 132, pp. 392–335, 1985.
    連結:
  20. [22] Stern, M. S., P. C. Kendall, and P. W. A. Mcllroy, “Analysis of the aspectral index method for vector modes of rib waveguides,” Inst. Elec. Eng. Proc. -J., vol. 137, pp. 21–26, 1990.
    連結:
  21. [23] Sudbo, A. S., “Why are accurate computations of mode fields in rectangular dielectric waveguides difficult?,” J. Lightw. Technol., vol. 10, pp. 418–419, 1992.
    連結:
  22. [24] Teng, C. H., B. Y. Lin, H. C. Chang, H. C. Hsu, C. N. Lin, and K. A. Feng, “A Legendre pseudospectral penalty scheme for solving time-domain Maxwell’s equations,” J. Sci. Comput., 2008. (in press)
    連結:
  23. [25] Thomas, N., P. Sewell, and T. M. Benson, “A new full-vectorial higher order finite-difference scheme for the modal analysis of rectangular dielectric waveguides,” J. Lightwave Technol., vol. 25, pp. 2563–2570, 2007.
    連結:
  24. [27] Yang, B., and J. S. Hesthaven, “A pseudospectral method for time-domain computation of electromagnetic scattering by bodies of revolution,” IEEE Trans. Antennas Propagat., vol. 47, pp. 132–141, 1999.
    連結:
  25. [28] Yu, C. P., and H. C. Chang, “A compact finite-difference frequency-domain method for the analysis of two-dimensional photonic crystals,” Opt. Express, vol. 12, pp. 1397–1408, 2004.
    連結:
  26. [4] Boyd, P. J., Chebyshev and Fourier Spectral Methods. New York, NY: Dover Publications, 1999.
  27. [6] Chiang, P.J., Development of Multidomain Pseudospectral Mode Solvers for Optical Waveguides and Photonic Crystals, M. S. Thesis, Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei, Taiwan, January 2007.
  28. [26] Wu, J. L., Development of Multidomain Pseudospectral Mode Solvers for Optical Waveguides and Photonic Crystals, M. S. Thesis, Graduate Institute of Communication Engineering, National Taiwan University, Taipei, Taiwan, June 2003.