About this Journal Submit a Manuscript Table of Contents 
International Journal of Antennas and Propagation
Volume 2012 (2012), Article ID 528037, 6 pages
doi:10.1155/2012/528037
Research Article

A High-Order Compact 2D FDFD Method for Waveguides

Gang Zheng1 and Bing-Zhong Wang2

1State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, State Oceanic Administration, Hangzhou 310012, China
2Institute of Applied Physics, University of Electronic Science and Technology of China, Chengdu 610054, China

Received 7 May 2012; Accepted 15 July 2012

Academic Editor: Zhongxiang Q. Shen

Copyright © 2012 Gang Zheng and Bing-Zhong Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

A high-order compact two-dimensional finite-difference frequency-domain (2D FDFD) method is proposed for the analysis of the dispersion characteristics of waveguides. A surface impedance boundary condition (SIBC) for the high-order compact 2D FDFD method is also given to model lossy metal waveguides. Four transverse field components are involved in the final eigenequation. Numerical examples are given, which show that this high-order compact 2D FDFD method is more efficient than the low-order compact 2D FDFD method and has a less storage cost.

1. Introduction

In practical engineering designs, it is very important to accurately and efficiently analyze the dispersion characteristics of waveguides. A two-dimensional finite-difference frequency-domain (2D FDFD) method was proposed for this purpose [1]. A compact 2D FDFD method was brought forward later in which only four transverse field components are involved in the final eigenequation [2]. In 2D FDFD methods, the cross-section of wave guides is discretized with compact Yee’s meshes, and corresponding difference equations are set up on the nodes. A matrix eigenequation can be established from them finally. The dispersion characteristics can be obtained by solving this eigenequation. Unlike in two-dimensional finite-difference time-domain (2D FDTD) methods [3, 4], in 2D FDFD methods the complex propagation constant can be found at a given frequency directly and there is no need of the discrete Fourier transform. Another advantage of 2D FDFD methods is that the dispersion characteristics of several modes at a given frequency can be analyzed at the same time.

In a low-order finite-difference method, fine meshes need to be used in order to obtain a high accuracy, and this is very time consuming. However, under the same spatial discretization, the accuracy can be largely improved by a high order finite-difference method.

In 2D FDFD methods, the main computational time is spent on solving the final matrix eigenequation, and it increases as the size of the final matrix eigenequation increases. In this letter, a high-order compact 2D FDFD method is proposed for calculating complex propagation constants of waveguides, and a corresponding surface impedance boundary condition (SIBC) is also given to model lossy metal waveguides for the high-order compact 2D FDFD method. The numerical results show that under the same accuracy, the number of the meshes in the high order compact 2D FDFD method is much less than in the low-order compact 2D FDFD method. Correspondingly, both the computational time and the number of nonzero elements largely decrease in the former, and the burdens of computation and storage are reduced.

2. Formulation

2.1. The High-Order Compact 2D FDFD Method

