International Journal of Antennas and PropagationVolume 2012 (2012), Article ID 490327, 9 pagesdoi:10.1155/2012/490327

Research Article

Modal Resonant Frequencies and Radiation Quality Factors of Microstrip Antennas

Jan Eichler, Pavel Hazdra, Miloslav Capek, and Milos Mazanek

Department of Electromagnetic Field, Faculty of Electrical Engineering, Czech Technical University in Prague, Technická 2, 166 27 Prague, Czech Republic

Received 9 August 2011; Revised 10 January 2012; Accepted 13 January 2012

Academic Editor: Charles Bunting

Copyright © 2012 Jan Eichler et al. 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.

The chosen rectangular and fractal microstrip patch antennas above an infinite ground plane are analyzed by the theory of characteristic modes. The resonant frequencies and radiation are evaluated. A novel method by Vandenbosch for rigorous evaluation of the radiation is employed for modal currents on a Rao-Wilton-Glisson (RWG) mesh. It is found that the resonant frequency of a rectangular patch antenna with a dominant mode presents quite complicated behaviour including having a minimum at a specific height. Similarly, as predicted from the simple wire model, the radiation exhibits a minimum too. It is observed that the presence of out-of-phase currents flowing along the patch antenna leads to a significant increase of the factor.

1. Introduction

Evaluation of the basic properties of microstrip patch antennas (MPA) has been numerously discussed in literature, see, for example, [1–3]. The two main MPA attributes are resonant frequency (or frequencies but we will deal mostly with the dominant mode) and the radiation factor. So far, only approximate results and semianalytic equations have been published. To our knowledge, this is the first time that these important characteristics have been studied in a rigorous way. The antennas are treated by using a modal approach (hence we do not a priori consider any feeding to be connected), namely, by the theory of characteristic modes (TCM), [4, 5]. Evaluation of the radiation is performed both by the TCM from the eigenvalues slope and by novel rigorous equations derived by Vandenbosch [6] and Vandenbosch and Volski [7].

2. The Theory of Characteristic Modes

For completeness, let us formulate the basics of the characteristic modes for perfectly conducting bodies of area . The scattered field is related to the electric surface currents by the electric field integral equation (EFIE) [8] Equation (1) is usually treated within the method of moments (MoMs) [8] framework and, due to the structure discretization, the operator is known as the “complex impedance matrix” .

Then the associated Euler’s equation to be solved is Equation (2) is a standard weighted eigenvalue equation leading to a set of real characteristic eigencurrents and associated eigenvalues . Properties of eigenvalues are described in [9], at this moment it is important to note that   reflects the amount of net reactive power (thus means resonance). Instead of eigenvalues, the so-called characteristic angles are introduced to show more visible behavior with frequency [9]. Characteristic currents form a complete orthogonal set, and hence the total current on a conducting body may be expressed as a linear combination of these mode currents [10].

2.1. Implementation of the Characteristic Modes Theory

Implementation of the modal decomposition process has been done in the MATLAB [11] environment using Makarov EFIE codes [12] with the RWG basis functions [13]. This usage is restricted to arbitrary 3D PEC structures with air dielectrics. Our developed TCM tool [14] has the following main advantages:(i)Comsol Multiphysics [15]/MATLAB’s PDE TooLbox mesh import,(ii)Optional Green’s function for infinite ground plane simulations,(iii)Single solver/multicore solver/distributed solver (within a computer network with installed MATLAB).

3. The Radiation Factor

In [6] a novel theory able to rigorously calculate radiated power and stored energies directly from currents flowing along the antenna has been presented. The radiation factor is then readily evaluated by the definition [16]: The equations for radiated power and stored electric and magnetic energies , are where where is a free-space wavenumber, is the surface current density, and is the distance between interacting current elements. The tilde denotes that the radiation contribution has been subtracted from the stored energies at every point in space [17]. It is assumed that the currents are flowing in a vacuum.

3.1. The Modal Radiation Factor

