# Electromagnetic scattering at the waveguide step between equilateral triangular waveguides.

1. IntroductionNowadays radiofrequency systems are widely used for very diverse applications. The GPS system for positioning, the LTE system for mobile communications, the DBS system for digital broadcasting, and ground penetrating radar or remote security screening are some examples of the usage of the electromagnetic spectrum that extends from some GHz to several hundred GHz [1, 2]. The increasing demand of the quality of service has led to very restrictive system specifications, which are translated into more complex microwave and millimeter-wave devices. In order to accomplish such designs, CAD (Computer Aided Design) tools have experienced a strong development over the last decades [3, 4].

These tools solve Maxwell equations using numerical approaches or quasi-analytical methods [5]. Commercial tools as CST Microwave Studio or Ansys HFSS focus on the numerical solutions to Maxwell equations and develop methods as the Finite Element Method or the Finite Differences Time Domain Method. This kind of methods has the advantage of being extremely flexible when coping with problems that involve arbitrary geometries. In contrast, they can be extremely time-consuming compared to the quasi-analytical methods. This paper focuses on the latter type of methods (i.e., quasi-analytical methods), particularly in the Mode-Matching method [6]. This method was developed some decades ago and has been applied to classical waveguides with analytic modal spectrum such as the rectangular waveguide, the circular waveguide, and the elliptical waveguide for many types of waveguide devices such as transformers, filters, and couplers [7-11]. However, it has not been applied to the equilateral triangular waveguide yet.

This paper addresses the Mode-Matching method applied to the equilateral triangular geometry, giving a flexible mathematical formulation based on well-known plane waves. This formulation simplifies the surface integrals of the method and allows a very easy generalization to other problems involving waveguides whose modes can be written as a plane wave expansion. In order to validate the formulation, several case studies are presented and compared with the results given by other different numerical methods from commercial software.

2. Modal Fields in the Equilateral Triangular Waveguide

The electromagnetic field enclosed in an arbitrary cross section waveguide can be described as an orthogonal series of modes. When the considered transmission medium exhibits translation symmetry along the propagation axis the electric ([??]) and magnetic ([??]) field can be expressed through their transversal and longitudinal components [12-14]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1b)

The terms [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are called the modal amplitudes at an arbitrary z plane, of which the complex variables [a.sub.p] and [b.sub.p] are determined by the excitation and load of the waveguide. The propagation constant [[gamma].sub.p] and the transversal fields [[??].sub.p] and [[??].sub.p] of each mode depend on the geometry of the waveguides under consideration. The geometries of a classical circular waveguide and a classical rectangular waveguide are shown in Figures 1(a) and 1(b), respectively. In this paper the equilateral triangular waveguide, shown in Figure 1(c), is analyzed. As it is known, for uniform cross-section waveguides made up of a single enclosing metallic conductor filled with a homogeneous dielectric material, the three-dimensional Maxwell equations simplify to the scalar 2D Helmholtz equation [13,14]:

[[DELTA].sub.t][[phi].sub.p] + [k.sup.2.sub.cp][[phi].sub.p] = 0, (2)

which must be held in the inner cross section of the waveguide under consideration. From this function, directly related to the z-field component of the electric or magnetic field, the full electromagnetic field is obtained. The expressions of the modal components [[??].sub.p], [[??].sub.p], [e.sub.zp], and [h.sub.zp] in (1a)-(1b) are included in Appendix. The orthogonal modes that constitute the complete electromagnetic solution for a single conductor waveguide can be divided into two families: the so-called TE (Transversal Electric) modes and TM (Transversal Magnetic) modes. The former are calculated imposing homogeneous Neumann-type boundary conditions in (2) along the contour C of the waveguide (see Figure 1(d)), whereas the latter are calculated imposing homogeneous Dirichlet-type boundary conditions:

[partial derivative][[phi].sub.p]/[partial derivative]n[|.sub.c] = 0, TE modes, (3a) [[phi].sub.p][|.sub.C] = 0, TM modes. (3b)

An initial characterization of the triangular waveguide can be found in [15,16]. The scalar function [[phi].sub.p] for the complete solution in an equilateral triangular waveguide is found in [17,18], where the complete set of solutions are divided not only into TE and TM modes but also into symmetries with respect to one of its medians. These functions referring to the coordinate system selected in this paper (see Figure 1(d)), using superscript S for symmetric solutions and A for asymmetric solutions, are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4d)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)