Electric field 𝐄 and magnetic field 𝐇 are normalized with the square root of the free-space wave impedance ( 𝜂 0 = √ 𝜇 0 / 𝜀 0 ). It is assumed that waveguides are uniform in 𝑧 -direction. Then the field of a mode in a waveguide can be expressed as 𝐄 𝐸 ( 𝑥 , 𝑦 , 𝑧 ) = 𝑥 ( 𝑥 , 𝑦 ) 𝐞 𝑥 + 𝐸 𝑦 ( 𝑥 , 𝑦 ) 𝐞 𝑦 + 𝐸 𝑧 ( 𝑥 , 𝑦 ) 𝐞 𝑧 𝑒 − 𝛾 𝑧 , 𝐻 𝐇 ( 𝑥 , 𝑦 , 𝑧 ) = 𝑥 ( 𝑥 , 𝑦 ) 𝐞 𝑥 + 𝐻 𝑦 ( 𝑥 , 𝑦 ) 𝐞 𝑦 + 𝐻 𝑧 ( 𝑥 , 𝑦 ) 𝐞 𝑧 𝑒 − 𝛾 𝑧 , ( 1 ) where 𝛾 is the complex propagation constant of the mode. After substituting (1) into Maxwell’s curl equations, the following equations are obtained: − 𝜕 𝐸 𝑧 𝜕 𝑥 + 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑦 = 𝛾 𝐸 𝑥 , − ( 2 ) 𝜕 𝐸 𝑧 𝜕 𝑦 − 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑥 = 𝛾 𝐸 𝑦 , − ( 3 ) 𝜕 𝐻 𝑧 𝜕 𝑥 − 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑦 = 𝛾 𝐻 𝑥 , − ( 4 ) 𝜕 𝐻 𝑧 𝜕 𝑦 + 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑥 = 𝛾 𝐻 𝑦 , ( 5 ) 𝜕 𝐻 𝑥 − 𝜕 𝑦 𝜕 𝐻 𝑦 𝜕 𝑥 + 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑧 = 0 , ( 6 ) 𝜕 𝐸 𝑦 − 𝜕 𝑥 𝜕 𝐸 𝑥 𝜕 𝑦 + 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑧 = 0 . ( 7 )

The cross-section of the waveguide is discretized with compact Yee’s meshes (see Figure 1). After the forth-order central difference is used to approximate the partial derivatives in (2)–(7), the following difference equations are obtained: 1 𝐸 2 4 Δ 𝑥 𝑧 ( 𝑖 + 2 , 𝑗 ) − 𝐸 𝑧 − 9 ( 𝑖 − 1 , 𝑗 ) 𝐸 8 Δ 𝑥 𝑧 ( 𝑖 + 1 , 𝑗 ) − 𝐸 𝑧 ( 𝑖 , 𝑗 ) + 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑦 ( 𝑖 , 𝑗 − 1 ) = 𝛾 𝐸 𝑥 1 ( 𝑖 , 𝑗 ) , ( 8 ) 𝐸 2 4 Δ 𝑦 𝑧 ( 𝑖 , 𝑗 + 2 ) − 𝐸 𝑧 − 9 ( 𝑖 , 𝑗 − 1 ) 𝐸 8 Δ 𝑦 𝑧 ( 𝑖 , 𝑗 + 1 ) − 𝐸 𝑧 ( 𝑖 , 𝑗 ) − 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑥 ( 𝑖 − 1 , 𝑗 ) = 𝛾 𝐸 𝑦 1 ( 𝑖 , 𝑗 ) , ( 9 ) 𝐻 2 4 Δ 𝑥 𝑧 ( 𝑖 + 2 , 𝑗 ) − 𝐻 𝑧 ( − 9 𝑖 − 1 , 𝑗 ) 𝐻 8 Δ 𝑥 𝑧 ( 𝑖 + 1 , 𝑗 ) − 𝐻 𝑧 ( 𝑖 , 𝑗 ) − 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑦 ( 𝑖 + 1 , 𝑗 ) = 𝛾 𝐻 𝑥 ( 1 𝑖 , 𝑗 ) , ( 1 0 ) 𝐻 2 4 Δ 𝑦 𝑧 ( 𝑖 , 𝑗 + 2 ) − 𝐻 𝑧 − 9 ( 𝑖 , 𝑗 − 1 ) 𝐻 8 Δ 𝑦 𝑧 ( 𝑖 , 𝑗 + 1 ) − 𝐻 𝑧 ( 𝑖 , 𝑗 ) + 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑥 ( 𝑖 , 𝑗 + 1 ) = 𝛾 𝐻 𝑦 1 ( 𝑖 , 𝑗 ) , ( 1 1 ) 𝐻 2 4 Δ 𝑦 𝑥 ( 𝑖 − 1 , 𝑗 + 1 ) − 𝐻 𝑥 − 9 ( 𝑖 − 1 , 𝑗 − 2 ) 𝐻 8 Δ 𝑦 𝑥 ( 𝑖 − 1 , 𝑗 ) − 𝐻 𝑥 − 1 ( 𝑖 − 1 , 𝑗 − 1 ) 𝐻 2 4 Δ 𝑥 𝑦 ( 𝑖 + 1 , 𝑗 − 1 ) − 𝐻 𝑦 + 9 ( 𝑖 − 2 , 𝑗 − 1 ) 𝐻 8 Δ 𝑥 𝑦 ( 𝑖 , 𝑗 − 1 ) − 𝐻 𝑦 ( 𝑖 − 1 , 𝑗 − 1 ) − 𝑗 𝑘 0 𝜀 𝑟 𝐸 𝑧 1 ( 𝑖 , 𝑗 ) = 0 , ( 1 2 ) 𝐸 2 4 Δ 𝑥 𝑦 ( 𝑖 + 2 , 𝑗 ) − 𝐸 𝑦 − 9 ( 𝑖 − 1 , 𝑗 ) 𝐸 8 Δ 𝑥 𝑦 ( 𝑖 + 1 , 𝑗 ) − 𝐸 𝑦 − 1 ( 𝑖 , 𝑗 ) 𝐸 2 4 Δ 𝑦 𝑥 ( 𝑖 , 𝑗 + 2 ) − 𝐸 𝑥 + 9 ( 𝑖 , 𝑗 − 1 ) 𝐸 8 Δ 𝑦 𝑥 ( 𝑖 , 𝑗 + 1 ) − 𝐸 𝑥 ( 𝑖 , 𝑗 ) − 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑧 ( 𝑖 , 𝑗 ) = 0 , ( 1 3 ) where Δ 𝑥 and Δ 𝑦 are mesh sizes in 𝑥 - and 𝑦 -directions, respectively. These high-order difference equations are set up on all corresponding nodes except the nodes near the interface between two materials. Mixed-order difference equations are set up on these nodes, in which the second-order central difference approximation is used for the normal derivative and the forth-order central difference approximation is used for the tangential derivative.

