Table of Contents
International Journal of Microwave Science and Technology
Volume 2008, Article ID 538216, 8 pages
Research Article

Characterization of Microstrip Ring with a Narrow Gap by an Iterative Method

1Laboratoire de Physique de la Matière Molle, Faculté des Sciences de Tunis, 2092 El Manar, Tunisia
2Laboratoire d'Electronique, Groupe Microonde ENSEEIHT, 31071 Toulouse, France

Received 6 December 2006; Revised 12 April 2007; Accepted 27 February 2008

Academic Editor: Safleddin Safavi-Naeini

Copyright © 2008 M. Yeddes 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.


This paper presents an implementation of an iterative method based on the wave's concept, for analysing a ring resonator. This method includes a two-dimensional fourier modal transformation (2D-FMT) in a wave guide environment. The method has the advantage of simplicity in that it does not involve basis functions and inversion of matrices, as used in other calculation methods. Therefore, it is capable of analyzing larger bodies than other classical techniques. An implementation of the iterative calculation is shown for the extraction of S parameters of microwave components. The good agreement between simulation results and experimental published data justifies the design procedure and validates the present analysis approach.

1. Introduction

In the Photonic Band Gap structures (PBG), an artificial lattice of dielectric can give rise to frequency regions where photons cannot propagate in any direction. Microwave circuits with a PBG structure are difficult to implement in compact form because such a structure requires at least five or six periods of units cells, and the size of one period is directly proportional to the wavelength of the center frequency of the PBG. In this paper, we show that a simple microstrip ring with a narrow gap can exhibit PBG characteristics.

An implementation of the iterative calculation is shown for extraction of ?????? parameters of a ring resonator. It consists in generating a recursive relationship between a given wave source and reflected waves at the interface containing the circuit. This discontinuity plane is divided into cells and characterized by a scattering operator matrix depending on boundary conditions. Then, a 2D-FMT algorithm is used to pass from spatial to spectral domain for each iteration of the recursive process [1]. The regular ring shows S parameters similar to those of a conventional transmission line; the introduction of a small gap in this regular ring induces two stop bands which can be tuned by a varactor incorporated into the narrow gap; this technique is detailed in the second part of this paper.

2. Formulation

The wave concept is introduced by writing the tangential electric field ???? and surface tangential current density ???? in terms of incident and reflected waves (see Figure 1) [1, 2]. It leads to the following set of equations: ????=12v???????????+?????????????,??=1,2,??=12v???????????-???????????,??=1,2,(1) where ????? is the surface tangential current density as ???????????, with ???? a unit vector normal to the interface O; ???? and ???? are the incident and reflected waves associated with the discontinuity interface O; ??0?? is the characteristic impedance of the middle ??, given by ??0??=v??0/??0??????, where ??????, is the relative permittivity of the region ??; ???? and ???? are the tangential electric field and magnetic field, respectively, on the plan O.

Figure 1: Definition of wave concept.

FMT: fourier modal transform.

On each region (region 1, 2), it is possible to definite a scattering matrix on spatial domain and reflection coefficient on spectral domain.

The interface O (Spatial domain) in which the circuit is defined is divided into pixels [3]; the interface contains three subdomains (dielectric, metal, and source).

By using the boundary condition in each subdomain (dielectric, metal, and source), it is possible to define the scattering matrix ????, who binds the incident waves to the reflected waves on the spatial domain [1, 3]: ????=????????????????????,(2) where ???=-????-((-1+??1+??2)/(1+??1+??2))????+((1-??2)/(1+??2))????, ????=(2??/(1+??2))????+(2??/(1+??1+??2))????, ?????=(2??/(1+??2))????+(2??/(1+??1+??2))????, and ??????=-????-((-1-??1+??2)/(1+??1+??2))????+((1-??2)/(1+??2))????, where ???=??01??02,??=??0v??01??02,??1=??0??01,??2=??0??02,(3)??????

