透過您的圖書館登入
IP:52.90.121.17
  • 學位論文

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

New Pseudospectral Penalty Models for Electromagnetic Mode Problems

指導教授 : 張宏鈞
若您是本文的作者,可授權文章由華藝線上圖書館中協助推廣。

摘要


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

關鍵字

譜方法 波導 光子晶體 模態

並列摘要


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.

參考文獻


[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] Berenger, J.-P., “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys., vol. 114, pp. 185–200, 1994.
[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.
[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.
[7] Collin, R., Field Theory of Guided Waves. New York: McGraw-Hill, 1960.

延伸閱讀