528037.fig.001
Figure 1: Compact Yee’s meshes.
2.2. The SIBC for the High-Order Compact 2D FDFD Method

For lossy metal, the tangential electric field 𝐄 t a n and the tangential magnetic field 𝐇 t a n on its surface satisfy 𝐄 t a n = 𝑍 𝑠 𝐧 × 𝐇 t a n , ( 1 4 ) where 𝐧 is a unit vector normal to the surface of lossy metal (see Figure 2). The normalized surface impedance 𝑍 𝑠 is 𝑍 𝑠 = ( 1 + 𝑗 ) 𝜎 𝛿 𝜂 0 , ( 1 5 ) where 𝜎 and 𝛿 are the conductivity and the skin depth of lossy metal, respectively.

fig2
Figure 2: (a) A surface of lossy metal. (b) A corner of lossy metal.

The surface in Figure 2(a) is considered here, and the other surfaces of lossy metal can be treated in the same way. It can be seen that the 𝐄 t a n nodes of meshes coincide with the surface of lossy metal, but the 𝐇 t a n nodes of meshes do not coincide with the surface. Therefore, the following approximation is taken 𝐻 𝑦 ( 0 , 𝑗 ) ≈ 𝐻 𝑦 𝐻 ( 1 , 𝑗 ) , ( 1 6 ) 𝑧 ( 0 , 𝑗 ) ≈ 𝐻 𝑧 ( 1 , 𝑗 ) . ( 1 7 )

The discrete form of (14) on the surface is 𝐸 𝑦 ( 1 , 𝑗 ) = − 𝑍 𝑠 𝐻 𝑧 𝐸 ( 0 , 𝑗 ) , ( 1 8 ) 𝑧 ( 1 , 𝑗 ) = 𝑍 𝑠 𝐻 𝑦 ( 0 , 𝑗 − 1 ) . ( 1 9 )

After substituting (16) into (19), the SIBC for 𝐸 𝑧 on the surface is obtained as follows 𝐸 𝑧 ( 1 , 𝑗 ) = 𝑍 𝑠 𝐻 𝑦 ( 1 , 𝑗 − 1 ) . ( 2 0 )

