# Complex Resonances of a Rectangular Patch in a Multilayered Medium: A New Accurate and Efficient Analytical Technique.

1. INTRODUCTIONFor several decades, the study of microwave antennas and resonators based on microstrip technology has been an important issue. This can be explained by their conformability and minimized shape, low weight and cost, and due to the high number of applications including, just for examples, mobile satellite communications, direct broadcast satellite services and non-destructive testing sensors in permittivity and porosity measurements [1-6].

Furthermore, the use of multilayered dielectric substrates allows to obtain many advantages as evidenced by the many works devoted to this subject [7-19]. It has been verified that multilayered dielectric substrates can be used to reduce the radiation losses by surface waves to enhance the efficiency. On the other hand, a dielectric substrate placed on top of a microstrip antenna can provide protection against environmental hazards and, appropriately located, can increase the gain of the antenna. Moreover, multilayered microstrip circuits allow for more versatile designs and offer the advantage of greater compactness.

What emerges from the literature devoted to analyze microstrip circuits and antennas in multilayered media, is that full-wave methods [20-22], taking into account the effects of the electromagnetic coupling, surface waves, and radiation loss, are more accurate than traditional quasi-static methods [23] and equivalent waveguide models [24]. Moreover, the most preferred full-wave techniques to analyze non-shielded structures are integral equation formulations, allowing to express the fields as functions of unknowns on finite support, discretized by means of the variational method of moments.

Among them, the method introduced in [14, 25-39] to analyse propagation, radiation and scattering by polygonal cross-section cylinders and one-dimensional perfectly conducting plates (disks, annular rings, ...) in a homogeneous or a layered medium has been shown to be a very accurate and efficient approach. The problem is formulated in terms of surface integral equations in the spectral domain and discretized by means of Galerkin's method with analytically Fourier transformable orthonormal expansion basis reconstructing the behaviour of the fields at the edges [40]. Moreover, suitable analytical acceleration techniques to fast evaluate the obtained improper integrals of oscillating and slowly decaying functions has been developed. Even if this method can be immediately extended to the analysis of the radiation and the scattering by a perfectly conducting rectangular plate in a homogeneous or a layered medium, the proposed acceleration techniques reveals to be not effective when applied to the resulting infinite double integrals of oscillating and slowly decaying functions. In order to overcome this problem, in [41] the range of integration has been subdivided into annular sectors and the region external to the maximum sector divided into three subregions over which the generic integrand has been approximated with a suitable asymptotic behaviour so to reduce the corresponding double integral to a single one. Recently, a new technique, drastically outperforming the acceleration technique presented in [41], has been introduced in [42]: a new quickly converging representation of the infinite double integrals resulting from the analysis of the scattering by a rectangular plate in a homogeneous medium has been obtained by means of a suitable integration procedure in the complex plane. Unfortunately, in case of a multilayered medium, such technique cannot be immediately applied due to the oscillating nature of the Green's functions of the problem.

The aim of this paper is the introduction of a new accurate and efficient analytical technique to study the complex resonances of a rectangular patch in a multilayered medium. The problem formulated as an EFIE in the spectral domain is discretized by means of products of Chebyshev polynomials of first and second kind multiplied by their orthogonal weights in a Galerkin's scheme. In order to fast evaluate the obtained infinite double integrals of oscillating and slowly decaying functions, suitable half-space contributions are pulled out of the kernels of such integrals obtaining exponentially decaying integrands. Moreover, the slowly converging integrals of the extracted contributions are expressed as combinations of proper and fast converging improper integrals of non-oscillating and exponentially decaying functions by means of suitable algebraic manipulations and using an appropriate integration procedure in the complex plane generalizing the one proposed in [42].

In Sections 2 the formulation of the problem and the discretization of the integral equations are presented. The new acceleration technique is illustrated in Section 3. Section 4 is devoted to show the accuracy and the efficiency of the presented technique and the conclusions are summarized in Section 5.

2. FORMULATION AND SOLUTION OF THE PROBLEM