The modal radiation factor may be evaluated from the slope of modal eigenvalues [18]: In [18], (10) is supposed to be an approximation of the radiation , but in resonance it is actually exact.

Since characteristic modes are normalized to radiate unit power [4], (3) reduces to For parallel or series RLC circuit (hence, for one mode), the “impedance ” equals the exact “current ” [6]: Inserting valid for lossless antennas [19] and using the fact that , (12) results in providing that It is therefore concluded that the modal equals the by definition, and it can be proven (using the reactance theorem [20, 21]) that in resonance it also equals the radiation defined from energies by (11).

3.2. Software Implementation

The above equations were implemented in MATLAB for the RWG triangular mesh where two different interaction situations occur:

(a) Distant ElementsWhen the triangular elements are not overlapping, current density on triangles may be simply approximated as point sources located at the centre of triangles [22], see Figure 1. No actual integration is then needed. This centroid approach is very fast with satisfactory accuracy as will be shown later (however it may fail for patches located very close to the ground plane).

Figure 1: Distance between nonoverlapping current elements [23].

(b) Overlapping (Self) ElementsAs known from the method of moments, the so-called “self” contributions are of great importance when dealing with calculations on discrete elements (meshes).Here, the self-interaction occurs when two triangles are overlapping each other. Due to the behavior of integral kernels, only rapidly varying term has to be carefully treated. Since ( being the longest side of the triangle ) is satisfied, one needs only to use the first term in the Taylor series expansion. The dominant singular static part is and the integral to be worked out is where is a triangular area. Using simplex coordinates transformation (Figure 2), the result is [23, 24] where is the triangle area, are the edge lengths (see Figure 2), and is the perimeter of the triangle.

Figure 2: Self-term evaluation. (a) Original problem, (b) simplex coordinates transformation [23].

4. Applications: Rectangular Patch Antenna

Let us first concentrate on a rectangular patch antenna of dimensions  mm and  mm (further noted as R50 × 30) placed in air at a height above an infinite ground plane. Only the dominant TM₀₁ mode will be studied. The reason for choosing a patch with is that we do not have to deal with degenerated modes.

Using the image theory, the radiator in the plane at height above an infinite electric ground plane is modelled as two patches separated by . The total number of triangular elements is 676. In the TCM analyser, a proper out-of-phase mode is selected (Figure 3).

Figure 3: Model of MPA above infinite ground plane for  mm, dominant mode TM₀₁ shown.

The resonant frequency of the dominant mode is shown as a function of height , see Figure 4. It has been evaluated from a modal resonant condition for eigenvalues employing an adaptive frequency sweep for each height. The behaviour is quite peculiar, especially for greater heights. For low heights ( mm or ), the resonant frequency decreases “regularly,” and quasianalytical formulas (see, e.g., [1, 3]) based on the fringing field concept are valid below this range. For mm () there is absolute minimum of the TM₀₁ resonant frequency. Further on, the resonant frequency rises to reach its maximum for mm (). Around this specific height the patch also shows the minimum of the radiation . The above described process repeats periodically. It is yet unclear to the authors as what is the physical background to the resonant frequency discontinuity around .

Figure 4: R50 × 30 resonant frequency of the dominant TM₀₁ mode. The dashed red curve is a quasianalytical equation from [1].

The terms , and obtained from (5)–(9) and eigenvalues are plotted at Figure 5 for  mm as a function of frequency. There is excellent agreement between the difference in stored energies and the eigenvalues, both obtained in a completely different manner.

Figure 5: Reactive energies and their differences for an R50 × 30 patch at height of 25 mm.

There is also very good agreement between the exact and confirming the validity of the proposed algorithm via (14), see Figure 6. Note that in (10) does not require the currents to be calculated on the structure while is evaluated in a rigorous way from modal currents (11).

Figure 6: The radiation for dominant mode of an R50 × 30 patch as a function of height .