In the preceding definitions, variable b in constants [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] represents the height of the equilateral cross section of the waveguide. It is related to the side e of the equilateral triangular waveguide through b = e[square root of 3]/2. Each modal solution that belongs to the series (1a)-(1b) is denoted by subindex p ([[??].sub.p], [[??].sub.p]), which identifies a single pair (m, n) and one of the four families that have been introduced (T[E.sup.S], T[E.sup.A], T[M.sup.S], and T[M.sup.A]). T[E.sup.S.sub.m] modes arise with subscripts m [greater than or equal to] n [greater than or equal to] 0, T[E.sup.A.sub.mn] modes with subscripts m > n [greater than or equal to] 0, T[M.sup.S.sub.m] modes with m [greater than or equal to] n > 0, and T[M.sup.A.sub.mn] modes with m > n > 0. It can be seen that (4a)-(4b) fulfill homogeneous Neumann boundary conditions (3a) and that (4c)-(4d) fulfill homogeneous Dirichlet boundary conditions (3b). Finally, the cut-off wavenumber [k.sub.cp] in (2) can be calculated for each mode using

[k.sub.cp] = (4[pi]/3e) [square root of ([m.sup.2] + [n.sup.2] + mn)], (6)

which shows that symmetric and asymmetric TE and TM modes are degenerated (share the same cut-off wavenumber value). It can be easily seen that the cut-off wavenumber [k.sub.cp] can also be derived as:

[k.sup.2.sub.cp] = [[??].sup.2.sub.xi] + [[??].sup.2.sub.yi], i = 1, 2, 3. (7)

The cut-off frequency associated with each particular mode, which indicates the minimum frequency where the propagation constant associated with that mode has null attenuation (i.e., it is purely imaginary) for lossless media [13, 14], is defined as

[f.sub.cp] = ([k.sub.cp]/2[pi] [square root of ([mu][epsilon])]), (8)

where [mu] is the magnetic permeability and [epsilon] is the electric permittivity of the transmission medium filling the waveguide (being [([square root of ([mu]e)]).sup.-1] = 2.9979 * [10.sup.8] m/s for the vacuum). The fundamental modes (i.e., those with lowest cut-off frequency) are T[E.sup.A.sub.10] and T[E.sup.S.sub.10]. The cut-off wavenumbers, the cut-off frequencies, and the associated families of the first fifty modes in a hollow waveguide of side e = 1 cm are given in Table 1. This table can be used for any other equilateral triangular waveguide since [k.sub.c] and [f.sub.c] only need to be scaled by the new side e'.

The electric and magnetic field patterns of the first twelve modes are shown in Figure 2. Electric fields (in blue arrows) and magnetic fields (in green arrows) are perpendicular to each other. Besides, the former are perpendicular to the sides of the equilateral triangle, whereas the latter are tangential, as it is required by the boundary conditions of a perfect conductor. Symmetric TE modes (T[E.sup.S.sub.mn]) and asymmetric TM modes (T[M.sup.A.sub.mn]) show perfect electric wall (PEW) symmetry with respect to x = 0. On the other hand, asymmetric TE modes (T[E.sup.A.sub.mn]) and symmetric TM modes (T[M.sup.S.sub.m]n) show perfect magnetic wall (PMW) symmetry at x = 0.

3. Overview of the Mode-Matching Method

The main goal of the Mode-Matching method is to solve the modal amplitudes of (1a)-(1b) and relate them through the Generalized Scattering Matrix (GSM) [9, 19], used to characterize the full response by cascading individual building blocks also in antenna problems [20, 21]. In this section an overview of this method is introduced regardless of the cross section geometry of the waveguides involved in the analyzed discontinuity. In Section 4 this method will be particularized to our geometry of interest (i.e., the equilateral triangular waveguide).

The Mode-Matching method is based on the fact that the total transversal fields ([[??].sub.t], [[??].sub.t]) must be continuous across the transversal surface in the waveguide step under consideration. This surface is shown in Figure 3(a), where the waveguide (w) has been ended with a perfect conductor wall except in the intersected area with the smaller waveguide (s), and both waveguides are depicted without conductor thickness. According to (1a) and (1b), the transversal electromagnetic field at both sides of a step or discontinuity placed at the plane z = 0 is represented by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9d)