The geometry of the problem is sketched in Figure 1. A planar layered medium of L+1 homogeneous and isotropic layers of dielectric permittivity [[epsilon].sub.q] = [[epsilon].sub.0] [[epsilon].sub.rq, magnetic permeability [[mu].sub.q] = [[mu].sub.0][[mu].sub.rq] and wave number [k.sub.q] = [omega] [squarer root of [[epsilon].sub.q][[mu].sub.q]] with q [member of] {1, ..., L + 1} is represented, where [[epsilon].sub.0] and [[mu].sub.0] are the dielectric permittivity and the magnetic permeability of the vacuum and [omega] is the angular frequency. Moreover, a coordinate system with the z axis orthogonal to the discontinuity surfaces is introduced, and at the interface z = -[d.sub.q], a perfectly conducting rectangular patch of dimensions 2a and 2b, with the sides parallel to the x and y axes and the centre at the point x = y = 0, is located.

Resonant solutions can be obtained by imposing the transversal component of the electric field with respect to the z axis to be vanishing on the patch's surface, i.e., [43]

[mathematical expression not reproducible] (1)

with [absolute value of x] [less than or equal to] a and [absolute value of y] [less than or equal to] b, where [[??].bar](*) is the double Fourier transform with respect to the x and y axes of the surface current density on the rectangular patch and [44, 45]

[mathematical expression not reproducible], (2a)

[mathematical expression not reproducible], (2b)

[mathematical expression not reproducible], (2c)

[Z.sup.TE.sub.q] = [omega][[mu].sub.q]/[k.sub.q,z], (2d)

[Z.sup.TM.sub.q] = [k.sub.q,z]/[omega][[epsilon].sub.q], (2e)

[[??].sup.TE,TM.sub.q[+ or -]] = [[??].sup.TE,TM.sub.L+1,+] = 0, (2f)

[mathematical expression not reproducible], (2g)

[R.sup.TE,TM.sub.1,-] = [R.sup.TE,TM.sub.L+1,+] = 0, (2h)

[R.sup.TE.sub.q,[+ or -]] = [[mu].sub.q[+ or -]1] [k.sub.q,z] - [[mu].sub.q][k.sub.q[+ or -]1,z]/ [[mu].sub.q[+ or -]1] [k.sub.q,z] + [[mu].sub.q][k.sub.q[+ or -]1,z], (2i)

[R.sup.TM.sub.q,[+ or -]] = [[epsilon].sub.q[+ or -]1] [k.sub.q,z] - [[epsilon].sub.q][k.sub.q[+ or -]1,z]/ [[epsilon].sub.q[+ or -]1] [k.sub.q,z] + [[epsilon].sub.q][k.sub.q[+ or -]1,z], (2j)

[mathematical expression not reproducible], (2k)

[k.sub.q,z] = [square root of ([k.sup.2.sub.q] - [k.sup.2.sub.y] - [k.sup.2.sub.x])]. (2l)

The system of integral Equation (1) can be reduced to a symmetric matrix equation by means of Galerkin's method with products of Chebyshev polynomials of first and second kind multiplied by their orthogonal weights, i.e., factorizing the behaviour of the surface current density at the edges [40], as expansion functions. Therefore, the elements of the coefficient matrix are proportional to the following double integrals

[mathematical expression not reproducible] (3)

with r, s [member of] {x, y} and n, h, m, k nonnegative integers, where

[mathematical expression not reproducible], (4a)

[mathematical expression not reproducible] (4b)

is proportional to the Fourier transform of the generic expansion function, [p.sub.x] = 1, [p.sub.y] = 0, and [J.sub.v](*) is the Bessel function of first kind and order v. Moreover, the complex resonant frequencies of the problem can be readily obtained by enforcing the determinant of the truncated coefficient matrix to be zero.

3. A NEW ACCELERATION TECHNIQUE TO EFFICIENTLY EVALUATE THE ELEMENTS OF THE COEFFICIENTS MATRIX

As a first task, by means of the change of variables [k.sub.x] = [rho][c.sub.[psi]] and [k.sub.y] = [rho][s.sub.[psi]], where [c.sub.[psi]] = cos [psi] and [s.sub.[psi]] = sin [psi], it is simple to rewrite (3) as follows

[mathematical expression not reproducible]. (5)

Observing that [46]

[mathematical expression not reproducible], (6)

it can be concluded that the integrand of the integral in (5) is an oscillating function which decays asymptotically as 1/[rho].sub.2]. Therefore, the computation time rapidly increases with the accuracy required for the solution.

In order to overcome this problem, a new acceleration technique is proposed here.

The expression in (2a) can be reduced to the following non-oscillating function by posing ([[epsilon].sub.l], [[mu].sub.l]) = ([[epsilon].sub.g], [[mu].sub.q]) for l [member of] {1, ..., q-1} and ([[epsilon].sub.l], [[mu].sub.l]) = ([[epsilon].sub.g+1], [[mu].sub.q+1]) for l [member of] {q + 2, ..., L + 1}, i.e., when only two half-spaces are involved,