The incident waves are reflected on the upper and lower parts of the metallic box. So the relationship between waves (?????? and (??=TE,TM)???TE????TM???=?GTE??00GTM??????TE????TM???,(5)) in spectral domain is [1, 3] G????=1-??0??????????,?????Coth????,??h???1+??0??????????,?????Coth????,??h???,??TE????,??=??????,????????0??0,??TM????,??=??????0??0??????,??,??????,??=?????????2+????????2-??20??????.(6) where h??

??: the ???? middle height.

This approach could then be considered as alternative between spatial and spectral domains. From the passage spatial to spectral domain, we use a 2D-FMT algorithm in wave guide environment (see Appendix B).

The iterative process continues to the nth iteration, when the electric field and the current density converge; the iterative process is stopped. It is possible to calculate ???? and ????=v???????????+?????,????=1v???????????-?????.(7) at the interface O: ??????

Finally, it is possible to determine the admittance matrix, from that, the scattering matrix ????=???????1-???1+??-1,(8) between the two ports can be obtained by the following equation: ???? [Y] is the admittance matrix.

The algorithm of the iterative process is given by Figure 2.

Figure 2: Iterative processes algorithm.

???? and O: the numbers of pixels in the plan ????,??.

????,?? and ????,??: incident wave and reflected wave, respectively, in the spectral domain.

????,?? and ??12: incident wave and reflected wave, respectively, in the spatial domain.

Figure 3 shows the ??12 convergence at 5?GHz for the structure given by Figure 4.

Figure 3: ??12 versus iterations number.
Figure 4: Ring resonator geometry.

We can see that the iterative process converges at 350 iterations. So we can calculate ??12 versus frequency, using at each frequency the same number of the iterations.

3. Results and Illustrations

Figure 5 shows the simulated and measured spectra of transmitted microwave (??12) of the microstip ring with the directly connected feedline.

Figure 5: Transmission coefficient of the ring resonator.

The geometry and dimensions of the regular ring are given by Figure 4; it is fabricated on a RT/Duriod 6010 substrate.

A comparison between the simulation results and measured results [4] is done; a good agreement is observed. The mismatch between the simulation and measurement data is inferior to 6?dB for the majority of frequencies. We note that the regular ring shows ??=1 parameter similar to those of conventional transmission line.

The introduction of a small gap in the closed loop of the regular ring induces strong attenuations around 6.5?GHz and 20?GHz corresponding to ??=3 and ????=????v2??????e?,(9) in (9): 2????=???? where ??e?=6.77 and ?? are given by simulation by [4], c is speed of light and ??=0.2 is wavelength.

The new geometry of the ring is given by Figure 6; the dimension of the gap is (??11)?mm. This structure induces multiple reflections with a fixed phase correlation.

Figure 6: Ring with gap.

The measured and simulated reflected coefficients (??21) are illustrated in Figure 7, and Figure 8 shows the transmitted coefficient ??=2. We note the absence of the stop band around 13?GHz correspond to 2????-??=???? in (9), because the microwave propagation along the ring is successively reflected with the fixed phase correlation when they satisfy ??11, where G is the gap dimension.

Figure 7: Reflection coefficient ??12 of the microstrip ring with a narrow gap.
Figure 8: Transmission coefficient ??11 of the ring with a narrow gap.

We see a satisfactory argument between, respectively, measured [4] and simulated reflected coefficient ??12 and transmitted coefficient (??11).

In the frequency range from 12?GHz to 17?GHz, the measurement data for 2???? show distinct standing wave behavior. This does not have any physical interest. They are disturbances because of the complexity of the structure.

Because a ring resonator can be represented by a parallel LC circuit of inductance L and capacitance C, introduction of a narrow gap in the ring slightly changes the values of C and L from the original ones. The value of C depends on the gap width G, the line width w, and the circumference 2????, while the value of L depends on the value of w and v??=1/????.

Thus the attenuations valleys around 6.5?GHz and 20?GHz, witch were explained to be due to multiple reflections occurring at the gap, correspond to the resonance frequencies of the equivalent LC circuit of the ring.

