Implementation of Zero-Thickness Impedance Boundary Condition for Method of Moments and Application to Microstrip Antennas

Article information

J. Electromagn. Eng. Sci. 2022;22(3):386-388
Publication date (electronic) : 2022 April 20
doi : https://doi.org/10.26866/jees.2022.3.l.4
1Deparment of Electronic Engineering, Inha University, Incheon, Korea
*Corresponding Author: Il-Suek Koh (e-mail: ikoh@inha.ac.kr)
Received 2021 May 28; Revised 2021 July 7; Accepted 2021 July 13.

Abstract

To analyze a microstrip antenna consisting of a lossy metal, we proposed a scheme for a method of moments (MoM) to completely implement a zero-thickness impedance boundary condition (IBC). This model is more accurate for lossy metal than other boundary conditions, such as a resistive sheet model. The accuracy of the proposed scheme is numerically demonstrated by comparing antenna parameters, such as the scattering parameter, directivity, gain, and efficiency, computed by two commercial tools, FEKO and HFSS. These tools and the proposed scheme were applied to two microstrip antennas made from lossy metals with various conductivity.

I. Introduction

Various numerical methods have been used to analyze microstrip antennas, such as the method of moments (MoM) and the finite element method (FEM) [1, 2]. The accuracy and convergence properties of the methods can deteriorate in a thin imperfect conductor (lossy metal). Therefore, the metal is usually assumed to be a perfect electrical conductor (PEC) with zero thickness. It can be efficiently modeled by the electric field integral equation (EFIE) [1]. A resistive (impedance) sheet model is used for a lossy metal, but its accuracy may not be satisfactory [2].

In general, the impedance boundary condition (IBC) may be more accurate for lossy metal with a wide range of conductivity than the resistive sheet model. The combined field integral equation (CFIE) for the IBC is given by:

(1) [Y0EincZ0Hinc]tan=[jkL+ηZIYK-n^×I-ZK+ηZn^×IjkL+I]tan[JeJm]

The parameters in Eq. (1) are found in a previous study [3]. Due to the K-operator in Eq. (1), the IBC surface should be closed. Therefore, the zero-thickness lossy metal cannot be analyzed based on the IBC. We proposed a technique for a zero-thickness lossy metal based on the complete IBC formula in Eq. (1). We compared the proposed scheme with two commercial tools, FEKO and HFSS, in terms of the scattering parameter (|S11|), directivity (D), gain (G), and efficiency (eff.) of a microstrip antenna consisting of lossy metal.

II. Zero-Thickness Impedance Plate

The feeding source is modeled by the delta gap [1]. Einc and Hinc in Eq. (1) are given by Einc=lm Vin and Hinc ≈ 0, respectively. Here, lm and Vin are the edge length of the port and the source voltage. Here, it is assumed that Vin = 1V.

When the thickness (t) of the metal becomes zero, mesh 1 and mesh 2 on the opposite surfaces coincide with each other, as shown in Fig. 1. The opposite meshes facing the free-space and dielectric can be separately assembled. In addition, due to the non-penetrating properties of the IBC, no direct interactions happen among the meshes on the opposite surfaces. Therefore, strong direct interactions can be avoided. As a result, the zero-thickness IBC surface is closed, and the K-operator in Eq. (1) can be calculated exactly.

Fig. 1

Zero-thickness metal configuration.

We considered a PEC antenna (η = 0) made with very thin metal (17 or 0 μm thick). The antenna geometry is shown in Fig. 2(a); the units are in millimeter, and it operates at 76 GHz. The dielectric substrate is RO3003, which has a dielectric constant and loss tangent of 3.0 and 0.0013, respectively. The reference impedance is assumed to be 50Ω. The delta gap port has a separation of 40 μm from the substrate, as shown in Fig. 2(b) [4]. Table 1 shows a comparison of the results of the proposed IBC formula and FEKO, which have excellent agreement for the zero-thickness case. Fig. 2(c) shows the convergence rate of the IBC formula. The low thickness can generate strong interactions on the metal surface. Therefore, the convergence stagnates. But zero thickness can prevent strong interactions, as explained earlier, resulting in fast convergence. Hence, the results for the 0 μm case agree well with those of FEKO, but do not for the 17-μm case.

Fig. 2

Microstrip patch antenna geometry and convergence rate of iterative method of IBC MoM.

Simulation results of Fig. 2(a) antenna