where superscripts (w) and (s) refer to each involved waveguide (see Figure 3(b) for clarification). It is important to highlight that, in order to formulate an affordable computational problem, in the preceding equations the number of modes in each waveguide has been truncated to [N.sub.w] and [N.sub.s], respectively. The boundary conditions applied across the aperture interface [A.sub.s] that preserve the transversal fields are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (10)

where EFBC stands for Electric Field Boundary Condition and MFBC stands for Magnetic Field Boundary Condition. The boundary conditions are imposed using a Galerkin method, where the EFBC is tested with the modal magnetic fields of waveguide (w) (i.e., [h.sup.(w).sub.j] and the MFBC is tested with the modal electric fields of waveguide (s) (i.e., [e.sup.(s).sub.i]) [11]. Arranging the modal amplitudes in column vectors ([a.sub.g] = [[*** [a.sup.(g).sub.p] ***].sup.T], [b.sub.g] = [[*** [b.sup.(g).sub.p] ***].sup.T], being g = w, s), this method leads to the following system of linear equations:

EFBC: ([a.sub.w] + [b.sub.w]) = [X.sup.t] ([a.sub.s] + [b.sub.s]) ([N.sub.w] eqs.), MFBC: X ([a.sub.w] - [b.sub.w]) = ([b.sub.s] - [a.sub.s]) ([N.sub.s] eqs.), (11)

where the term X is the following matrix:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

The previous integral is expressed in terms of the electric field, but it is possible to use the magnetic field as well since [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [Y.sub.j] and [Z.sub.i] stand for the mode admittance and the mode impedance, respectively, in Appendix. The relationship between the modal amplitudes obtained solving the previous linear equation system is given by the GSM of the discontinuity:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

where F = 2[([I.sub.S] + X[X.sup.t]).sup.-1] and I is the identity matrix of the corresponding size ([N.sub.w] or [N.sub.s]). The analysis of two or more discontinuities or steps is carried out cascading GSM matrices as in [10, 11, 19].

4. Mode-Matching Applied to Equilateral Triangular Steps

Section 2 presented the basis functions [[phi].sup.S/A.sub.TE,TM] needed to formulate the modal fields in any equilateral triangular waveguide. The modal fields [[??].sup.(g).sub.p] or [[??].sup.(g).sub.p] are required to compute the inner cross product matrix X seen in Section 3, where the Mode-Matching method was introduced. This section addresses the resolution of the surface integrals associated with the inner cross products and aims at providing a suitable formulation to solve discontinuities where one of the waveguide surfaces is completely included into the other one (see Figure 4). Towards that end, the functions given in (4a)-(4d) must be rewritten in a more convenient form based on plane waves:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

The resulting parameters [A.sub.u], [k.sub.xu], and [k.sub.yu] of the preceding equation are given in Table 2. In order to highlight the benefits of the mathematical formulation used in this paper, it is sometimes convenient to define the wavenumber vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and the position vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and then rewrite [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Once the functions of the four mode families are formulated using (14), the surface integrals involved in the Mode-Matching method are greatly simplified. The first surface integral that must be solved is the one associated with the normalization constant [N.sub.p] (see Appendix). This surface integral is solved in [17, 18] and its values are directly given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

The other surface integrals that must be solved in the Mode-Matching method are those needed to calculate the inner cross product matrix (X). The electric fields in (12) are computed according to Appendix, using the transversal gradient of the function in (14). Each term [X.sub.ij] of the inner cross product matrix X is computed as

[X.sub.ij] = ([X.sub.ij,x] + [X.sub.ij,y]) ([square root of ([Y.sup.(w).sub.j])]/[square root of ([Y.sup.(s).sub.i])], (16)

where the terms [X.sub.ij,x] and [X.sub.ij,y] are independent of the frequency and the material filling the waveguide and refer to the surface integrals associated with the transversal fields in [??] and [??]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

In the above expressions [A.sub.s] is the intersected area between both equilateral waveguides as it is shown in Figure 4. The variables [k.sub.xr], [k.sub.xq], [k.sub.yr], and [k.sub.yq] are the wavenumbers in (14), for each respective waveguide. Only the new amplitudes [A.sub.sr], [A.sub.sq], [A.sub.tr], and [A.sub.tq] must be computed according to the involved gradients shown in Appendix. The results of these amplitudes are given in Table 4.

These inner cross products can be arranged according to the families to which the modes in (12) belong, with matrix X then becoming a set of submatrices:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)

Some of the surface integrals associated with these submatrices can be discarded since their values are zero [22]. Therefore (18) becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

Moreover, having the TE and TM modes divided into symmetries allows discarding even more surface integrals when both waveguides have their symmetry axis aligned as in Figure 4. In that case, symmetric TE modes do not excite asymmetric TE modes, and symmetric TM modes do not excite asymmetric TM modes, nor does this happen in TE and TM modes of different perfect walls. Thus, in those conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

One significant advantage of the mathematical formulation used in this paper to describe the problem of the electromagnetic scattering between equilateral triangular waveguides is that any offset introduced in the reference coordinate system results again in plane waves with the same wavenumber vector and just a change in their associated amplitudes. Figure 5 shows an arbitrary offset between two coordinate systems. According to this representation, the basis functions in (14) referring to the new coordinate system are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (21)

which can be rewritten using the same parameters [k.sub.xu] and [k.sub.yu] already calculated in (14) and therefore introducing a simple change in its amplitudes:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (22)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Another advantage of this formulation is that all the surface integrals involved in the calculation of the inner cross product matrix X can be reduced to the following simple exponential integral:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (23)

where [B.sub.1] and [B.sub.2] are related to the additions of the corresponding wavenumbers and [A.sub.s] is the inner surface of the triangular intersection.

5. Results and Discussion

In order to prove the efficiency of the formulation presented in this paper, three cases are solved. For each case, the GSM is compared at different frequencies with the results given by the commercial software CST Microwave Studio. It is important to bear in mind that the modal series that describe the electromagnetic field in each waveguide must be truncated to a finite number of modes. Nevertheless, it is not advisable to use a very high number of modes from the beginning since it increases the computational cost and therefore the efficiency of the Mode-Matching technique is compromised. Besides, in this scenario, where two modal series have been matched, the problem of relative convergence arises.

One of the well-known strategies to address the problem is to set a reference frequency [f.sup.ref.sub.c] and select the modes whose cut-off frequency is below it in each waveguide. The reference frequency [f.sup.ref.sub.c] is increased until the results decrease their variation with the number of selected modes. It will be seen that at that point they show a good agreement with the commercial software CST Microwave Studio.

According to the selected stopping criterion based on a cut-off frequency threshold, three pairs of mode relationships are used: 34-10, 106-34, and 214-58, where the former number indicates the number of modes considered in the largest waveguide and the latter the ones considered in the smallest one. The three case studies analyzed are depicted in Figure 6.

5.1. Single Step between Aligned Waveguides. In the first case study a single step between two equilateral triangular waveguides is considered, where both of them are centered at our coordinate reference system (see Figure 6(a)). The side of the waveguide (s) is [square root of 3] mm, which is port 2 in the generalized scattering parameters shown later on, whereas the side of the waveguide (w), which is port 1, is 2[square root of 3] mm. Figures 7 and 8 show the comparison of the two first degenerated modes (T[E.sup.S.sub.10] and T[E.sup.A.sub.1]) using the presented formulation for the Mode-Matching method and the commercial software CST Microwave Studio.

The frequency range (i.e., 117-122 GHz) is selected in order to ensure the propagation of the pair of fundamental modes in each waveguide and guarantee a better comparison with the commercial software.

Figures 7(a) and 7(b) show the magnitude and phase of the reflection parameters of the T[E.sup.A.sub.10] mode, that is, the relationship between the amplitudes of the incident and the reflected T[E.sup.A.sub.10] mode in each waveguide. Figures 7(c) and 7(d) represent the magnitude and phase of the transmission parameters of the T[E.sup.A.sub.10] mode, that is, the relationship between the amplitudes of the same type of mode going through both waveguides.

In this case the energy transfer between the T[E.sup.A.sub.10] mode and the T[E.sup.S.sub.10] mode is not represented, since it is known in advance that this transference is zero due to the dual symmetries that these modes present. For the same reason, thanks to the division of the modes into symmetries the number of modes used is half the number used in the two following cases, leading to a faster convergence.

Figure 8 follows the same scheme with the mode T[E.sup.S.sub.10]. It can be seen that, compared to the T[E.sup.A.sub.10] mode, the T[E.sup.s.sub.10] mode requires more modes in the series to achieve the results given by the numerical method in CST in general. In the T[E.sup.A.sub.10] mode only the [absolute value of [S.sub.11]] parameter and the phase of the transmission parameters require a larger number of modes to converge.

5.2. Single Step between Misaligned Waveguides. In the second case study the previous waveguides are selected again, but a vertical offset of 0.08 mm and a horizontal offset of 0.02 mm are introduced (see Figure 6(b)). Figures 9-11 depict the results obtained using both the Mode-Matching technique and the CST Microwave Studio solver again for the two fundamental degenerate modes.

Figures 9 and 10 show that in this case both the T[E.sup.A.sub.10] and the T[E.sup.S.sub.10] modes require a higher number of modes to be similar to the results obtained by the CST commercial software.

It is interesting to include the energy transfer between the modes of different families, which is nonzero in this case study since the symmetry axis between the considered waveguides is no longer aligned. Figure 11 represents the relationships between the modal amplitudes of mode T[E.sup.A.sub.10] in one of the waveguides and mode T[E.sup.S.sub.10] in the other one.

5.3. Double Step between Misaligned Waveguides. In the third and last case study, two steps between the same waveguides are analyzed (see Figure 6(c)). In this case the notation associated with the number of modes used refers to the three waveguides, 34-10-34 being the modes selected in the three waveguides involved. In order to test the evanescent modes (i.e., those whose cut-off frequency is below the analyzed range of frequencies), a range of frequencies is selected where none of the modes in the smallest waveguide propagates (i.e., 90-95 GHz).

The thickness of waveguide (s) is 0.3 mm to limit the attenuation of its first three pairs of modes to 3.9 dB, 10 dB, and 11.5 dB, respectively, having significant higher-order mode interaction between the two discontinuities. In Figures 12-14 it can be seen that all the results achieved using a very low number of modes (i.e., 34-10-34) and those with a higher number show greater differences than in the previous two cases. In any case, the results finally converge once the number of modes is increased.

6. Conclusion

The complete mathematical formulation of the Mode-Matching method applied to equilateral triangular discontinuities has been presented. The mode spectrum for the equilateral triangular waveguide has been shown. A flexible formulation based on plane waves has been presented, simplifying the surface integrals involved. Moreover, this mathematical formulation can be easily extended to discontinuities that involve triangular equilateral waveguides with other classical waveguides such as the rectangular, circular, or elliptical waveguides, where their modes can be expanded in plane waves as well. Finally, three case studies have been analyzed and compared with the results obtained with the commercial software CST Microwave Studio. A good agreement between responses has been proved, with Mode-Matching having the advantage of a quasi-analytical method with high efficiency.

Appendix

Mode Calculation

The electromagnetic field of a TE mode, at the angular frequency [omega] = 2[pi]f, can be expressed using the following equations [13, 14]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.1a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.1b)

where the wave impedance (admittance) [Z.sub.p] ([Y.sub.p]) is

[Z.sub.p] = [Y.sup.-1.sub.p] = j[omega][mu]/[[gamma].sub.p], (A.2)

and the normalization constant for unit complex power [N.sub.p] and the propagation constant [[gamma].sub.p] are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.3a)

[[gamma].sup.2.sub.p] = [k.sup.2.sub.cp] - [[omega].sup.2] [mu][EPSILON]. (A.3b)

The electromagnetic field of a TM mode can be expressed using the following equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.4a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.4b)