If we change the values of G, we can change the capacitance C, consequently the resonance frequency corresponding to ??12 takes another values.

Figure 9 shows the simulated transmitted coefficient (??1=??11??1+??12??2+??13??3,??(10a)2=??21??1+??22??2+??23??3,??(10b)3=??31??1+??32??2+??33??3.(10c)) of different values of G; we can observe the shift of the two attenuations valleys.

Figure 9: Transmission characteristics of the ring resonator for different value of G.

Another method can be used for varying the capacitance C; we can mount a varactor diode at the gap of the ring (see Figure 10), when a varactor is connected parallel to the gap, the total capacitance at the discontinuity is increased to the sum of the varactor capacitance and the gap capacitance.

Figure 10: Ring with varactor diode.

In simulation, the varactor diode is considered an auxiliary source that is used to study in a disjoined way two distinct parts of a circuit [5].

The expression of the density of current according to electric field for the two principal’s source and the auxiliary source will be according to following equations: (??1,??1)(??2,??2) and (??3,??3): density of current and electric field in the tow principal’s source.

??????: density of current and electric field in the auxiliary source.

????=(????????(0)????????(????+?????????)+???1-???/??+????????)????.(11) are the coupling admittance between different sources determinate with the iterative method.

Figure 11 shows the equivalent circuit of the varactor diode.

Figure 11: Equivalent circuit of the varactor diode.

According to this circuit, we can deduce the following relation: ????=??3?? With ????=??3?? and ??, ?? and ??=?? are the diode dimensions.