From Figure 6 it is seen that the radiation has a minimum for a specific height. It is deduced that the reason lies in the cancelling of the radiated power between the two out-of-phase currents. Similar behaviour has been observed in the case of two half-wave thin-wire dipoles with opposite sinusoidal currents, separated by , see [25] for details. Actually these two out-of-phase dipoles may serve as a very simple model for a patch antenna with a dominant mode. When the dipoles are reduced to elementary (Hertzian) ones, an approximate analytical solution is available and in [25] we showed that the is led by the function After deriving (18), the condition is worked-out and the first nontrivial root of (19) could be approximated as [25] For sinusoidal currents on dipoles the minimum (evaluated numerically) occurs for .

The minimum of the patch under study is obtained at , a value that is remarkably close to the simple dipole model.

4.1. Algorithm Convergence

Since no other methods for calculating modal are available, is taken as a reference, and the relative error percentage is defined as: where is calculated from the currents using (11). Four different heights were chosen,  mm (),  mm (),  mm (), and  mm (), and the relative error was evaluated as a function of total triangular elements (including the mirror), see Figure 7. All quality factors were evaluated at the resonant frequency of the dominant mode for the R50 × 30 patch. As discussed earlier, the centroid approximation became more inaccurate with low heights . However, even for the lowest analyzed value , the relative error is in the order of a few percent for reasonable mesh density (hundreds of elements). Further improvements to the integration routine are considered for the future.

Figure 7: Relative error of the factor as a function of triangular elements (mesh density).

4.2. Fractional Bandwidth of the R50 × 30 Patch Antenna

It is known that the fractional bandwidth (FBW) is related to the unloaded factor and the desired matching VSWR level. For VSWR < we have [26] Using a full-wave simulator CST-MWS [27], an R50 × 30 patch has been simulated and the FBWCST for VSWR < 2 was calculated as: where and are margins for VSWR < 2 and is the centre frequency. Only very low heights were studied since we used a simple probe feed which introduces an inductance component to the total input impedance. The comparison in Figure 8 shows good agreement of both fractional bandwidths.

Figure 8: Fractional bandwidth FBW (VSWR < 2) for a R50 × 30 patch.

5. Applications: Fractal Antennas

In this section, a bit more complex structures will be studied. The first one (the “Self Affine U” fractal, SAU), has been described in [28] and further analyzed in [29]. This kind of radiating motif is employed as a dual-band radiator with mutually orthogonal radiation patterns at both bands. Therefore we are analyzing the first two modes, where the currents are orthogonal. These are depicted in Figures 9 and 10 for first (SAU1) and second (SAU2) fractal iteration, respectively. The current of the first (lower) mode has two out-of-phase components (see Figure 11 for schematic current paths) while the second mode comprises inphase currents only. As we know from previous studies, opposite currents contribute to a rapid increase of the radiation , and it is expected that will have a much higher than .

Figure 9: The first two characteristic modes (currents and charges) for the SAU1 structure.

Figure 10: The first two characteristic modes (currents and charges) for the SAU2 structure.

Figure 11: The main current paths for the first two modes of the SAU1/2 structure.

Figure 11 presents a very simple concept showing the main current paths for the and modes discussed above including the mirroring effect of the infinite ground plane. It could be simply stated that more opposing current paths lead to significant increase in .

We show detailed behaviour only for SAU2 (the situation is similar for SAU1)—see Figure 12 that confirms high for the mode. Characteristic angles are calculated for  mm, the actual height for which the dual-band antenna was designed [29].

(a)

(b)

Figure 12: Characteristic angles (left) and radiation for the SAU2.

5.1. The FCL-2 Fractal Antenna

The second presented structure is the so-called fractal clover leaf (FCL) of the second iteration, [14]. The antenna is fed by an L-probe [30] that excites its dominant mode and is located at height  mm. Actually, the dominant mode is composed of two degenerated modes and (Figure 13). The second higher mode is shown at Figure 14 for completeness.

(a)

(b)

(c)

(d)

Figure 13: Degenerated dominant mode , and of the FCL2 antenna (currents and charges).

(a)

(b)

Figure 14: Second higher mode (currents and charges).