where the normalization constant for unit power [N.sub.p] and the wave number [k.sub.cp] hold the same definition as in (A.3a) and (A.3b), respectively, and the wave impedance (admittance) [Z.sub.p] ([Y.sub.p]) is

[Z.sub.p] = [Y.sup.-1.sub.p] = [[gamma].sub.p]/j[omega][EPSILON]. (A.5)

http://dx.doi.org/10.1155/2016/2974675

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the Spanish government under Grant TEC2013-47106-C3-1/2-R and the Comunidad de Madrid Program S2013/ICE-3000 (SPADERADARCM).

References

[1] A. E. Grant and J. H. Meadows, Communication Technology Update and Fundamentals, Focal Press, Oxford, UK, 14th edition, 2014.

[2] D. Raychaudhuri and M. Gerla, Emerging Wireless Technologies and the Future Mobile Internet, Cambridge University Press, New York, NY, USA, 2011.

[3] G. A. Vandenbosch and A. Vasylchenko, A practical guide to 3D electromagnetic software tools, InTech, Riieka, Croatia, 2011.

[4] Q. S. Cheng, J. W. Bandler, S. Koziel, M. H. Bakr, and S. Ogurtsov, "The state of the art of microwave CAD: EM-based optimization and modeling," International Journal of RF and Microwave Computer-Aided Engineering, vol. 20, no. 5, pp. 475-491, 2010.