If ??3=(????????(0)?????????????+?????????+v1-??3/??+????????)??3.(12), we obtained ??3/?? The value of ??3=(????????(0)?????????????+?????????+(1-??3/2??)+????????)??3.(13) tends towards zero, therefore we will have ??31??1+??32??2=?????????-??33???3+??????????3?????????????+?????????+(1-??3./2??)(14) Using (13) and (10b), we have ??????????3?????????????+?????????+(1-??3˜/2??)??????????3?????????????+?????????.-1(15) Well, ??31??1+??32??2=?????????-??33+????????????????(????+???????????)-13???3=??31??1+??32??2????????-??33+????????/??????????????+??????????.-1(16) Consequently, (14) becomes ??3 Replacing ??1=???11+??13??1???1+???12+??13??2???2,??2=???21+??23??1???1+???22+??23??2???2.(17) in (10a) and (10b), we obtain ??1=??31????????-??33+????????/(????????(????+?????????,??)-12=??32????????-??33+????????/(????????(????+?????????.)-1(18) With ?????=11+??13??1??12+??13??2??21+??23??1??22+??23??2?.(19) The admittance matrix of the ring with varactor diode is given by (19) ??TE?(??,??)=??,????TE|||??(??,??)TE?,??(??,??)TM?(??,??)=??,????TM|||??(??,??)TM??.(??,??(A.1)

The simulated results for the transmission coefficient of ring with a narrow gap and the ring with the varactor mounted at the gap are shown in Figure 12; we can see the new valley obtained with the varactor.

Figure 12: Transmission coefficient with and without varactor diode.

4. Conclusion

In this paper, we have proposed a microstrip ring that exhibits PBG characteristics. The stop band of the ring resonator is caused by multiple reflections with a fixed phase correlation at the gap; we saw different technique for the control of those stopped bands; the results are determined by using an iterative technique. The comparison of numerical results with the measurement published data verified the validation of the proposed technique. This technique takes the advantage of simplicity and it is in conjunction with the 2D-FMT which allows a high computational speed and memory consumption.


A. Amplitude of the Mode in Metallic Box

The metallic box modes are developed as the following equations: ???TE??TM?=??????????????,????????-??=????????????????????.(A.2) The magnitude function is written as 1??=?????/??2+????/??2v2????????/????,????=?1ifmin=02ifmin?0.(A.3) With ??1 Using ??2 and ??TE=?(??,??)??,??????????|||??(??,??)???+?(??,??)??,??????????|||??(??,??)???,??(??,??)TM=?(??,??)??,??????????|||??(??,??)???+?(??,??)??,??????????|||??(??,??)???.(??,??)(A.4), it is possible to establish ?????=cos??????????sin?????????,?????=sin??????????cos?????????,????=????,????=????.(A.5) The base function for metallic box is given by (??3) Moreover, ??TE?(??,??)=??,???????????(??,??)cos??????????sin?????????+???,???????????(??,??)sin??????????cos?????????,??TM?(??,??)=??,???????????(??,??)sin??????????cos?????????+???,???????????(??,??)cos??????????sin?????????.(A.6) becomes ??(??,??)

B. Fast Modal Transform FMT

From the passage spatial to spectral domain in a rectangular wave guide environment, we use a 2D-FMT (Tow-Dimensions Fast Modal Transform), deduce from the two-dimensions Fourier transform.

The two-dimensions usual Fourier transform of the function ???(??,??)=2??,2?????,??????,????-2???????????/?????-2???????????/???.(B.1) can be written by the following equation: ??

According to ??(????,????), ??(????,????)=-??(2??-????,????) is an odd function, so ???(??,??)=??,2?????,??????,?????-2???????????/???-??-4????????2???????????/??????-2???????????/???,???(??,??)=??,2?????,??????,????-2??????×???2????????-2???????????/???-??-2????????2???????????/??????-2???????????/???,???(??,??)=??,2?????,??????,????-2??????×???-2???????????/??-????)-??2??????????????/???-????????-2???????????/???,???(??,??)=??,2?????,??????,??×??-2?????????-2???????????/???-??2??????????????/???????-2???????????/???,???(??,??)=-2????,2?????,??????,????-2???????sin2?????????????-2???????????/???.(B.2), ??(????,????) Of the other share ??(????,????)=??(????,2??-????) is an even function according to y, so, ???(??,??)=-2????,?????,??????,????-2???????sin???????????×???-2???????????/?????-4??????+??2???????????/????,???(??,??)=-2????,?????,??????,????-2???????sin?????????????-2??????×???-2????????-2???????????/???+??2????????2???????????/????,???(??,??)=-2????,?????,??????,????-2???????sin?????????????-2??????×???-2???????+????????/???+??2???????+????????/????,???(??,??)=-2????,?????,??????,????-2???????sin?????????????-2??????×???-2???????????/???+??2???????????/????,???(??,??)=-4????,?????,??????,????-2????????-2???????×sin2????????????cos2???????????.(B.3)????(??,??)=4??2????????2?????????(??,??).(B.4) Consequently, we can write the FMT using the usual Fourier transform:


  1. A. Mami, H. Zairi, A. Gharsallah, and H. Baudrand, “Analysis of microstrip spiral inductor by using iterative method,” Microwave and Optical Technology Letters, vol. 35, no. 4, pp. 302–306, 2002. View at Publisher · View at Google Scholar
  2. A. Gharsallah, A. Gharbi, L. Desclos, and H. Baudrand, “Analysis of interdigital capacitor and quasi-lumped miniaturized filters using iterative method,” International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, vol. 15, no. 2, pp. 169–179, 2002. View at Publisher · View at Google Scholar
  3. M. Yeddes, M. Kaddour, A. Gharsallah, and A. Gharbi, “Analysis of DGS structures by using the iterative method,” International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, vol. 18, no. 4, pp. 297–311, 2005. View at Publisher · View at Google Scholar
  4. C.-S. Kee, M.-Y. Jang, I. Park et al., “Photonic band gap formation by microstrip ring: a way to reduce the size of microwave photonic band gap structures,” Applied Physics Letters, vol. 80, no. 9, pp. 1520–1522, 2002. View at Publisher · View at Google Scholar
  5. M. Yeddes, H. Zairi, A. Gharsallah, A. Gharbi, and H. Baudrant, “A global electromagnetic modeling of an active microstrip structure,” International Journal of RF and Microwave Computer-Aided Engineering, vol. 17, no. 4, pp. 425–434, 2007. View at Publisher · View at Google Scholar