Figure 15 shows the main current paths of these modes, and we can again deduce that the dominant mode will exhibit lower compared to . This is confirmed by Figure 18— has more than 200x higher radiation .

Figure 15: Schematic depiction of the dominant current paths for the dominant () and the second higher modes together with their modal radiation patterns.

6. Resonant Properties of Studied Antennas

The properties of studied antennas are summarized in this section. At first we observed that microstrip antenna could support different kinds of modes regarding their factors (see Figure 16):(a)low modes with the current flowing in one direction and not changing its phase (dominant modes of simple shapes like rectangular, circular patch, and so forth.)(b)high modes with part of the currents flowing in the opposite direction. These modes exist even on simple “U” shaped patch (Figure 9 left) and on complex (fractal) geometries.

Figure 16: Radiation s for different antennas/modes.

Secondly, it has been observed that resonant frequency is quite a complicated function of height. Unfortunately we do not yet have any physical explanation as to why some modes present minimum values of .

Looking at Figure 17, it is clear (and interesting) that the resonant frequency behaves quite differently for low- and high- modes. The resonant frequency of low- modes is much more sensitive to the height, whereas high- modes exhibit almost constant when the height is varied. The proposed explanation is that the opposite currents (responsible for high ) keep reactive fields very close to the radiating structure so the effect of a fringing field coupled to the ground plane becomes almost negligible.

Figure 17: Resonant frequencies for different antennas/modes.

Figure 18: Characteristic angles for the FCL2 structure at  mm.

7. Conclusions

Modal resonant properties of selected microstrip patch antennas have been studied with the help of characteristic modes and the novel theory published by Vandenbosch. It has been found that the resonant frequency of a simple rectangular patch antenna is quite a complicated function of its height above the infinite ground. Moreover, the dependency of resonant frequency is also found to be a function of the radiation factor (which is now possible to calculate in a rigorous way). Due to the complexity of the problem, no physical explanation for the resonant frequency behaviour has yet been found.

It is observed that the radiation factor decreases for “standard” heights (<~), however there exists an absolute minimum value of that has already been predicted by simple modeling of two elementary out-of-phase dipoles. Using proper feeding techniques (like the L-probe) allows us to design wideband compact antennas.

The theory now puts current distribution and the radiation factor into objective context. Whenever the current mode exhibits opposite components, high may appear.

Future work is needed to connect the presented theory with parameter sweeps or even optimization, so we will be able to design novel wideband/multimode compact antennas.

Acknowledgments

This work was supported by the Grant Agency of the Czech Technical University in Prague, Grant no. SGS11/065/OHK3/1T/13 and by the Project COST 1102. The authors would like to thank professor Vandenbosch for fruitful discussions, Neil Bell for his comments, and two anonymous reviewers who suggested some improvements to the paper.

References

P. Bhartia, I. Bahl, R. Garg, and A. Ittipiboon, Microstrip Antenna Design Handbook, Artech House, Norwood, Mass, USA, 2000.

K. F. Lee and W. Chen, Advances in Microstrip and Printed Antennas, John Wiley & Sons, New York, NY, USA, 1997.

J. R. James, P. S. Hall, and C. Wood, Microstrip Antenna Theory and Design, Peter Peregrinus, London, UK, 1986.

R. F. Harrington and J. Mautz, “Theory of characteristic modes for conducting bodies,” IEEE Transactions on Antennas and Propagation, vol. 19, no. 5, pp. 622–628, 1971.

E. A. Daviu, Analysis and design of antennas for wireless communications using modal methods, Ph.D. dissertation, Universidad Politécnica de Valencia, Valencia, Spain, 2008.

G. A. E. Vandenbosch, “Reactive energies, impedance, and Q factor of radiating structures,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 4, Article ID 5398856, pp. 1112–1127, 2010.