Now the SIBC for the remaining tangential component of the electric field on the surface is considered. After substituting (17) into (18), we have 𝐸 𝑦 ( 1 , 𝑗 ) = − 𝑍 𝑠 𝐻 𝑧 ( 1 , 𝑗 ) . ( 2 1 ) The mixed-order difference equation is set up on the node 𝐻 𝑧 ( 1 , 𝑗 ) as 1 𝐸 Δ 𝑥 𝑦 ( 2 , 𝑗 ) − 𝐸 𝑦 + 1 ( 1 , 𝑗 ) 𝐸 2 4 Δ 𝑦 𝑥 ( 1 , 𝑗 + 2 ) − 𝐸 𝑥 − 9 ( 1 , 𝑗 − 1 ) 𝐸 8 Δ 𝑦 𝑥 ( 1 , 𝑗 + 1 ) − 𝐸 𝑥 ( 1 , 𝑗 ) + 𝑗 𝑘 0 𝜇 𝑟 𝐻 𝑧 ( 1 , 𝑗 ) = 0 . ( 2 2 ) Substituting (22) into (21) to eliminate 𝐻 𝑧 ( 1 , 𝑗 ) , then the SIBC for 𝐸 𝑦 on the surface is obtained as 𝐸 𝑦 𝑍 ( 1 , 𝑗 ) = 𝑠 Δ 𝑥 𝑗 𝑘 0 𝜇 𝑟 Δ 𝑥 + 𝑍 𝑠 × ⎧ ⎪ ⎨ ⎪ ⎩ 𝐸 𝑦 ( 2 , 𝑗 ) − 9 Δ 𝑥 𝐸 8 Δ 𝑦 𝑥 ( 1 , 𝑗 + 1 ) − 𝐸 𝑥 + 1 ( 1 , 𝑗 ) 𝐸 2 4 Δ 𝑦 𝑥 ( 1 , 𝑗 + 2 ) − 𝐸 𝑥 ⎫ ⎪ ⎬ ⎪ ⎭ . ( 1 , 𝑗 − 1 ) ( 2 3 )

Now the SIBC for 𝐸 𝑥 ( 1 , 1 ) and 𝐸 𝑦 ( 1 , 1 ) in the corner (see Figure 2(b)) is considered. They can be expressed as 𝐸 𝑦 ( 1 , 1 ) = − 𝑍 𝑠 𝐻 𝑧 𝐸 ( 1 , 1 ) 𝑥 ( 1 , 1 ) = 𝑍 𝑠 𝐻 𝑧 ( 1 , 1 ) . ( 2 4 )

After combining (24) with the low-order difference equation on the node 𝐻 𝑧 ( 1 , 1 ) , the SIBC in this corner is obtained 𝐸 𝑥 ( 𝑍 1 , 1 ) = 𝑠 𝑗 𝑘 0 𝜇 𝑟 Δ 𝑥 Δ 𝑦 + 𝑍 𝑠 × ( Δ 𝑥 + Δ 𝑦 ) Δ 𝑥 𝐸 𝑥 ( 1 , 2 ) − Δ 𝑦 𝐸 𝑦 𝐸 ( 2 , 1 ) 𝑦 𝑍 ( 1 , 1 ) = 𝑠 𝑗 𝑘 0 𝜇 𝑟 Δ 𝑥 Δ 𝑦 + 𝑍 𝑠 × ( Δ 𝑥 + Δ 𝑦 ) Δ 𝑦 𝐸 𝑦 ( 2 , 1 ) − Δ 𝑥 𝐸 𝑥 . ( 1 , 2 ) ( 2 5 )

2.3. The Final Eigenequation