[mathematical expression not reproducible]. (7)

It is simple to verify that the difference between the functions in (2a) and (7) is

[mathematical expression not reproducible]. (8)

where

[mathematical expression not reproducible], (9a)

[mathematical expression not reproducible], (9b)

and the upper sign has to be taken for the TE case whilst the lower sign has to be taken for the TM case.

Therefore, by rewriting the double integral in (5) as

[mathematical expression not reproducible] (10)

and remembering (6), it can be concluded that the integrand of the second double integral (integral of the accelerated contribution) has an exponential asymptotic decay of the kind [mathematical expression not reproducible] where

[mathematical expression not reproducible]. (11)

On the other hand, as can be shown to follow, the slowly converging first double integral in (10) (integral of the extracted contribution) can be rewritten as a combination of very quickly converging integrals taking advantage of the non-oscillating nature of the kernel.

Firstly, by means of algebraic manipulations, the recurrence formula [46]

2v[J.sub.v](z) = z [[J.sub.v-1](z) + [J.sub.v+1(z)], (12)

and the property [46]

[J.sub.v] (z[e.sup.jp[pi]]) = [e.sup.jpv[pi]] [J.sub.v](z) (13)

with p integer, it is possible to rewrite the first integral in (10) as follows

[mathematical expression not reproducible], (14a)

[mathematical expression not reproducible], (14b)

[mathematical expression not reproducible], (14c)

where

[mathematical expression not reproducible], (15a)

[mathematical expression not reproducible], (15b)

[mathematical expression not reproducible] (15c)

with [eta] [member of] {[epsilon], [mu]}. It is immediate to establish that [mathematical expression not reproducible] that allow to consider only the cases for n + m and h + k even, with n [greater than or equal to] m and h [greater than or equal to] k (henceforth, these assumptions will be implicitly done).

Let us focus on the single integral over the semi-infinite range in (15a) that can be rewritten as follows

[mathematical expression not reproducible], (16)

where the choice of d > max{[k.sub.q], [k.sub.q+1]} has been discussed in [42]. The first integral at the second member of the previous formula is a proper integral. Moreover, supposing that a[c.sub.[psi]] [greater than or equal to] b[s.sub.[psi]] (analogous considerations can be done for a[c.sub.[psi]] < b[s.sub.[psi]] by reversing the role of the sides), the following alternative expression for the second integral in (16) can be found by generalizing the integration procedure in the complex plane detailed in [42]

[mathematical expression not reproducible], (17)

where

[mathematical expression not reproducible], (18a)

[h.sub.n,h,m,k] ([alpha], [beta]) = ([h.sup.(1).sub.n,h,m,k] ([alpha], [beta]) + [h.sup.(2).sub.n,m,k] ([alpha], [beta])/2, (18b)

[h.sup.(1).sub.n,m,k] ([alpha], [beta]) = [J.sub.n] ([alpha]) [J.sub.h] ([beta]) [H.sup.(1).sub.m] ([alpha]) [H.sup.(1).sub.k] ([beta]), (18c)

[mathematical expression not reproducible], (18d)

[I.sub.v](*) and [K.sub.v](*) are the modified Bessel functions of first and second kind and order v respectively, [Y.sub.v](*) is the Bessel function of second kind and order v, p[H.sup.(1).sub.v](*) = [J.sub.v](*) + j[Y.sub.v](*) is the Hankel function of first kind and order v, and R{*} denotes the real part of a complex number. Observing that [46]

[mathematical expression not reproducible], (19a)

[mathematical expression not reproducible], (19b)

it is possible to conclude that the first integral at the second member of (17) is an improper integral of a non-oscillating function which decays asymptotically as [mathematical expression not reproducible], while the last one is a proper integral.

Therefore, the integral at the first member of (16) (and, then, the double integral in (15a)) has been expressed as a combination of quickly converging integrals.

4. NUMERICAL RESULTS

The aim of this section is to show the accuracy and the efficiency of the presented technique in terms of calculation time and storage requirement.

It is worth noting that more than 20 integrals per second can be computed by means of an adaptive Gauss-Legendre cubature routine on a laptop equipped with an Intel Core 2 Duo CPU T9600 2.8 GHz, 3 GB RAM, running Windows XP.

The overall number of matrix coefficients is 4[N.sup.2][M.sup.2], where N and M are the number of expansion functions used along the x and y axes respectively for each component of the surface current density. However, due to the properties detailed above, the number of integrals of the form [mathematical expression not reproducible] that has to be computed is respectively reduced to NM(2NM + 1), [N.sup.+][N.sup.-][M.sup.+][M.sup.-]/16 and [N.sup.+]([N.sup.-] + 2)[M.sup.+]([M.sup.-] + 2)/16, where [P.sup.[+ or -]] = P + 2 [+ or -] mod(P, 2) with P [member of] {N, M} and mod(*,*) is the modulus operation.

Moreover, few expansion functions are needed to achieve accurate results. It can be appreciated in Table 1 where the complex resonant frequency of the square patch in a grounded one-layer dielectric medium sketched in Figure 2(a) is shown as a function of the number of expansion functions used supposing that N = M. It is clear as the accuracy on the real part of the resonant frequency reaches 3, 4, 5 and 6 significant figures by using, respectively, only 2, 7, 12 and 14 expansion functions along the x and y axes for each current component.

To conclude, in Tables 2 and 3, the resonant frequencies of the square patches in grounded two-layer and three-layer dielectric media sketched in Figures 2(b) and 2(c) are shown. The results obtained with this method (by using only 6 expansion functions along the x and y axes for each current component) and by means of a suitable application of the conformal mapping technique [18] are compared with the experimental data and the results achieved by using the commercial software HFSS. As can be seen, the presented method results to be more accurate than the formulas in [18]. Indeed, the resonant frequencies obtained with the presented method differ always less than 1% from the ones obtained by using HFSS, and only in two cases differ more than 2% and 3% from the experimental data (although it should be noted that, in such cases, the agreement with HFSS is excellent). Conversely, the percentage deviation of the resonant frequencies obtained with the formulas in [18] from those provided by HFSS and the experimental data is strictly dependent on the dielectric permittivities of the layers (just for an example, the percentage deviation showed in Table 2 increases more and more as [[epsilon].sub.r2] increases).

5. CONCLUSION

In this paper, a new accurate and efficient technique to analyze the complex resonances of a rectangular patch in a multilayered medium has been introduced. Future perspective is the generalization of the method to cases in which anisotropic media are involved.

ACKNOWLEDGMENT

This work was supported in part by the Italian Ministry of University (MIUR) under a Program for the Development of Research of National Interest (PRIN Grant 20093CJEJ5_002).

REFERENCES

[1.] Carver, K. R. and J. W. Mink, "Microstrip antenna technology," IEEE Trans. Antennas Propag., Vol. 29, No. 1, 2-24, Jan. 1981.

[2.] James, J. R. and P. S. Hall, Handbook of Microstrip Antennas, Stevenage, Peregrinus, UK, 1989.

[3.] Michalski, K. A. and D. Zheng, "Analysis of microstrip resonators of arbitrary shape," IEEE Trans. Microw. Theory Tech., Vol. 40, No. 1, 112-119, Jan. 1992.

[4.] Gupta, K. C. and M. D. Abouzahra (eds.), Analysis and Design of Planar Microwave Components, IEEE Press, Piscataway, NJ, 1994.

[5.] Bogosanovich, M., "Microstrip patch sensor for measurement of the permittivity of homogeneous dielectric materials," IEEE Trans. Instrum. Meas., Vol. 49, No. 5, 1144-1148, Oct. 2000.

[6.] Zucchelli, A., M. Chimenti, and E. Bozzi, "Application of a coaxial-fed patch to microwave nondestructive porosity measurements in low-loss dielectrics," Progress In Electromagnetics Research M, Vol. 5, 1-14, 2008.

[7.] Bahl, J., P. Bahartia, and S. S. Stuchly, "Design of microstrip antennas covered with a dielectric layer," IEEE Trans. Antennas Propag., Vol. 30, No. 2, 314-318, Mar. 1982.

[8.] Bhattacharayya, A. and T. Tralman, "Effects of dielectric superstrate on patch antennas," Electron. Lett., Vol. 24, No. 6, 356-358, 1988.

[9.] Harokopus, W. P. and P. B. Katehi, "Characterization of microstrip discontinuities on multilayer dielectric substrates including radiation losses," IEEE Trans. Microw. Theory Tech., Vol. 37, No. 12, 2058-2066, Dec. 1989.

[10.] Schwab, W. and W. Menzel, "On the design of planar microwave components using multilayer structures," IEEE Trans. Microw. Theory Tech., Vol. 40, No. 1, 67-72, Jan. 11992.

[11.] Afzalzadeh, R. and R. N. Karekar, "X-band directive single microstrip patch antenna using dielectric parasite," Electron. Lett., Vol. 28, No. 1, 17-19, 1992.

[12.] Yeung, E. K. L., J. C. Beal, and Y. M. M. Antar, "Multilayer microstrip structure analysis with matched load simulation," IEEE Trans. Microw. Theory Tech., Vol. 43, No. 1, 143-149, Jan. 1995.

[13.] Tsai, M.-J., F. De Flaviis, O. Fordham, and N. G. Alexopoulos, "Modeling planar arbitrarily shaped microstrip elements in multilayered media," IEEE Trans. Microw. Theory Tech., Vol. 45, No. 3, 330-337, Mar. 1997.

[14.] Losada, V., R. R. Boix, and M. Horno, "Resonant modes of circular microstrip patches in multilayered substrates," IEEE Trans. Microw. Theory Tech., Vol. 47, No. 4, 488-498, Apr. 1999.

[15.] Ling, F., D. Jiao, and J.-M. Jin, "Efficient electromagnetic modeling of microstrip structures in multilayer media," IEEE Trans. Microw. Theory Tech., Vol. 47, No. 9, 1810-1818, Sep. 1999.

[16.] Xia, L., C.-F. Wang, L.-W. Li, P.-S. Kooi, and M.-S. Leong, "Fast characterization of microstrip antenna resonance in multilayered media using interpolation/extrapolation methods," Microw. Opt. Technol. Lett., Vol. 28, No. 5, 342-346, Mar. 2001.

[17.] Sharma, A. and G. Singh, "Design of single pin shorted three-dielectric-layered substrates rectangular patch microstrip antenna for communication systems," Progress In Electromagnetics Research Letters, Vol. 2, 157-165, 2008.

[18.] Li, Y. and N. Bowler, "Resonant frequency of a rectangular patch sensor covered with multilayered dielectric structures," IEEE Trans. Antennas Propag., Vol. 58, No. 6, 1883-1889, Jun. 2010.

[19.] Aouabdia, N., N. E. Belhadj-Tahar, G. Alquie, and F. Benabdelaziz, "Theoretical and experimental evaluation of superstrate effect on rectangular patch resonator parameters," Progress In Electromagnetics Research B, Vol. 32, 129-147, 2011.

[20.] Taflove, A., Computional Electromagentics: The Finite-difference Time-domain Method, Artech House, Norwood, MA, 1995.

[21.] Jin, J. M., The Finite Element Method in Electromagnetics, Wiley, New York, 1993.

[22.] Harrington, R. F., Field Computation in Electromagnetics, Wiley, New York, 1993.

[23.] Silvester, P. and P. Benedek, "Equivalent capacitance of microstrip open circuits," IEEE Trans. Microw. Theory Tech., Vol. 20, No. 8, 511-516, Aug. 1972.

[24.] Kompa, G. and R. Mehran, "Planar waveguide model for calculating microstrip components." Electron. Lett., Vol. 11, No. 19, 459-460, 1975.

[25.] Eswaran, K., "On the solutions of a class of dual integral equations occurring in diffraction problems," Proceedings of the Royal Society of London, Series A, Vol. 429, No. 1877, 399-427, 1990.

[26.] Bagby, J. S., C.-H. Lee, Y. Yuan, and D. P. Nyquist, "Entire-domain basis MOM analysis of coupled microstrip transmission lines," IEEE Trans. Microw. Theory Tech., Vol. 40, No. 1, 49-57, Jan. 1992.

[27.] Veliev, E. I. and V. V. Veremey, "Numerical-analytical approach for the solution to the wave scattering by polygonal cylinders and flat strip structures," Analytical and Numerical Methods in Electromagnetic Wave Theory, M. Hashimoto, M. Idemen, and O. A. Tretyakov (eds.), Science House, Tokyo, 1993.

[28.] Park, S. and C. A. Balanis, "Dispersion characteristics of open microstrip lines using closed-form asymptotic extraction," IEEE Trans. Microw. Theory Tech., Vol. 45, No. 3, 458-460, Mar. 1997.

[29.] Park, S. and C. A. Balanis, "Closed-form asymptotic extraction method for coupled microstrip lines," IEEE Microw. Guid. Wave Lett., Vol. 7, No. 3, 84-86, Mar. 1997.

[30.] Amari, S., R. Vahldieck, and J. Bornemann, "Using selective asymptotics to accelerate dispersion analysis of microstrip lines," IEEE Trans. Microw. Theory Tech., Vol. 46, No. 7, 1024-1027, Jul. 1998.

[31.] Tsalamengas, J. L., "Rapidly converging direct singular integral-equation techniques in the analysis of open microstrip lines on layered substrates," IEEE Trans. Microw. Theory Tech., Vol. 49, No. 3, 555-559, Mar. 2001.

[32.] Lucido, M., G. Panariello, and F. Schettino, "Analysis of the electromagnetic scattering by perfectly conducting convex polygonal cylinders," IEEE Trans. Antennas Propag., Vol. 54, No. 4, 1223-1231, Apr. 2006.

[33.] Lucido, M., G. Panariello, and F. Schettino, "Electromagnetic scattering by multiple perfectly conducting arbitrary polygonal cylinders," IEEE Trans. Antennas Propag., Vol. 56, No. 2, 425-436, Feb. 2008.

[34.] Lucido, M., G. Panariello, and F. Schettino, "TE scattering by arbitrarily connected conducting strips," IEEE Trans. Antennas Propag., Vol. 57, No. 7, 2212-2216, Jul. 2009.

[35.] Lucido, M., G. Panariello, and F. Schettino, "Scattering by polygonal cross-section dielectric cylinders at oblique incidence," IEEE Trans. Antennas Propag., Vol. 58, No. 2, 540-551, Feb. 2010.

[36.] Coluccini, G., M. Lucido, and G. Panariello, "TM scattering by perfectly conducting polygonal cross-section cylinders: A new surface current density expansion retaining up to the second-order edge behavior," IEEE Trans. Antennas Propag., Vol. 60, No. 1, 407-412, Jan. 2012.

[37.] Lucido, M., "An analytical technique to fast evaluate mutual coupling integrals in spectral domain analysis of multilayered coplanar coupled striplines," Microw. Opt. Technol. Lett., Vol. 54, No. 4, 1035-1039, Apr. 2012.

[38.] Coluccini, G., M. Lucido, and G. Panariello, "Spectral domain analysis of open single and coupled microstrip lines with polygonal cross-section in bound and leaky regimes," IEEE Trans. Microw. Theory Tech., Vol. 61, No. 2, 736-745, Feb. 2013.

[39.] Lucido, M., "An efficient evaluation of the self-contribution integrals in the spectral-domain analysis of multilayered striplines," IEEE Antennas Wireless Propag. Lett., Vol. 12, 360-363, Mar. 2013.

[40.] Meixner, J., "The behaviour of electromagnetic fields at edges," IEEE Trans. Antennas Propag., Vol. 20, No. 4, 442-446, Jul. 1972.

[41.] Hongo, K. and H. Serizawa, "Diffraction of electromagnetic plane wave by rectangular plate and rectangular hole in the conducting plate," IEEE Trans. Antennas Propag., Vol. 47, No. 6, 1029-1041, Jun. 1999.

[42.] Coluccini, G. and M. Lucido, "A new high efficient analysis of the scattering by a perfectly conducting rectangular plate," IEEE Trans. Antennas Propag., Vol. 61, No. 5, 2615-2622, May 2013.

[43.] Chew, W. C. and Q. Liu, "Resonance frequency of a rectangular microstrip patch," IEEE Trans. Antennas Propag., Vol. 36, No. 8, 1045-1056, Aug. 1988.

[44.] Chew, W. C., Waves and Fields in Inhomogeneous Media, IEEE Press, New York, 1995.

[45.] Chew, W. C. and S. Y. Chen, "Response of a point source embedded in a layered medium," IEEE Antennas Wireless Propag. Lett., Vol. 2, No. 1, 254-258, 2003.

[46.] Abramowitz, M. and I. A. Stegun, Handbook of Mathematical Functions, Verlag Harri Deutsch, 1984.

Mario Lucido *

Received 2 February 2014, Accepted 21 February 2014, Scheduled 11 March 2014

* Corresponding author: Mario Lucido (lucido@unicas.it).

The author is with the D.I.E.I., Universita degli Studi di Cassino e del Lazio Meridionale, Cassino 03043, Italy.

Caption: Figure 1. Geometry of the problem.

Caption: Figure 2. Geometries analyzed throughout the paper. Square patches in a grounded (a) one-layer, (b) two-layer, and (c) three layer dielectric media.

Table 1. Complex resonant frequency of the square patch in a grounded one-layer medium sketched in Figure 2(a) for different numbers of basis functions used. 2a = 2b = 10 mm, [[epsilon].sub.r1] = 1, [[epsilon].sub.r2] = 1-046, [d.sub.2]-[d.sub.1] = 0.98 mm. N = M [f.sub.r] (GHz) 1 -- 2 12.50867 + j0.416032 3 12.49978 + j0.441532 4 12.49636 + j0.441182 5 12.49576 + j0.441112 6 12.49522 + j0.441061 7 12.49491 + j0.441032 8 12.49474 + j0.441015 9 12.49463 + j0.441004 10 12.49456 + j0.440997 11 12.49451 + j0.440992 12 12.49447 + j0.440989 13 12.49445 + j0.440987 14 12.49444 + j0.440986 15 12.49444 + j0.440986 Table 2. Resonant frequencies of the square patch in a grounded two-layer dielectric medium sketched in Figure 2(b). 2a = 2b = 10 mm, [[epsilon].sub.r1] = 1, [[epsilon].sub.r3] = 1.046, [d.sub.3]-[d.sub.2] = 0.98 mm. Li et al. [18] [[epsilon].sub.r2] [d.sub.2]- [d.sub.1] [f.sub.r] (GHz) Deviation (%) (mm) HFSS Meas. 1 -- 12.333 -0.54 0.06 3 0.508 11.259 -0.58 -0.58 3.27 0.508 11.195 0.18 -0.04 4.5 0.508 10.921 1.59 -0.26 6 0.635 10.502 4.76 1.96 6.15 0.635 10.471 4.45 3.67 9.2 0.635 9.930 7.35 4.25 9.8 0.635 9.838 7.52 3.83 10.2 0.635 9.780 7.77 7.77 This method [f.sub.r] (GHz) Deviation (%) HFSS Meas. 12.495 (+j0.441) 0.77 1.38 11.356 (+j0.354) 0.27 0.27 11.245 (+j0.347) 0.63 0.40 10.800 (+j0.317) 0.47 -1.37 10.099 (+j0.279) 0.74 -1.95 10.055 (+j0.277) 0.30 -0.45 9.292 (+j0.233) 0.45 -2.45 9.165 (+j0.226) 0.16 -3.27 9.085 (+j0.221) 0.11 0.11 Table 3. Resonant frequencies of the square patch in a grounded three-layer dielectric medium sketched in Figure 2(c). 2a = 2b = 10 mm, [[epsilon].sub.r1] = 1, [[epsilon].sub.r4] = 1.046, [d.sub.4]-[d.sub.3] = 0.99 mm. [[epsilon].sub.r2] [d.sub.2]- [[epsilon].sub.r3] [d.sub.3]- [d.sub.1] [d.sub.2] (mm) (mm) 6 0.635 3.27 0.508 9.2 0.635 3.27 0.508 3.27 0.508 6 0.635 9.2 0.635 6 0.635 3.27 0.508 9.2 0.635 6 0.635 9.2 0.635 Li et al. [18] [f.sub.r] (GHz) Deviation (%) HFSS Meas. 9.963 -3.74 -2.83 9.408 -2.25 -0.44 9.437 -2.46 -2.71 8.494 -3.75 -4.56 8.922 -0.31 -0.87 8.447 -1.2 -1.78 This method [f.sub.r] (GHz) Deviation (%) HFSS Meas. 10.060(+j0.340) -0.1 0.85 9.550(+j0.345) -0.78 1.06 9.687(+j0.273) 0.12 -0.13 8.745(+j0.277) -0.91 -1.74 8.952(+j0.226) 0.02 -0.53 8.519(+j0.224) -0.36 -0.94

Printer friendly Cite/link Email Feedback | |

Author: | Lucido, Mario |
---|---|

Publication: | Progress In Electromagnetics Research |

Date: | Feb 1, 2014 |

Words: | 4580 |

Previous Article: | MIMO Reference Antennas with Controllable Correlations and Total Efficiencies. |

Next Article: | A Novel Array Error Estimation Method for Azimuth Multichannel SAR. |

Topics: |