[5] A. A. Oliner, "Historical perspectives on microwave field theory," IEEE Transactions on Microwave Theory and Techniques, vol. 32, no. 9, pp. 1022-1045, 1984.

[6] A. Wexler, "Solution of waveguide discontinuities by modal analysis," IEEE Transactions on Microwave Theory and Techniques, vol. 15, no. 9, pp. 508-517, 1967.

[7] H. Patzelt and F. Arndt, "Double-plane steps in rectangular waveguides and their application for transformers, irises, and filters," IEEE Transactions on Microwave Theory and Techniques, vol. 30, no. 5, pp. 771-776, 1982.

[8] J. A. Ruiz-Cruz, M. A. El Sabbagh, K. A. Zaki, J. M. Rebollar, and Y. Zhang, "Canonical ridge waveguide filters in LTCC or metallic resonators," IEEE Transactions on Microwave Theory and Techniques, vol. 53, no. 1, pp. 174-182, 2005.

[9] J. Esteban and J. M. Rebollar, "Generalized scattering matrix of generalized two-port discontinuities: application to four-port and nonsymmetric six-port couplers," IEEE Transactions on Microwave Theory and Techniques, vol. 39, no. 10, pp. 1725-1734, 1991.

[10] G. Conciauro, M. Guglielmi, and R. Sorrentino, Advanced Modal Analysis: CAD Techniques for Waveguide Components and Filters, John Wiley & Sons, 1999.

