Radiation from a Cavity-Backed Circular Aperture Array Antenna Enclosed by an FSS Radome
Radiation from a Cavity-Backed Circular Aperture Array Antenna Enclosed by an FSS Radome
Kim, Jihyung;Lee, Sangsu;Shin, Hokeun;Jung, Kyung-Young;Choo, Hosung;Park, Yong Bae
2018-11-22 00:00:00
applied sciences Article Radiation from a Cavity-Backed Circular Aperture Array Antenna Enclosed by an FSS Radome 1 , † 2 , † 2 3 4 Jihyung Kim , Sangsu Lee , Hokeun Shin , Kyung-Young Jung , Hosung Choo 2 , and Yong Bae Park * Hanwha Systems, Yongin 17121, Korea; jihyung.kim@hanwha.com Department of Electrical and Computer Engineering, Ajou University, Suwon 16499, Korea; lss1507@ajou.ac.kr (S.L.); hokeun0305@ajou.ac.kr (H.S.) Department of Electronic Engineering, Hanyang University, Seoul 04763, Korea; kyjung3@hanyang.ac.kr School of Electronic and Electrical Engineering, Hongik University, Seoul 04066, Korea; hschoo@hongik.ac.kr * Correspondence: yong@ajou.ac.kr; Tel.: +82-31-219-2358 † Equally contributed first authors. Received: 29 October 2018; Accepted: 19 November 2018; Published: 22 November 2018 Abstract: Radiation from a cavity-backed circular aperture array antenna enclosed by a frequency selective surface (FSS) radome is studied using the hybrid analysis method, by combining the mode matching method, the ray tracing technique, and Huygens’s principle. The equivalent magnetic surface currents on the apertures are derived from the aperture electromagnetic fields, which are calculated based on the mode matching method. The rays are generated from the equivalent magnetic surface currents and used to analyze the FSS radome based on the ray tracing technique. After being obtained from both the mode matching method and the ray tracing technique, electromagnetic fields on an outermost radome are transformed into the equivalent electric and magnetic surface currents using Huygens’s principle. The radiated fields are computed from the equivalent surface currents and compared with the measured data. Keywords: frequency selective surface radome; cavity-backed circular aperture array antenna; mode matching method; ray tracing technique 1. Introduction A radome, a portmanteau of radar and dome, is a structural and weatherproof enclosure, and thus it is used to protect microwave antennas. The radome usually comprises dielectric layers that minimally affect electromagnetic signals to be transmitted or received by the antenna. Due to the effects of the radome on the electromagnetic signals, a thorough analysis of the radome is needed. To understand electromagnetic properties of the radomes, there have been extensive studies on various radome configurations [1–10]. The radiation from a circular aperture surrounded by a hemisphere radome was predicted based on the dyadic Green’s function technique and physical optics (PO) method [2]. The Von Karman radome, one of the renowned radome structures, has been analyzed using the method of moments (MoM) [3], the coupled surface integral equation [4], and the aperture integration-surface integration (AI-SI) [5]. Furthermore, the radiation characteristic of hemisphere, tangent-ogive, and cone-shaped radomes with multiple sources has been investigated based on the multilevel fast multipole algorithm (MLFMA) [6]. The iterative physical optics-boundary integral-finite element method (IPO-BI-FEM) was used to analyze the sandwich tangent-ogive radome [7]. These studies focused on the radome consisting of dielectrics; however, the dielectric radome is generally responsible for the increase in radar cross section (RCS) of aircraft due to its broadband transmission characteristic. Instead, a frequency selective surface (FSS) radome, which Appl. Sci. 2018, 8, 2346; doi:10.3390/app8122346 www.mdpi.com/journal/applsci Appl. Sci. 2018, 8, 2346 2 of 9 is a curved FSS within a multi-layer radome, is employed to reduce RCS over a wide frequency range because it exhibits bandpass or bandstop characteristics. Previous studies have investigated electromagnetic characteristics of single- and multi-layer FSS radomes based on the pole residue matching (PRM) [8] and the ray tracing technique [9,10]. The FSS are one kind of the metasurfaces showing a bandpass or a bandstop characteristic. These characteristics are not found in natural materials. Metasurfaces, metamaterials, graphene, and plasmonics are gaining increasing popularity among researchers due to their ability to adjust permittivity or permeability purposely and manipulate electromagnetic waves passing through the materials. Due to numerous potential applications, a variety of studies have been conducted for sensing technology, plasmonics, graphene, and photonics [11–16]. There have also been studies on analyzing the properties of metamaterials and metasurfaces [17,18]. However, the investigation of electromagnetic properties for the complex structures which consist of a radome with the metasurfaces such as FSS and an aperture array antenna has not been presented. This is because it is hard to use the full-wave methods to analyze them, for a given memory usage and time, when an object is electrically large in size and complicated. Therefore, hybrid techniques combining the full-wave methods (such as MoM, FEM, finite-difference time-domain (FDTD), and the mode matching method) and the asymptotic methods (such as the ray tracing technique and geometric optics (GO)) are required, but related studies on radiation properties from array antennas enclosed by an FSS radome seem to be lacking. Therefore, it is of great significance to investigate electromagnetic properties of the aperture array antennas with the FSS radome using a hybrid technique combining the mode matching method for modeling aperture fields and the ray tracing technique for analyzing the radome. In this paper, we analyze the radiation from a cavity-backed circular aperture array antenna enclosed by an FSS radome using the hybrid analysis method combining the mode matching method, the ray tracing technique, and Huygens’s principle. Towards this purpose, in particular, three different ways are carried out step by step. Firstly, electromagnetic characteristics of the cavity-backed circular aperture array antenna is predicted based on the mode matching method [19]. The obtained tangential electric fields are transformed into the equivalent magnetic current sources which in turn become rays to be used in the subsequent step. While using the rays, the ray tracing technique is employed to analyze the multi-layer FSS radome [10], and reflection and transmission properties of the FSS layer come from FEM simulations on HFSS. After obtained from the mode matching method and the ray tracing technique, electromagnetic fields on the outermost radome are transformed into the equivalent electric and magnetic surface currents using Huygens’s principle, and radiation fields in the far-field region are computed from the equivalent surface currents. In brief, we take an advantage of the hybrid analysis method combining the mode matching method, the ray tracing technique, and HFSS, a commercial EM simulator, based on FEM for an analysis of the aperture array antenna enclosed by the FSS radome. To verify our formulation, the FSS radome enclosing the cavity-backed circular aperture array antenna is fabricated and our computation results are compared with the measured data. 2. Field Analysis Figure 1 shows the problem geometry. Assume that the z-oriented electric point source is located in a circular cavity with multiple circular apertures in a conducting plane. e and d are dielectric n n constant and loss tangent of each layer (n = 1, 2, 3, 4, 5). Electromagnetic properties and size of each layer are shown in Table 1. Appl. Sci. 2018, 8, 2346 3 of 9 Appl. Sci. 2018, 8, x FOR PEER REVIEW 3 of 9 Figure 1. Problem geometry. Figure 1. Problem geometry. Table 1. Design parameters of the frequency selective surface (FSS) radome. Table 1. Design parameters of the frequency selective surface (FSS) radome. Parameters Value Parameters Value Parameters Value Parameters Value D 284.8 mm L 290.5 mm 1 D1 284.8 mm L11 290.5 mm D 286.4 mm L 291.5 mm 2 2 D2 286.4 mm L2 291.5 mm D 292.3 mm L 295.2 mm 3 3 D3 292.3 mm L3 295.2 mm D 292.4 mm L 295.3 mm 4 4 D4 292.4 mm L4 295.3 mm D 298.4 mm L 299.1 mm 5 5 D 300 mm L 300 mm D5 298.4 mm L5 299.1 mm 6 6 e 4.35 d 0.0032 r1 1 D6 300 mm L6 300 mm e 1 d 0.0038 r2 2 ϵr1 4.35 δ1 0.0032 e 4.4 d 0.02 r3 3 ϵr2 1 δ2 0.0038 e 1 d 0.0038 r4 4 e ϵr3 4. 4.354 δd 3 0.02 0.0032 r5 5 ϵr4 1 δ4 0.0038 ϵr5 4.35 δ5 0.0032 Figure 2 shows an entire analysis procedure for the radiation from a cavity-backed circular aperture array antenna with the FSS radome based on the hybrid analysis method combining the Figure 2 shows an entire analysis procedure for the radiation from a cavity-backed circular mode matching method, the ray tracing technique, and HFSS, a commercial EM simulator, based aperture array antenna with the FSS radome based on the hybrid analysis method combining the on FEM. Firstly, for analysis of the antenna, we solved the electromagnetic boundary-value problem mode matching method, the ray tracing technique, and HFSS, a commercial EM simulator, based on of the circular cavity with the circular aperture array antenna based on the mode matching method FEM. Firstly, for analysis of the antenna, we solved the electromagnetic boundary-value problem of in our previous study [19]. The mode matching method provides an advantage in the fact that the circular cavity with the circular aperture array antenna based on the mode matching method in its solution is rigorous and theoretically robust. Moreover, it is time-efficient in an analysis of our previous study [19]. The mode matching method provides an advantage in the fact that its open-boundary problems compared to other numerical techniques such as FEM, FDTD, and MoM, solution is rigorous and theoretically robust. Moreover, it is time-efficient in an analysis of to name a few. An analysis using the mode matching method takes the following steps. Above all, open-boundary problems compared to other numerical techniques such as FEM, FDTD, and MoM, the whole region to be solved is divided into sub-regions on a basis of boundary. Electromagnetic to name a few. An analysis using the mode matching method takes the following steps. Above all, fields are defined in each region. Afterwards, boundary conditions are enforced to obtain a set the whole region to be solved is divided into sub-regions on a basis of boundary. Electromagnetic of simultaneous equations for modal coefficients. Matrix calculation enables an evaluation of the fields are defined in each region. Afterwards, boundary conditions are enforced to obtain a set of modal coefficients and electromagnetic fields in all regions can be calculated based on the obtained simultaneous equations for modal coefficients. Matrix calculation enables an evaluation of the modal coefficients. Note that the mode matching method should be used in the separable coordinate modal coefficients and electromagnetic fields in all regions can be calculated based on the obtained systems where eigen-modes can be defined in all regions. By using the tangential electromagnetic modal coefficients. Note that the mode matching method should be used in the separable coordinate field established from the mode matching method and the surface equivalence theorem [20], we can systems where eigen-modes can be defined in all regions. By using the tangential electromagnetic derive the equivalent magnetic surface currents on the apertures from the aperture fields, which will field established from the mode matching method and the surface equivalence theorem [20], we can be used in a formation of rays. Meanwhile, the FSS layer is a cross-loop dipole FSS which has a derive the equivalent magnetic surface currents on the apertures from the aperture fields, which will passband resonant frequency at 10 GHz. A detailed configuration and size of the designed FSS can be used in a formation of rays. Meanwhile, the FSS layer is a cross-loop dipole FSS which has a be found in the reference [10]. Reflection and transmission coefficients of the designed FSS layer are passband resonant frequency at 10 GHz. A detailed configuration and size of the designed FSS can obtained via numerous simulations using HFSS, an electromagnetic full-wave simulator. Reflection and be found in the reference [10]. Reflection and transmission coefficients of the designed FSS layer are transmission coefficients between dielectric layers, not involved with FSS layer, can be determined on obtained via numerous simulations using HFSS, an electromagnetic full-wave simulator. Reflection a polarization basis, as described in the reference [10]. Afterwards, to apply the ray tracing technique and transmission coefficients between dielectric layers, not involved with FSS layer, can be for the analysis of the radiation from the cavity-backed circular aperture antenna enclosed by the FSS determined on a polarization basis, as described in the reference [10]. Afterwards, to apply the ray tracing technique for the analysis of the radiation from the cavity-backed circular aperture antenna Appl. Sci. 2018, 8, 2346 4 of 9 Appl. Sci. 2018, 8, x FOR PEER REVIEW 4 of 9 enclosed by the FSS radome, we determined intercept points of the rays and the surfaces of the FSS radome, we determined intercept points of the rays and the surfaces of the FSS radome using the radome using the iterative method [1]. Then, we can generate the rays from the equivalent magnetic iterative method [1]. Then, we can generate the rays from the equivalent magnetic surface currents surface currents on the apertures and trace the rays passing through the radome. When analyzing on the apertures and trace the rays passing through the radome. When analyzing the multi-layer FSS the multi-layer FSS radome, we used the ray tracing technique and the reflection and transmission radome, we used the ray tracing technique and the reflection and transmission coefficients at the FSS coefficients at the FSS layer from the simulation of HFSS based on FEM [21]. Therefore, our hybrid layer from the simulation of HFSS based on FEM [21]. Therefore, our hybrid method provides a more method provides a more time-efficient method for analyzing the aperture array antenna enclosed by time-efficient method for analyzing the aperture array antenna enclosed by the multi-layer FSS radome the multi-layer FSS radome than the method using the full-wave analysis only. On the other hand, it than the method using the full-wave analysis only. On the other hand, it is noted that locally flat is noted that locally flat condition should be satisfied for the accuracy of the ray tracing technique. condition should be satisfied for the accuracy of the ray tracing technique. Since our radome structure Since our radome structure has smooth surfaces, the ray tracing technique can be employed. Rough has smooth surfaces, the ray tracing technique can be employed. Rough surfaces, ripples, and defects surfaces, ripples, and defects on the surface may decrease the accuracy. Additionally, the ray tracing on the surface may decrease the accuracy. Additionally, the ray tracing technique is applicable to the technique is applicable to the region where the far-field condition is satisfied. In our analysis, each region where the far-field condition is satisfied. In our analysis, each aperture is divided into hundreds aperture is divided into hundreds of small cells to apply the ray tracing technique. Then, we of small cells to apply the ray tracing technique. Then, we calculated the radiation field from each calculated the radiation field from each small cell, respectively, and summed up all the computed small cell, respectively, and summed up all the computed fields based on the superposition principle. fields based on the superposition principle. In this case, the far-field condition of each aperture is In this case, the far-field condition of each aperture is given by 2D / = 0.0267 m. Our ray tracing given by 2D /λ = 0.0267 m. Our ray tracing technique satisfies the far-field condition because the technique satisfies the far-field condition because the distance from the aperture center to the radome distance from the aperture center to the radome inner surface is larger than 0.142 m (4.75 λ). Lastly, inner surface is larger than 0.142 m (4.75 ). Lastly, we calculated the electromagnetic fields and the we calculated the electromagnetic fields and the equivalent electric and magnetic surface currents equivalent electric and magnetic surface currents over the radome’s outer surface via using the results over the radome’s outer surface via using the results from the ray tracing technique. Then, the from the ray tracing technique. Then, the equivalent currents on the radome’s outer surface can be equivalent currents on the radome’s outer surface can be used to calculate the radiation fields in the used to calculate the radiation fields in the far-field region. far-field region. Figure 2. Analysis procedure. Figure 2. Analysis procedure. 3. Numerical Results and Measurement 3. Numerical Results and Measurement Before proceeding with the analysis of radiation pattern of the 3 3 circular aperture array Before proceeding with the analysis of radiation pattern of the 3 × 3 circular aperture array antenna enclosed by the multi-layer FSS radome as shown in Figure 3, it is important to check the antenna enclosed by the multi-layer FSS radome as shown in Figure 3, it is important to check the accuracy of our mode matching formulation and to analyze the antenna properties in detail. Firstly, accuracy of our mode matching formulation and to analyze the antenna properties in detail. Firstly, the number of modes used is m = 19 ( direction) and n = 10 (z direction). In order to ensure the number of modes used is m= 19 (ϕ direction) and n = 10 (z direction). In order to ensure that the convergence was made, we tabulated the modal coefficients in the circular cavity in Table 2. that the convergence was made, we tabulated the modal coefficients in the circular cavity in Table 2. From Table 2, we figured out that the contributing modes are TM , TM , TM , TM , TM , TM , From Table 2, we figured out that the contributing mod 02 03es are 04 TM 05 06, ( 4)5 -4 5 TM , TM , TM , TM , TM , TM , TM , and TE . Also, we plotted normalized magnetic fields 02 03 04 05 06 07 45 02 Appl. Sci. 2018, 8, x FOR PEER REVIEW 5 of 9 Appl. Sci. 2018, 8, 2346 5 of 9 Appl. Sci. 2018, 8, x FOR PEER REVIEW 5 of 9 Appl. Sci. 2018, 8, x FOR PEER REVIEW 5 of 9 in the circular cavity at 10 GHz in Figure 4. It is seen that z-oriented electric point source affects the in the circular cavity at 10 GHz in Figure 4. It is seen that z-oriented electric point source affects the TM , TM , and TE . Also, we plotted normalized magnetic fields in the circular cavity at 10 GHz in TM and TE modes simultaneously. Figure 5 depicts the magnitude of electric fields on each aperture. 07 45 02 in the circular cavity at 10 GHz in Figure 4. It is seen that z-oriented electric point source affects the TM and TE modes simultaneously. Figure 5 depicts the magnitude of electric fields on each aperture. Figure 4. It is seen that z-oriented electric point source affects the TM and TE modes simultaneously. Note that the peak value occurs at the center of the middle aperture and the radiation becomes peak TM and TE modes simultaneously. Figure 5 depicts the magnitude of electric fields on each aperture. Note that the peak value occurs at the center of the middle aperture and the radiation becomes peak Figure 5 depicts the magnitude of electric fields on each aperture. Note that the peak value occurs at in a direction normal to the aperture array antenna. If the location of the electric point source in the Note that the peak value occurs at the center of the middle aperture and the radiation becomes peak in a direction normal to the aperture array antenna. If the location of the electric point source in the the center of the middle aperture and the radiation becomes peak in a direction normal to the aperture cavity is changed, the electric fields on the apertures and radiation patterns are also changed. It in a direction normal to the aperture array antenna. If the location of the electric point source in the cavity is changed, the electric fields on the apertures and radiation patterns are also changed. It array antenna. If the location of the electric point source in the cavity is changed, the electric fields on means that the equivalent magnetic currents on the apertures can affect the radiation properties of cavity is changed, the electric fields on the apertures and radiation patterns are also changed. It means that the equivalent magnetic currents on the apertures can affect the radiation properties of the apertures and radiation patterns are also changed. It means that the equivalent magnetic currents the aperture array antenna enclosed by the multi-layer FSS radome and the electric and magnetic means that the equivalent magnetic currents on the apertures can affect the radiation properties of the aperture array antenna enclosed by the multi-layer FSS radome and the electric and magnetic on currents on the apertur the radome surfac es can affect thee. radiation properties of the aperture array antenna enclosed by the the aperture array antenna enclosed by the multi-layer FSS radome and the electric and magnetic currents on the radome surface. multi-layer FSS radome and the electric and magnetic currents on the radome surface. currents on the radome surface. Figure 3. The 3 × 3 cavity-backed circular aperture array antenna. Figure 3. The 3 × 3 cavity-backed circular aperture array antenna. Figure 3. The 3 × 3 cavity-backed circular aperture array antenna. Figure 3. The 3 3 cavity-backed circular aperture array antenna. (a) (b) (a) (b) (a) (b) Figure 4. Normalized internal (within the cavity) magnetic field in dB scale at 10 GHz. (a) H Figure 4. Normalized internal (within the cavity) magnetic field in dB scale at 10 GHz. (a) H Figure 4. Normalized internal (within the cavity) magnetic field in dB scale at 10 GHz. (a) H (xz-plane), y y (xz-plane), (b) H (yz-plane). Figure 4. Normalized internal (within the cavity) magnetic field in dB scale at 10 GHz. (a) H (b (xz ) H -plane), ( (yz-plane). b) H (yz-plane). (xz-plane), (b) H (yz-plane). Figure 5. The magnitude of electric fields on each aperture. Figure 5. The magnitude of electric fields on each apertur e. Figure 5. The magnitude of electric fields on each aperture. Figure 5. The magnitude of electric fields on each aperture. Appl. Sci. 2018, 8, x FOR PEER REVIEW 6 of 9 Appl. Sci. 2018, 8, 2346 6 of 9 Table 2. Convergence behaviors of the modal coefficients of the cavity. | | A . m= 19 ⋯ m= 4 ⋯ m = 0 ⋯ m = 4 ⋯ m = 19 mn −14 −1 −1 −1 −14 n = 1 1.34 × T 1 able 0 2. Conver ⋯ −gence 3.89 ×behaviors 10 ⋯ of the 1.40 modal × 10 coef ⋯ ficients −3.59 of the × 10 cavity ⋯ . −1.36 × 10 −15 −2 −2 −15 n = 2 ⋯ ⋯ ⋯ ⋯ 2.25 × 10 −2.89 × 10 −2.05 −1.88 × 10 −2.29 × 10 jA j m = 19 m = 4 m = 0 m = 4 m = 19 mn −16 −2 −2 −16 n = 3 5.57 × 10 ⋯ 2.48 × 10 ⋯ 1.53 ⋯ 3.01 × 10 ⋯ −5.7 × 10 14 1 1 1 14 n = 1 1.34 10 3.89 10 1.40 10 3.59 10 1.36 10 −16 −1 −1 −1 −17 n = 4 ⋯ ⋯ ⋯ ⋯ 1.18 × 10 2.87 × 10 9.92 × 10 2.92 × 10 −1.75 × 10 15 2 2 15 n = 2 2.05 2.25 10 2.89 10 1.88 10 2.29 10 −17 −19 16 2 2 16 n = 5 1.64 × 10 ⋯ −2.41 ⋯ 1.82 ⋯ −2.44 ⋯ −1.59 × 10 n = 3 5.57 10 2.48 10 1.53 3.01 10 5.7 10 16 1 1 1 17 −20 −3 −3 −19 n = 4 1.18 10 2.87 10 9.92 10 2.92 10 1.75 10 n = 6 ⋯ ⋯ −1.82 ⋯ ⋯ −6.85 × 10 2.89 × 10 2.92 × 10 7.03 × 10 17 19 n = 5 1.64 10 2.41 1.82 2.44 1.59 10 −19 −6 −6 −19 n = 7 −7.44 × 10 ⋯ 6.69 × 10 ⋯ −3.42 ⋯ 7.02 × 10 ⋯ 1.99 × 10 20 3 3 19 n = 6 6.85 10 2.89 10 1.82 2.92 10 7.03 10 ⋮ ⋮ 19 ⋰ ⋮ 6 ⋰ ⋮ ⋱ ⋮ 6 ⋱ ⋮ 19 n = 7 7.44 10 6.69 10 3.42 7.02 10 1.99 10 −22 −8 −5 −8 −22 n = . 10 −7.33. × 10 ⋯ −9.2 × 10 . ⋯ −2 × 10 . ⋯ −9.3 × 10 . ⋯ 6.36 × . 10 . . . . . . . . . . . . . . . . . . |B | m= 19 ⋯ m= 4 ⋯ m = 0 ⋯ m = 4 ⋯ m = 19 22 8 5 8 22 mn n = 10 7.33 10 9.2 10 2 10 9.3 10 6.36 10 −15 −1 −2 −1 −15 n = 1 2.06 × 10 ⋯ −4.98 × 10 ⋯ −8.15 × 10 ⋯ −4.97 × 10 ⋯ −1.77 × 10 jB j m = 19 m = 4 m = 0 m = 4 m = 19 mn −16 −1 −1 −16 n = 2 3.56 × 10 ⋯ −1.7 × 10 ⋯ −1.52 ⋯ −1.66 × 10 ⋯ −2.91 × 10 15 1 2 1 15 n = 1 2.06 10 4.98 10 8.15 10 4.97 10 1.77 10 −17 −1 −1 −2 −17 n = 3 9.16 × 10 ⋯ −1 × 10 ⋯ 6.52 × 10 ⋯ −9.1 × 10 ⋯ −7.06 × 10 16 1 1 16 n = 2 3.56 10 1.7 10 1.52 1.66 10 2.91 10 17−17 1 −2 1−1 2−2 − 17 17 n = 3 n = 4 9.16 2.07 ×10 10 ⋯ −