G. A. E. Vandenbosch and V. Volski, “Lower bounds for radiation Q of very small antennas of arbitrary topology,” in Proceedings of the 4th European Conference on Antennas and Propagation (EuCAP '10), pp. 1–4, Barcelona, Spain, April 2010.

R. F. Harrington, Field Calculation by the Method of Moments, IEEE Press, New York, NY, USA, 1993.

M. Cabedo-Fabres, E. Antonino-Daviu, A. Valero-Nogueira, and M. F. Bataller, “The theory of characteristic modes revisited: a contribution to the design of antennas for modern applications,” IEEE Antennas and Propagation Magazine, vol. 49, no. 5, pp. 52–68, 2007.

P. Hazdra and P. Hamouz, “On the modal superposition lying under the MoM matrix equations,” Radioengineering, vol. 17, no. 3, pp. 42–46, 2008.

Mathworks, http://www.mathworks.com/

S. N. Makarov, Antenna and EM Modeling with Matlab, John Wiley & Sons, New York, NY, USA, 2002.

S. M. Rao, D. R. Wilton, and A. W. Glisson, “Electromagnetic scattering by surfaces of arbitrary shape,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 3, pp. 409–418, 1982.

M. Capek, P. Hazdra, P. Hamouz, and M. Mazanek, “Software tools for efficient generation, modelling and optimisation of fractal radiating structures,” IET Microwaves, Antennas and Propagation, vol. 5, no. 8, pp. 1002–1007, 2011.

Comsol Multiphysics, http://www.comsol.com/

J. L. Volakis, Ch. Chen, and K. Fujimoto, Small Antennas: Miniaturization Techniques & Applications, chapter 2, McGraw-Hill, New York, NY, USA, 1st edition, 2010.

J. S. McLean, “A re-examination of the fundamental limits on the radiation q of electrically small antennas,” IEEE Transactions on Antennas and Propagation, vol. 44, no. 5, pp. 672–676, 1996.

R. F. Harrington and J. R. Mautz, “Control of radar scattering by reactive loading,” IEEE Transactions on Antennas and Propagation, vol. 20, no. 4, pp. 446–454, 1972.

R. F. Harrington, Time-Harmonic Electromagnetic Fields, IEEE Press, New York, NY, USA, 2001.

D. R. Rhodes, “Observable stored energies of electromagnetic systems,” Journal of the Franklin Institute, vol. 302, no. 3, pp. 225–237, 1976.

D. R. Rhodes, “Reactance Theorem,” Proceedings of the Royal Society A, vol. 353, no. 1672, pp. 1–10, 1977.

J. Shaeffer, “MOM3D method of moments code theory manual,” Research report 189594, Lockheed Advanced Development Company, 1992.

M. Capek, P. Hazdra, and J. Eichler, “A method for the evaluation of radiation Q based on modal approach,” IEEE Transactions on Antennas and Propagation. In press

P. Arcioni, M. Bressan, and L. Perregrini, “On the evaluation of the double surface integrals arising in the application of the boundary integral method to 3-D problems,” IEEE Transactions on Microwave Theory and Techniques, vol. 45, no. 3, pp. 436–439, 1997.

P. Hazdra, M. Capek, and J. Eichler, “Radiation Q-factors of thin-wire dipole arrangements,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 556–560, 2011.

A. D. Yaghjian and S. R. Best, “Impedance, bandwidth, and Q of antennas,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 4, pp. 1298–1324, 2005.

CST Gmbh, http://www.cst.com/

S. N. Sinha and M. Jain, “A self-affine fractal multiband antenna,” IEEE Antennas and Wireless Propagation Letters, vol. 6, Article ID 891519, pp. 110–112, 2007.

J. Eichler, P. Hazdra, M. Čapek, T. Kořínek, and P. Hamouz, “Design of a dual-band orthogonally polarized l-probe-fed fractal patch antenna using modal methods,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 1389–1392, 2011.

C. L. Mak, K. M. Luk, K. F. Lee, and Y. L. Chow, “Experimental study of a microstrip patch antenna with an L-shaped probe,” IEEE Transactions on Antennas and Propagation, vol. 48, no. 5, pp. 777–783, 2000.