[11] J. A. Ruiz-Cruz, J. M. Rebollar, and J. R. Montejo-Garai, Computer Aided Design of Waveguide Devices by Mode-Matching Methods, InTech, Rijeka, Croatia, 2010.

[12] P. M. Morse and H. Feshbach, Methods of Theoretical Physics, McGraw-Hill, New York, NY, USA, 1953.

[13] R. E. Collin, Field Theory of Guided Waves, IEEE Press, New York, NY, USA, 1991.

[14] D. M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.

[15] S. A. Schelkunoff, Electromagnetic Waves, van Nostrand, New York, NY, USA, 1943.

[16] P. L. Overfelt and D. J. White, "TE and TM modes of some triangular cross-section waveguides using superposition of plane waves," IEEE Transactions on Microwave Theory and Techniques, vol. 34, no. 1, pp. 161-167, 1986.

[17] B. J. McCartin, "Eigenstructure of the equilateral triangle, Part I: the Dirichlet problem," SIAM Review, vol. 45, no. 2, pp. 267-287, 2003.

[18] B. J. McCartin, "Eigenstructure of the equilateral triangle, part II: the Neumann problem," Mathematical Problems in Engineering, vol. 8, no. 6, pp. 517-539, 2002.

[19] T. Itoh, Numerical Techniques for Microwave and Millimeter-Wave Passive Structures, John Wiley & Sons, New York, NY, USA, 1989.

[20] J. Corcoles, M. A. Gonzalez De Aza, and J. Zapata, "Full-wave analysis of finite periodic cylindrical conformal arrays with Floquet spherical modes and a hybrid finite element -generalized scattering matrix method," Journal of Electromagnetic Waves and Applications, vol. 28, no. 1, pp. 102-111, 2014.

[21] J. Rubio, J. Corcoles, and M. Gonzalez de Aza, "Inclusion of the feeding network effects in the generalized-scattering-matrix formulation of a finite array," IEEE Antennas and Wireless Propagation Letters, vol. 8, pp. 819-822, 2009.

[22] G. Figlia and G. G. Gentili, "On the line-integral formulation of mode-matchingtechnique," IEEE Transactions on Microwave Theory and Techniques, vol. 50, no. 2, pp. 578-580, 2002.