For a larger structure, a 1-by-7 single-series-fed antenna is considered and operates at 10.1 GHz. It consists of PEC or copper plates with zero thickness. The detailed dimensions of the antenna can be found in another study [5]. The conductivity (σ) of the copper is 5.96 × 107 [S/m] [6]. The substrate of the antenna is Duroid5880, which has a dielectric constant and loss tangent of 2.2 and 0.0009, respectively. When the metal thickness is greater than the skin depth, the normalized impedance is calculated as η=πfε0/σ·(1+j). Here, f and ɛ0 are the frequency and permittivity of free space, respectively. For the copper, η / (1+ j) is given by 0.00007 at 10.1 GHz. Fig. 3(b) shows |S11| obtained from the proposed scheme and FEKO over a frequency range of 9.8–10.2 GHz. Table 3 shows a comparison of the antenna parameters computed by the proposed scheme and FEKO at 10.1 GHz. For the copper plate, the gain and efficiency have slightly larger discrepancies due to the different algorithms used to deal with the lossy metal.

Fig. 3

1-by-7 single-series-fed taper antenna array and |S11|.

Simulation results of Fig. 3(a) antenna

We investigated the effect of the conductivity on the final antenna simulation accuracy by considering the conductivities of several lossy metals, of which the type, and conductivity are given in [6]. The conductivity range includes that of the transparent conductor indium tin oxide (ITO), σ ≈1.3×106. The antenna in Fig. 2(a) is considered.

FEM can provide more accurate results for this kind of lossy conductor. Therefore, the results of HFSS were also compared. Fig. 4 shows |S11| and efficiency at 76 GHz, and for a clear comparison, the results of σ =1×108 and PEC are added. The PEC current can be used to estimate the absorbed power by the imperfect conductor, based on which the efficiency and |S11| can be modified [7]. This approximation, the proposed scheme, HFSS and FEKO results are compared in Fig. 4. In the FEKO simulation, the conductor is modeled by the impedance sheet. The accuracy of the proposed scheme is far over the whole conductivity range. For high conductivities, the results of FEKO and the proposed scheme are very similar, but slightly different from that of HFSS, because of the usage of different port models and the delta gap being located slightly away from the substrate in Fig. 2(b).

Fig. 4

|S11| and efficiency vs. conductivity.

III. Conclusion

A MoM scheme was proposed to deal with zero-thickness lossy metal based on the IBC and was applied to the analysis of a microstrip antenna. A lossy metal can usually be modeled by the IBC, so the advantage of the proposed scheme can provide more accurate results over a wide range of conductivity, which was numerically verified by comparing the results of the proposed scheme, FEKO, and HFSS. However, the major disadvantage of the proposed scheme is that the number of meshes increases, because meshes are required for the closed IBC metal, resulting in more memory required and computation time.

Acknowledgments

This work was supported by Institute of Information & communications Technology Planning & Evaluation (IITP) grant funded by the Korea government (MSIT) (No. 2019-0-00098, Advanced and Integrated Software Development for Electromagnetic Analysis).

References

1. Gibson WC. The Method of Moments in Electromagnetics 2nd Editionth ed. Chapman and Hall/CRC; 2014.
2. Jin JM, Volakis JL, Yu CL, Woo AC. Modeling of resistive sheet in finite element solutions. IEEE Transactions on Antennas and Propagation 40:727–731. Jun. 1992;
3. Lee H, Koh I, Yoo J. Modified SDIE MLFMM method for inhomogeneous impedance material containing PEC or PMC Region. The Journal of Korean Institute of Electromagnetic Engineering and Science 29(10):805–810. 2018;
4. Przesmycki R, Bugaj M, Nowosielski L. Broadband microstrip antenna for 5G wireless systems operating at 28GHz. Electronics 10(1):1. Jan. 2021;
5. Yuan T, Yuan N, Li LW, Leong MS. Design and analysis of phased antenna array with low sidelobe by fast algorithm. Progress in Electromagnetics Research 87:131–147. 2008;
7. Kiminki S, Ylä-Oijala P, Korkeakoulu T, Nimi T. Computation of antenna parameters using electromagnetic surface integral equation. MS dissertation Helsinki University of Technology; 2009.

Article information Continued

Fig. 1

Zero-thickness metal configuration.

Fig. 2

Microstrip patch antenna geometry and convergence rate of iterative method of IBC MoM.

Fig. 3

1-by-7 single-series-fed taper antenna array and |S11|.

Fig. 4

|S11| and efficiency vs. conductivity.

Table 1

Simulation results of Fig. 2(a) antenna

Thickness Simulator D [dB] G [dB]
17μm FEKO 7.6544 6.9486
proposed 5.3447 3.7413
0μm FEKO 7.6135 7.3448
proposed 7.5947 7.3343

Table 3

Simulation results of Fig. 3(a) antenna

Metal Simulator D [dB] G [dB] eff. [%]
PEC FEKO 14.9783 10.4067 34.9
proposed 14.8872 9.6751 30.12
Copper FEKO 14.9784 10.3979 34.83
proposed 14.8805 8.7992 24.65