From (8) to (13) and the SIBC, (26) is obtained as ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ 0 𝐴 1 , 2 𝐴 1 , 3 0 0 0 𝐴 1 , 7 𝐴 2 , 1 0 0 𝐴 2 , 4 0 0 0 0 𝐴 3 , 2 𝐴 3 , 3 𝐴 0 0 0 0 4 , 1 0 0 𝐴 4 , 4 𝐴 4 , 5 𝐴 4 , 6 0 𝐴 5 , 1 0 0 0 𝐴 5 , 5 𝐴 5 , 6 0 𝐴 6 , 1 0 0 0 0 𝐴 6 , 6 0 0 𝐴 7 , 2 0 0 0 0 𝐴 7 , 7 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ 𝑥 1 𝑥 2 𝑥 3 𝑥 4 𝑥 5 𝑥 6 𝑥 7 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ 𝑥 = 𝛾 1 𝑥 2 0 0 0 0 0 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ , ( 2 6 ) where 𝑥 1 , 𝑥 2 , 𝑥 3 , 𝑥 4 , 𝑥 5 , 𝑥 6 , and 𝑥 7 are all column vectors. The vector 𝑥 1 represents the transverse electric field inside the FDFD domain, the vector 𝑥 2 represents the transverse magnetic field inside the FDFD domain, the vector 𝑥 3 represents the 𝑧 -direction electric field inside the FDFD domain, the vector 𝑥 4 represents the 𝑧 -direction magnetic field inside the FDFD domain, the vector 𝑥 5 represents the transverse tangential component of the electric field on the surface of lossy metal, the vector 𝑥 6 represents the transverse tangential component of the electric field in the corners, and the vector 𝑥 7 represents the 𝑧 -direction electric field on the surface of lossy metal.

After eliminating 𝑥 3 , 𝑥 4 , 𝑥 5 , 𝑥 6 , and 𝑥 7 in (26), the following eigenequation is derived: ⎡ ⎢ ⎢ ⎣ 0 𝐵 1 , 2 𝐵 2 , 1 0 ⎤ ⎥ ⎥ ⎦ ⎡ ⎢ ⎢ ⎣ 𝑥 1 𝑥 2 ⎤ ⎥ ⎥ ⎦ ⎡ ⎢ ⎢ ⎣ 𝑥 = 𝛾 1 𝑥 2 ⎤ ⎥ ⎥ ⎦ , ( 2 7 ) where 𝐵 1 , 2 = 𝐴 1 , 2 − 𝐴 1 , 3 𝐴 − 1 3 , 3 𝐴 3 , 2 − 𝐴 1 , 7 𝐴 − 1 7 , 7 𝐴 7 , 2 , 𝐵 2 , 1 = 𝐴 2 , 1 − 𝐴 2 , 4 𝐴 − 1 4 , 4 ⎡ ⎢ ⎢ ⎣ 𝐴 4 , 5 𝐴 − 1 5 , 5 𝐴 5 , 6 𝐴 − 1 6 , 6 𝐴 6 , 1 − 𝐴 5 , 1 − 𝐴 4 , 6 𝐴 − 1 6 , 6 𝐴 6 , 1 + 𝐴 4 , 1 ⎤ ⎥ ⎥ ⎦ . ( 2 8 ) The matrix blocks 𝐴 3 , 3 , 𝐴 4 , 4 , 𝐴 5 , 5 , 𝐴 6 , 6 and 𝐴 7 , 7 are full rank and have nonzero elements only on their diagonal lines. So their inverses are very easy to obtain and the matrix in (27) is still very sparse. After solving (27), the complex propagation constants of the modes in the waveguide can be obtained.

3. Numerical Results

The dispersion characteristics of a rectangular metal waveguide with the size of 19.05 mm × 9 mm are respectively analyzed by the low-order compact 2D FDFD method and the high order compact 2D FDFD method, in order to compare their computational times and the numbers of nonzero elements of the matrices in the final eigenequations. Both of the two methods are programmed with Matlab 7.0 and running in a P4 desktop computer (CPU: 3.2 GHz).