Ana Moran-Lopez, (1) Juan Corcoles, (1) Jorge A. Ruiz-Cruz, (1) Jose R. Montejo-Garai, (2) and Jesus M. Rebollar (2)

(1) Escuela Politecnica Superior, Universidad Autonoma de Madrid, C/Francisco Tomasy Valientell, 28049 Madrid, Spain

(2) ETSI de Telecomunicacion, Universidad Politecnica de Madrid, Avenida Complutense 30, 28040 Madrid, Spain

Correspondence should be addressed to Ana Moran-Lopez; ana.moran@uam.es

Received 18 May 2016; Accepted 19 June 2016

Academic Editor: Remi Leandre

Caption: Figure 1: (a) Classical circular waveguide, (b) classical rectangular waveguide, (c) equilateral triangular waveguide of this paper, and (d) reference coordinate system used. Conductivity [sigma] = [infinity] indicates a perfect electric wall boundary condition.

Caption: Figure 2: Electric and magnetic field pattern (in blue and green arrows, resp.) for the first twelve modes in the equilateral triangular waveguide.

Caption: Figure 3: Generic waveguide step: (a) surface views and (b) modal amplitudes.

Caption: Figure 4: Step between equilateral triangular waveguides.

Caption: Figure 5: Change in the coordinate system.

Caption: Figure 6: Case studies: single step (a) without offset and (b) with offset and (c) double step with offset.

