×

A fourth order finite difference method for waveguides with curved perfectly conducting boundaries. (English) Zbl 1231.78045

Summary: A novel high order finite difference method is introduced for optical waveguides with smoothly curved perfectly electric conducting (PEC) boundaries. The proposed method shares some similarities with our previous matched interface and boundary (MIB) methods developed for treating dielectric interfaces of optical waveguides, such as the use of a simple Cartesian grid, the standard finite difference schemes, and fictitious values. However, the PEC boundary conditions have a physical nature quite different from that of the jump conditions at the dielectric interfaces, i.e., all six electric and magnetic field components are prescribed in the jump conditions, while only three of them are known at the PEC walls. Consequently, the previously developed MIB methods are not applicable to deal with the perfectly conducting boundaries. To overcome this difficulty, a novel ray-casting fictitious domain method is constructed to enforce the PEC conditions along the normal direction. Such a boundary implementation couples the transverse magnetic field components so that the resulting ray-casting MIB method is a full vectorial approach for the modal analysis of optical waveguides. The new MIB method is validated by considering both homogeneous and inhomogeneous waveguides. Numerical results confirm the designed fourth order of accuracy.

MSC:

78M20 Finite difference methods applied to problems in optics and electromagnetic theory
78A50 Antennas, waveguides in optics and electromagnetic theory

Software:

SRRIT
Full Text: DOI

References:

[1] Hadley, G. R., Numerical simulation of waveguides of arbitrary cross-section, Int. J. Electron. Commun., 58, 86-92 (2004)
[2] Chiang, P.-J.; Wu, C.-L.; Teng, C.-H.; Yang, C.-S.; Chang, H.-C., Full-vectorial optical waveguide mode solvers using multidomain pseudospectral frequency-domain (PSFD) formulations, IEEE J. Quantum Electron., 44, 56-66 (2008)
[3] Koshiba, M.; Tsuji, Y., Curvilinear hybrid edge/nodal elements with triangular shape for guide-wave problems, J. Lightwave Technol., 18, 737-743 (2000)
[4] Dossou, K.; Fontaine, M., A high order isoparametric finite element method for the computation of waveguide modes, Comput. Method Appl. Mech. Engrg., 194, 837-858 (2005) · Zbl 1073.78009
[5] Marais, N.; Davidson, D. B., Numerical evaluation of hierarchical vector finite elements on curvilinear domains in 2D, IEEE Trans. Antenn. Propag., 54, 734-738 (2006)
[6] Notaros, B. M., Higher order frequency-domain computational electromagnetics, IEEE Trans. Antenn. Propag., 56, 2251-2276 (2008) · Zbl 1369.78486
[7] Dey, S.; Mittra, R., A locally conformal finite-difference time-domain (FDTD) algorithm for modeling three-dimensional perfectly conducting objects, IEEE Microwave Guided Wave Lett., 7, 273-275 (1997)
[8] Yu, W.; Mittra, R., A conformal FDTD algorithm for modeling perfectly conducting objects with curved-shaped surfaces and edges, Microwave Opt. Tech. Lett., 27, 136-138 (2000)
[9] Xiao, T.; Liu, Q. H., Enlarged cells for the conformal FDTD method to avoid the time step reduction, IEEE Microwave Compon. Lett., 14, 551-553 (2004)
[10] Benkler, S.; Chavannes, N.; Kuster, N., A new 3D conformal PEC FDTD scheme with user-defined geometric precision and derived stability criterion, IEEE Trans. Antenn. Propag., 54, 1843-1849 (2006)
[11] Bettencourt, M. T., Flux limiting embedded boundary technique for electromagnetic FDTD, J. Comput. Phys., 227, 3134-3158 (2008) · Zbl 1134.78017
[12] Tornberg, A.-K.; Engquist, B., Consistent boundary conditions for the Yee scheme, J. Comput. Phys., 227, 6922-6943 (2008) · Zbl 1154.65354
[13] Zhao, S.; Wei, G. W., High-order FDTD methods via derivative matching for Maxwell’s equations with material interfaces, J. Comput. Phys., 200, 60-103 (2004) · Zbl 1050.78018
[14] Hadley, G. R., High-accuracy finite-difference equations for dielectric waveguide analysis II: Dielectric corners, J. Lightwave Technol., 20, 1219-1231 (2002)
[15] Thomas, N.; Sewell, P.; Benson, T. M., A new full-vectorial higher order finite-difference scheme for the modal analysis of rectangular dielectric waveguides, J. Lightwave Technol., 25, 2563-2570 (2007)
[16] Zhao, S., Full-vectorial matched interface and boundary (MIB) method for the modal analysis of dielectric waveguides, J. Lightwave Technol., 26, 2251-2259 (2008)
[17] Chiang, Y. C.; Chiou, Y. P.; Chang, H. C., Improved full-vectorial finite-difference mode solver for optical waveguides with step-index profiles, J. Lightwave Technol., 20, 1609-1618 (2002)
[18] Zhao, S., High order vectorial analysis of waveguides with curved dielectric interfaces, IEEE Microwave Compon. Lett., 19, 266-268 (2009)
[19] Zhao, S., High order matched interface and boundary methods for the Helmholtz equation in media with arbitrarily curved interfaces, J. Comput. Phys., 229, 3155-3170 (2010) · Zbl 1187.78044
[20] Zhao, S., On the spurious solutions in the high-order finite difference methods, Comput. Method Appl. Mech. Engrg, 196, 5031-5046 (2007) · Zbl 1173.74462
[21] Zhao, S.; Wei, G. W., Matched interface and boundary (MIB) method for the implementation of boundary conditions in high-order central finite differences, Int. J. Numer. Method Engrg., 77, 1690-1730 (2009) · Zbl 1161.65080
[22] Zhou, Y. C.; Zhao, S.; Feig, M.; Wei, G. W., High order matched interface and boundary (MIB) schemes for elliptic equations with discontinuous coefficients and singular sources, J. Comput. Phys., 213, 1-30 (2006) · Zbl 1089.65117
[23] Harrington, R. F., Time-Harmonic Electromagnetic Fields (2001), Wiley-IEEE Press
[24] Bai, Z.; Stewart, G. W., SRRIT: a Fortran subroutine to calculate the dominant invariant subspace of a nonsymmetric matrix, ACM Trans. Math. Software, 23, 494-513 (1997) · Zbl 0912.65029
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.