Firstly, the metal is assumed to be lossless. The complex propagation constant of the TE10 mode is calculated at 141 frequency points by the low-order compact 2D FDFD method with 40 × 40 compact Yee’s meshes in the cross-section, the low-order compact 2D FDFD method with 32 × 32 compact Yee’s meshes and the high order compact 2D FDFD method with 16 × 16 compact Yee’s meshes, respectively. The analytical formula of the TE10 mode is used as a standard and the relative errors of the two methods are shown in Figure 3. The curve of the high order compact 2D FDFD method with the coarsest meshes is in a very good agreement with the curve of the low-order compact 2D FDFD method with the finest meshes. The total computational times of solving the final eigenequations at all frequency sample points and the numbers of nonzero elements of the matrices in the final eigenequations are shown in Table 1 for the two methods, respectively. It can be seen that both of them are the least in the high order compact 2D FDFD method without loss of accuracy. And these mean high computational efficiency and low storage cost. The complex propagation constants of other modes are also calculated, and the relative errors are shown in Figure 3, too.

tab1
Table 1: The computational times and the number of nonzero elements.
528037.fig.003
Figure 3: The relative errors of the complex propagation constants obtained by the two methods, respectively.

Secondly, the metal is assumed with a conductivity σ = 5.8 × 107 S/m. The complex propagation constants of mode 1 and mode 4 are calculated by the two methods, which are, respectively, corresponding to the TE10 mode and the TM11 mode of the lossless metal waveguide. The results of the two methods are compared with the results of Ansoft-HFSS, and they are shown in Figures 4 and 5 from the cutoff areas to the propagation areas. It can be seen that although the coarsest meshes are used in the high order compact 2D FDFD method, the results of it and Ansoft-HFSS are still in a very good agreement in both cases of mode 1 and mode 4. However, the curves of the low-order compact 2D FDFD method with coarser meshes are not consistent with the others, but its curves tend to the others as the meshes become finer.

fig4
Figure 4: (a) The phase constant of mode 1. (b) The attenuation constant of mode 1.
fig5
Figure 5: (a) The phase constant of mode 4. (b) The attenuation constant of mode 4.

4. Conclusion

In this letter, a high-order compact two-dimensional finite-difference frequency-domain (2D FDFD) method is proposed for analyzing the dispersion characteristics of waveguides. A corresponding surface impedance boundary condition (SIBC) is given to model lossy metal waveguides. Some low-order approximations are used in the SIBC. However, the high order compact 2D FDFD method with the SIBC should be more accurate than the low-order compact 2D FDFD method, because the forth-order central difference is used in the most part of computational domain. In the numerical examples, the high order compact 2D FDFD method is compared with the low-order compact 2D FDFD method, and both of the cases of lossless and lossy metal are tested. The results show that both phase constants and attenuation constants can be accurately calculated through the high order compact 2D FDFD method, but both of the computational time and the storage cost are largely reduced.

Acknowledgment

This work was supported by the Doctoral Program of Higher Education of China (no. 20060614005).

References

  1. M. L. Lui and Z. Chen, “A direct computation of propagation constant using compact 2-D full-wave eigen-based finite-difference frequency-domain technique,” in Proceedings of the International Conference on Computational Electromagnetics and Its Applications (ICCEA'99), pp. 78–81, Beijing, China, November 1999. View at Publisher · View at Google Scholar
  2. Y.-J. Zhao, K.-L. Wu, and K.-K. M. Cheng, “A compact 2-D full-wave finite-difference frequency-domain method for general guided wave structures,” IEEE Transactions on Microwave Theory and Techniques, vol. 50, no. 7, pp. 1844–1848, 2002. View at Publisher · View at Google Scholar · View at Scopus
  3. S. Xiao, R. Vahldieck, and H. Jin, “Full-wave analysis of guided wave structures using a novel 2-D FDTD,” IEEE Microwave and Guided Wave Letters, vol. 2, no. 5, pp. 165–167, 1992. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Asi and L. Shafai, “Dispersion analysis of anisotropic inhomogeneous waveguides using compact 2D-FDTD,” Electronics Letters, vol. 28, no. 15, pp. 1451–1452, 1992. View at Scopus