Caption: Figure 7: Case study a: mode T[E.sup.A.sub.10] in waveguide 1 to mode T[E.sup.A.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 8: Case study a: mode T[E.sup.S.sub.10] in waveguide 1 to mode T[E.sup.S.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 9: Case study b: mode T[E.sup.A.sub.10] in waveguide 1 to mode T[E.sup.A.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 10: Case study b: mode T[E.sup.S.sub.10] in waveguide 1 to mode T[E.sup.S.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 11: Case study b: mode T[E.sup.A.sub.10] in waveguide 1 to mode T[E.sup.S.sub.10] in waveguide 2 and mode T[E.sup.S.sub.10] in waveguide 1 to mode T[E.sup.A.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 12: Case study c: mode T[E.sup.A.sub.10] in waveguide 1 to mode T[E.sup.A.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 13: Case study c: mode T[E.sup.S.sub.10] in waveguide 1 to mode T[E.sup.S.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Caption: Figure 14: Case study c: mode T[E.sup.A.sub.10] in waveguide 1 to mode T[E.sup.S.sub.10] in waveguide 2 and mode T[E.sup.S.sub.10] in waveguide 1 to mode T[E.sup.A.sub.10] in waveguide 2: (a) reflection magnitude, (b) reflection phase, (c) transmission magnitude, and (d) transmission phase.

Table 1: Cut-off wavenumbers, cut-off frequencies, and associated families for the first fifty modes in a hollow equilateral triangular waveguide of side e = 1 cm. (m, n) [k.sub.c] [f.sub.c] Mode designation [[cm.sup.-1]] [GHz] (1, 0) 4.189 19.99 [TE.sup.A], [TE.sup.S] (1, 1) 7.255 34.62 [TE.sup.S], [TM.sup.S] (2, 0) 8.378 39.97 [TE.sup.A], [TE.sup.S] (2, 1) 11.082 52.88 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (3, 0) 12.566 59.96 [TE.sup.A], [TE.sup.S] (2, 2) 14.510 69.23 [TE.sup.S], [TM.sup.S] (3, 1) 15.103 72.06 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (4, 0) 16.755 79.94 [TE.sup.A], [TE.sup.S] (3, 2) 18.259 87.12 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (4, 1) 19.195 91.59 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (5, 0) 20.944 99.93 [TE.sup.A], [TE.sup.S] (3, 3) 21.766 103.85 [TE.sup.S], [TM.sup.S] (4, 2) 22.165 105.76 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (5, 1) 23.322 111.28 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (6, 0) 25.133 119.92 [TE.sup.A], [TE.sup.S] (4, 3) 25.479 121.57 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] (5, 2) 26.159 124.81 [TE.sup.A], [TE.sup.S], [TM.sup.A], [TM.sup.S] Table 2: Values of parameters [k.sub.xu], [k.sub.yu], and [A.sub.u], where the parameters denoted as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are the ones described immediately after (4a)-(4d), and Y and O are defined in Table 3. u [k.sub.xu] [k.sub.yu] [A.sub.u] 1 [[??].sub.x1] [[??].sub.y1] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 2 -[[??].sub.x1] [[??].sub.y1] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 3 [[??].sub.x1] -[[??].sub.y1] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 4 -[[??].sub.x1] [[??].sub.y1] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 5 [[??].sub.x2] [[??].sub.y2] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 6 -[[??].sub.x2] [[??].sub.y2] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 7 [[??].sub.x2] -[[??].sub.y2] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 8 -[[??].sub.x2] [[??].sub.y2] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 9 [[??].sub.x3] [[??].sub.y3] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 10 -[[??].sub.x3] [[??].sub.y3] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 11 [[??].sub.x3] -[[??].sub.y3] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 12 -[[??].sub.x3] [[??].sub.y3] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Table 3: Values of parameters [PSI] and [OMEGA] for Table 2. Mode designation [PSI] [OMEGA] [TE.sup.S] 0 0 [TE.sup.A] 0 [pi]/2 [TM.sup.S] [pi]/2 0 [TM.sup.A] [pi]/2 [pi]/2 Table 4: Values of parameters A sr, A , Atr, and Atq. r, q [A.sub.sr,sq] [TE] [A.sub.sr,sq] [TM] 1 j[[??].sub.y1][A.sub.u] j[[??].sub.x1][A.sub.u] 2 j[[??].sub.y1][A.sub.u] -j[[??].sub.x1][A.sub.u] 3 -j[[??].sub.y1][A.sub.u] j[[??].sub.x1][A.sub.u] 4 -j[[??].sub.y1][A.sub.u] -j[[??].sub.x1][A.sub.u] 5 j[[??].sub.y2][A.sub.u] j[[??].sub.x2][A.sub.u] 6 j[[??].sub.y2][A.sub.u] -j[[??].sub.x2][A.sub.u] 7 -j[[??].sub.y2][A.sub.u] j[[??].sub.x2][A.sub.u] 8 -j[[??].sub.y2][A.sub.u] -j[[??].sub.x2][A.sub.u] 9 j[[??].sub.y3][A.sub.u] j[[??].sub.x3][A.sub.u] 10 j[[??].sub.y3][A.sub.u] -j[[??].sub.x3][A.sub.u] 11 -j[[??].sub.y3][A.sub.u] j[[??].sub.x3][A.sub.u] 12 -j[[??].sub.y3][A.sub.u] -j[[??].sub.x3][A.sub.u] r, q [A.sub.tr,tq] [TE] [A.sub.tr,tq] [TM] 1 -j[[??].sub.x1][A.sub.u] j[[??].sub.y1][A.sub.u] 2 j[[??].sub.x1][A.sub.u] j[[??].sub.y1][A.sub.u] 3 -j[[??].sub.x1][A.sub.u] -j[[??].sub.y1][A.sub.u] 4 j[[??].sub.x1][A.sub.u] -j[[??].sub.y1][A.sub.u] 5 -j[[??].sub.x2][A.sub.u] j[[??].sub.y2][A.sub.u] 6 j[[??].sub.x2][A.sub.u] j[[??].sub.y2][A.sub.u] 7 -j[[??].sub.x2][A.sub.u] -j[[??].sub.y2][A.sub.u] 8 j[[??].sub.x2][A.sub.u] -j[[??].sub.y2][A.sub.u] 9 -j[[??].sub.x3][A.sub.u] j[[??].sub.y3][A.sub.u] 10 j[[??].sub.x3][A.sub.u] j[[??].sub.y3][A.sub.u] 11 -j[[??].sub.x3][A.sub.u] -j[[??].sub.y3][A.sub.u] 12 j[[??].sub.x3][A.sub.u] -j[[??].sub.y3][A.sub.u]

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Moran-Lopez, Ana; Corcoles, Juan; Ruiz-Cruz, Jorge A.; Montejo-Garai, Jose R.; Rebollar, Jesus M. |

Publication: | Advances in Mathematical Physics |

Article Type: | Report |

Date: | Jan 1, 2016 |

Words: | 6140 |

Previous Article: | Inverse uniqueness in interior transmission problem and its eigenvalue tunneling in simple domain. |

Next Article: | A nonlinear Schrodinger equation resonating at an essential spectrum. |

Topics: |