# A New Approach for Calculation of Metamaterial Printed Fractal Antenna Using Galerkin’s Method

Authors: Amjad S. Saedy; Al-Samawal. Saleh; Sadiq. Ali
Journal: International Journal of Antennas and Propagation
 (2022)

---

## Abstract

This work provides a new approach for computing the impedance of a proposed multiband printed fractal antenna for wireless applications. Galerkin’s method is applied to deduce the impedance relationship of the proposed structure and then compute the return loss verses frequency by converting the impedance matrix of the proposed antenna [Z] to the scattering matrix [S]. This model is developed in order to study the impedance of the proposed antenna after adding a metamaterial structure in the antenna substrate. The obtained model is able to determine the resonant frequencies and the return loss of the proposed antenna. The model is also able to define the changes in these values when the dimensions of the proposed structure change. The proposed antenna provides multiband wireless applications in the (1–10) GHz frequency band, and the return loss of the proposed fractal antenna has been improved using negative permittivity and negative permeability metamaterial structure.


---

## Body

## 1. Introduction

The technical development in wireless communications creates multiapplications devices.These devices work on various frequencies or multiband frequencies for applications like Bluetooth, WiFi, WLAN, and WiMAX [1]. These applications require a new type of printed antenna that works on wide band frequencies or multiband frequencies. Fractal, multiscale antenna and multidimension antennas are the common proposed types [1].Fractal geometry has been used to design multiband antennas that resonate at different frequencies according to the effective lengths of the fractal structure.These antennas are characterized by different effective dimensions or several slots with a specific arrangement, printed on a substrate consisting of a material with a specified dielectric constant [2, 3].Various fractal structures, such as Cantor Bar, Koch, Tree, Minkowski, and Hilbert, were used for designing multiband antennas [4–6].Dual-band and triple-band fractal patch antennas were proposed in [7]. The structure was based on the Sierpinski carpet fractal, while a defected ground structure and a reflector plane were utilized for enhancing the gain of the antenna [7].Abdaljabar et al. [8] proposed a microstrip fractal patch antenna with a COVID-19 shape in order to create a miniature antenna for dual-band wireless, satellite, and radar applications.The proposed antenna operated and resonated on two frequencies of 7.5 GHz and 17 GHz within C and Ku bands, respectively.Souza et al. [9] presented a dual-band small-size microstrip antenna. Koch fractal contour structure was used to reduce the antenna size by 70%.Some studies suggested fractal patch antenna shapes which have different slots like triangular, rectangular, square and hexagonal. For example, a rectangular Minkowski fractal slot antenna was proposed in [10]. The proposed fractal antenna was used to increase the effective area and bandwidth of the antenna. This antenna had a compact size where the number of operating frequencies and the gain were higher than the reference antennas.Wide slot triband antenna for wireless applications was proposed in [11]. The proposed structure consisted of a wide circular and rectangular slot. Those slots were placed in the design to excite multiple resonating modes, at different frequencies.Recently, great interest has emerged in the interaction of electromagnetic waves with manufactured structures that provide negative electrical permittivity and negative magnetic permeability at certain frequencies. These structures are called metamaterials [12].Metamaterial structure has the property that controls the direction of electromagnetic radiation to gather the energy in a small angular domain [12, 13]. Besides that, the metamaterials antennas have been used to improve the performance of the antenna because it has unique band gap features and periodic structures. The metamaterials antenna has a structure that saves and re-radiates energy making its size small and behaves as a larger antenna [13]. For those previous reasons, the metamaterial structures have been used to increase the radiation power of the printed antenna.A metamaterial, inspired tapered patch antenna operating in 2.5 GHz, 3.5 GHz, 4.6 GHz, and 5.8 GHz, was proposed in [14]. Circular split-ring resonator (CSRR) metamaterial structure was placed near the radiating edges of the antenna to enhance the return loss of the antenna.Researchers in [15] introduced an antenna which has a rectangular patch on the front side along with two CSRRs and a metamaterial cover on the back side. Their approach was intended to improve the performance of the antenna for 2.4 GHz WiFi applications.In [16], metamaterial fractal defected ground structure was presented to obtain multiband wireless applications (2.8 GHz, 4.1–4.45 GHz, and 5.6 GHz). Koch snowflake fractal was etched from the four edges of the resonating patch to improve return loss. A double negative metamaterial unit cell was used to increase the gain and radiation efficiency of a compact patch antenna at the resonance frequencies.Researchers in [17] proposed a design of an ultrawideband fractal antenna. Metamaterial was located in the ground plane in order to enhance the gain and directivity of the antenna.Paper [18] presents a dual-band planar antenna for GSM and WiMAX applications. The researcher suggested adding two SRR cells to improve the directivity and gain and a partial ground plane to improve bandwidth.In paper [19], the ground plane was engineered by the placement of an SRR array to induce additional resonance due to the occurrence of magnetic dipole moment.Square-shaped SRR was proposed in [20] to obtain a dual-band fractal printed antenna operation with wideband operation at the higher frequency band.Most of the previous studies used simulators only to design the proposed antennas, and these antennas operated at specific frequencies.The main objective of this research is to present a mathematical model in order to design a multiband metamaterial printed fractal antenna. The proposed mathematical model gives wide design powers so that the antenna resonant frequencies can be changed according to the required application.


## 2. Materials and Methods

### 2.1. Galerkin’s Method

Galerkin’s method aims to study the discontinuity in the waveguide and calculate the equivalent impedance. This is done by equating the discontinuity in the waveguide in an electrical circuit and then calculating the equivalent impedance of that circuit using Kershof laws [21, 22].The discontinuity in the waveguide occurs due to a sudden change in waveguide dimensions. However, it can be caused due to a change in the type of the material. This discontinuity affects the propagation of modes in waveguide [21].


### 2.2. Design of Printed Antenna Using Galerkin’s Method

We assume that the patch antenna is placed in a waveguide, as shown in Figure1(a), and thus a discontinuity will occur in this waveguide. The impedance of the equivalent circuit (Figure 1(b)) is calculated using the Galerkin method. After that, the antenna impedance versus frequency (Z = R + j X) is studied. The resonant frequency of the antenna is the frequency at which the true part of the impedance is the greatest and the imaginary part is the least [21].Figure 1
Discontinuity in the waveguide (a) and equivalent circuit (b). a, b: the dimensions of the first waveguide (the substrate of the antenna). c, d: the dimensions of the second waveguide (patch antenna).Y1, Y2: the equivalent permittivity expressing the first and second waveguide. J1: the current resulting from the propagating modes in the first waveguide. J2: the current resulting from the propagating modes in the second waveguide. E1: the voltage source resulting from the propagating modes in the first waveguide. E2: The voltage source resulting from the propagating modes in the second waveguide. E: the voltage source expressing the discontinuity in the waveguide.
(a)(b)The return lossS11 versus frequency is studied by converting the impedance matrix [Z] of the antenna to the scattering matrix [S], to determine resonant frequencies at which the antenna can radiate or receive energy [23, 24].Figure1 shows the discontinuity in the waveguide and the equivalent electrical circuit.


### 2.3. The Impedance Relation

From the equivalent electrical circuit, equations (1) and (2) can be derived as follows:(1)E=E1=E2,(2)J=−J1−J2+Y1+Y2∗E.Equations (1) and (2) can be written as follows:(3)E1E2J=001001−1−1Y1+Y2J1J2E,(4)E1=V1f1′⇒J1=I1f1′E2=V2f1⇒J2=I2f1,E=∑qXqfq.V1 and V2 are the voltages resulting from the propagating modes in the first and second waveguide. Xq is the discontinuity impedance. fq is the propagating modes in the discontinuity region.Therefore, equation (3) becomes as follows:(5)V1V20=00f1′/fp00f1/fp−<fp/f1′>−<fp/f1>fpY1+Y2fqI1I2∑qXq.f1′/fp is the propagating modes in the first waveguide. f1/fp is the propagating modes in the second waveguide. fpY1+Y2fq is the permittivity of the circuit and all propagating modes in this circuit. The matrix equation (5) can be rewritten into the following matrix equation:(6)V1V20=00A00B−At−BtDI1I2XAtI1+BtI2=DX⇒X=D−1AtBtI1I2⇒V1V2=ABX⇒V1V2=ABD−1AtBtI1I2⇒Z=ABD−1AtBt.(1)
[A] is a matrix related to the frequencies propagated in the first waveguide.(2)
[B] is a matrix related to the frequencies propagated in the second waveguide.(3)
[D] is a matrix related to the frequencies propagated in the first and second waveguide and the permittivity of the first and second waveguide, it is given by equations (7)–(10).(7)D=∑fp/fTE2y1TE+y2TE+∑fp/fTm2y1TM+y2TM,(8)y1TE=γ1jwμ0cotγ1l1,y2TE=γ2jwμ0cotγ2l2,(9)y1TM=jwε0εrγ1cotγ1l1,y2TM=jwε0γ2cotγ2l2,(10)γ22=m2π2a2+n2π2b2−ko2,γ12=m2π2a2+n2π2b2−k02εr,(11)w=2πf,k0=2πλ,μ0=4π10−7,ε0=136π109.λm is the wavelength. cm/sec is the speed of light in a vacuum. fHZ is the frequency of the antenna. l1m is the thickness of the substrate, l2≫l1. εr is the relative dielectric constant of the substrate. γ1andγ2 are the propagation constants.


### 2.4. Antenna Geometry and Design

The proposed fractal antenna is designed according to the following stages.


#### 2.4.1. First Stage: Design of a Reference Antenna Operating at a Frequency of 2.4 GHz

According to [24], the dimensions of reference patch antenna using transmission line theory are calculated and shown in Table 1.Table 1
Dimensions of the proposed fractal antenna in mm.
LWεrSub hSub ySubxL1W149.441.42.21.57123.483.638.0754.85The antenna is designed on substrate (εr = 2.2) with dimensions 123.4 × 83.6 × 1.57 mm3.Table1 shows the dimensions of the reference antenna in mm.The reference antenna is simulated using HFSS software, and the return loss versus frequency is depicted in Figure2.Figure 2
The return loss versus frequency of the reference antenna (f = 2.4 GHz and S11 = −22 dB).


#### 2.4.2. Second Stage: Design of Multiband Fractal Antenna

The proposed fractal antenna is characterized by different effective dimensions or several slots with a specific arrangement as shown in Figure3(b).Figure 3
The proposed fractal antenna.
(a)(b)Fractal antenna structure is proposed using the Cantor Bar towards the two axes (X, Y). The antenna patch is divided into four smaller self-similar segments as shown in Figure 3(a), with a scaling ratio of 1/3. Then, a metal strip in the form of a sign (+) is added within the four segments connecting the antenna feeding line to this strip [25].The dimensions of this metal strip areX1 = L/3.5 = 10.97 mm and Y1=W/3.5 = 9.2 mm.


## 2.1. Galerkin’s Method

Galerkin’s method aims to study the discontinuity in the waveguide and calculate the equivalent impedance. This is done by equating the discontinuity in the waveguide in an electrical circuit and then calculating the equivalent impedance of that circuit using Kershof laws [21, 22].The discontinuity in the waveguide occurs due to a sudden change in waveguide dimensions. However, it can be caused due to a change in the type of the material. This discontinuity affects the propagation of modes in waveguide [21].


## 2.2. Design of Printed Antenna Using Galerkin’s Method

We assume that the patch antenna is placed in a waveguide, as shown in Figure1(a), and thus a discontinuity will occur in this waveguide. The impedance of the equivalent circuit (Figure 1(b)) is calculated using the Galerkin method. After that, the antenna impedance versus frequency (Z = R + j X) is studied. The resonant frequency of the antenna is the frequency at which the true part of the impedance is the greatest and the imaginary part is the least [21].Figure 1
Discontinuity in the waveguide (a) and equivalent circuit (b). a, b: the dimensions of the first waveguide (the substrate of the antenna). c, d: the dimensions of the second waveguide (patch antenna).Y1, Y2: the equivalent permittivity expressing the first and second waveguide. J1: the current resulting from the propagating modes in the first waveguide. J2: the current resulting from the propagating modes in the second waveguide. E1: the voltage source resulting from the propagating modes in the first waveguide. E2: The voltage source resulting from the propagating modes in the second waveguide. E: the voltage source expressing the discontinuity in the waveguide.
(a)(b)The return lossS11 versus frequency is studied by converting the impedance matrix [Z] of the antenna to the scattering matrix [S], to determine resonant frequencies at which the antenna can radiate or receive energy [23, 24].Figure1 shows the discontinuity in the waveguide and the equivalent electrical circuit.


## 2.3. The Impedance Relation

From the equivalent electrical circuit, equations (1) and (2) can be derived as follows:(1)E=E1=E2,(2)J=−J1−J2+Y1+Y2∗E.Equations (1) and (2) can be written as follows:(3)E1E2J=001001−1−1Y1+Y2J1J2E,(4)E1=V1f1′⇒J1=I1f1′E2=V2f1⇒J2=I2f1,E=∑qXqfq.V1 and V2 are the voltages resulting from the propagating modes in the first and second waveguide. Xq is the discontinuity impedance. fq is the propagating modes in the discontinuity region.Therefore, equation (3) becomes as follows:(5)V1V20=00f1′/fp00f1/fp−<fp/f1′>−<fp/f1>fpY1+Y2fqI1I2∑qXq.f1′/fp is the propagating modes in the first waveguide. f1/fp is the propagating modes in the second waveguide. fpY1+Y2fq is the permittivity of the circuit and all propagating modes in this circuit. The matrix equation (5) can be rewritten into the following matrix equation:(6)V1V20=00A00B−At−BtDI1I2XAtI1+BtI2=DX⇒X=D−1AtBtI1I2⇒V1V2=ABX⇒V1V2=ABD−1AtBtI1I2⇒Z=ABD−1AtBt.(1)
[A] is a matrix related to the frequencies propagated in the first waveguide.(2)
[B] is a matrix related to the frequencies propagated in the second waveguide.(3)
[D] is a matrix related to the frequencies propagated in the first and second waveguide and the permittivity of the first and second waveguide, it is given by equations (7)–(10).(7)D=∑fp/fTE2y1TE+y2TE+∑fp/fTm2y1TM+y2TM,(8)y1TE=γ1jwμ0cotγ1l1,y2TE=γ2jwμ0cotγ2l2,(9)y1TM=jwε0εrγ1cotγ1l1,y2TM=jwε0γ2cotγ2l2,(10)γ22=m2π2a2+n2π2b2−ko2,γ12=m2π2a2+n2π2b2−k02εr,(11)w=2πf,k0=2πλ,μ0=4π10−7,ε0=136π109.λm is the wavelength. cm/sec is the speed of light in a vacuum. fHZ is the frequency of the antenna. l1m is the thickness of the substrate, l2≫l1. εr is the relative dielectric constant of the substrate. γ1andγ2 are the propagation constants.


## 2.4. Antenna Geometry and Design

The proposed fractal antenna is designed according to the following stages.


### 2.4.1. First Stage: Design of a Reference Antenna Operating at a Frequency of 2.4 GHz

According to [24], the dimensions of reference patch antenna using transmission line theory are calculated and shown in Table 1.Table 1
Dimensions of the proposed fractal antenna in mm.
LWεrSub hSub ySubxL1W149.441.42.21.57123.483.638.0754.85The antenna is designed on substrate (εr = 2.2) with dimensions 123.4 × 83.6 × 1.57 mm3.Table1 shows the dimensions of the reference antenna in mm.The reference antenna is simulated using HFSS software, and the return loss versus frequency is depicted in Figure2.Figure 2
The return loss versus frequency of the reference antenna (f = 2.4 GHz and S11 = −22 dB).


### 2.4.2. Second Stage: Design of Multiband Fractal Antenna

The proposed fractal antenna is characterized by different effective dimensions or several slots with a specific arrangement as shown in Figure3(b).Figure 3
The proposed fractal antenna.
(a)(b)Fractal antenna structure is proposed using the Cantor Bar towards the two axes (X, Y). The antenna patch is divided into four smaller self-similar segments as shown in Figure 3(a), with a scaling ratio of 1/3. Then, a metal strip in the form of a sign (+) is added within the four segments connecting the antenna feeding line to this strip [25].The dimensions of this metal strip areX1 = L/3.5 = 10.97 mm and Y1=W/3.5 = 9.2 mm.


## 2.4.1. First Stage: Design of a Reference Antenna Operating at a Frequency of 2.4 GHz

According to [24], the dimensions of reference patch antenna using transmission line theory are calculated and shown in Table 1.Table 1
Dimensions of the proposed fractal antenna in mm.
LWεrSub hSub ySubxL1W149.441.42.21.57123.483.638.0754.85The antenna is designed on substrate (εr = 2.2) with dimensions 123.4 × 83.6 × 1.57 mm3.Table1 shows the dimensions of the reference antenna in mm.The reference antenna is simulated using HFSS software, and the return loss versus frequency is depicted in Figure2.Figure 2
The return loss versus frequency of the reference antenna (f = 2.4 GHz and S11 = −22 dB).


## 2.4.2. Second Stage: Design of Multiband Fractal Antenna

The proposed fractal antenna is characterized by different effective dimensions or several slots with a specific arrangement as shown in Figure3(b).Figure 3
The proposed fractal antenna.
(a)(b)Fractal antenna structure is proposed using the Cantor Bar towards the two axes (X, Y). The antenna patch is divided into four smaller self-similar segments as shown in Figure 3(a), with a scaling ratio of 1/3. Then, a metal strip in the form of a sign (+) is added within the four segments connecting the antenna feeding line to this strip [25].The dimensions of this metal strip areX1 = L/3.5 = 10.97 mm and Y1=W/3.5 = 9.2 mm.


## 3. Results and Discussion

### 3.1. Design of Proposed Fractal Antenna Using Galerkin Method

Galerkin’s method is applied to the proposed fractal antenna. The value of the equivalent impedance of the antenna versus frequency is then calculated.The resonant frequency of the antenna is the frequency at which the true part of the impedance is the greatest and the imaginary part is the least.Figures4(a)–4(c) show the resonant frequencies of the proposed fractal antenna 2.68 GHz, 5.15 GHz, and 7.72 GHz.Figure 4
(a) The first resonant frequency of the proposed fractal antenna (f = 2.68 GHz) . (b) The second resonant frequency of the proposed fractal antenna (f = 5.15 GHz). (c) The third resonant frequency of the proposed fractal antenna (f = 7.72 GHz).
(a)(b)(c)


### 3.2. Return Loss

After computing the impedance matrix [Z] of the proposed antenna, equation (12) is used to compute the return loss S11 versus frequency [23, 24].(12)s11=z11z22+z02−z12z21z11−z01z22+z02−z12z21.Figures5(a)–5(c) show the value of the return loss of the proposed fractal antenna.Figure 5
(a) Return loss of the proposed antenna (f = 2.68 GHz and S11 = −13 dB). (b) Return loss of the proposed antenna (f = 5.15 GHz and S11 = −18 dB). (c) Return loss of the proposed antenna (f = 7.72 GHz and S11 = −24 dB).
(a)(b)(c)The return loss of the proposed fractal antenna is −13 dB, −18 dB, and −24 dB at 2.68 GHz, 5.15 GHz, and 7.72 GHz, respectively. As a result, the proposed fractal antenna resonates at three frequency bands, thus increasing the number of resonant frequencies in which the antenna can radiate or receive energy (i.e., improving the bandwidth of the antenna).


### 3.3. Metamaterial Structure Design

The metamaterial structure consists of Split-Ring Resonators (SRRs) with straight metallic wires. The SRR cell consists of two rings made of a non-magnetic metal that is folded concentrically with a slit in each of them, as shown in Figure6 [12].Figure 6
Metamaterial structure.Table2 shows the dimensions of the proposed metamaterial structure which is used in this work. This structure gives negative electrical permittivity and negative magnetic permeability in the (4.8–5.8) GHz frequency band. The thickness of the substrate (h2) is 0.7 mm.Table 2
The dimensions of the proposed metamaterial structure in mm.
agwda1r1h2εr50.60.50.55.20.40.72.2TheS-parameters (reflection coefficient S11 and transmission coefficient S21) of the metamaterial structure are analyzed using HFSS software; then, the effective parameters of the metamaterial structure (such as the electrical permittivity (ε), magnetic permeability (μ), impedance (z), and the refractive index (n)) are calculated using Matlab program [16], as shown in Figure 7. As it can be seen, the real parts of the electrical permittivity, magnetic permeability, and refractive index reach a negative value in the (4.8–5.8) GHz band.Figure 7
(a) The impedance, (b) electrical permittivity, (c) magnetic permeability, and (d) refractive index, of the proposed metamaterial structure.
(a)(b)(c)(d)


### 3.4. Metamaterial Fractal Antenna Design

Figure8 shows the proposed metamaterial fractal antenna.Figure 8
The proposed metamaterial fractal antenna.The proposed metamaterial structure is placed in the substrate of the proposed antenna under the patch, at a height ofh2 = 0.7 mm above the ground plane, as shown in Figure 9.Figure 9
The substrate of the proposed metamaterial fractal antenna. (a) Front view. (b) Top view.
(a)(b)


### 3.5. Effective Permittivity and Effective Permeability of the Metamaterial Substrate

The addition of metamaterial structure in the substrate of the proposed fractal antenna will produce a negative electrical permittivity (ε2 < 0) and negative magnetic permeability (μ2 <0) in the (4.8–5.8) GHz frequency band, as shown in Figure 7.Therefore, the effective electrical permittivity (εeff) and effective magnetic permeability (μeff) of the equivalent substrate must be calculated using following equations [26, 27], as shown in Figure (10):(13)1εeff=∑n=1N1εnhnht,1μeff=∑n=1N1μnhnht.Figure 10
The double substrate and its equivalent single layer model.After the development of Galerkin’s method for computing the impedance and return loss, we obtain the frequencies as shown in Figure11.Figure 11
Frequencies of the metamaterial fractal antenna.The return loss versus frequency of the proposed metamaterial fractal antenna is depicted in Figure12; the figure shows that the return loss is −35 dB at 5 GHz.Figure 12
Return loss of the fractal antenna with loaded metamaterial structure.Figure13 shows a comparison of the return loss of the proposed antenna with and without metamaterial structure. The return loss has been improved by about 17 dB, using the proposed metamaterial structure.Figure 13
Comparison of return loss of the proposed antenna with and without metamaterial.To prove the obtained results of the designed model, HFSS simulator is used to design the proposed antenna.Figure14(a) illustrates the return loss of the proposed antenna using HFSS simulator.Figure 14
(a) Return loss of the metamaterial fractal antenna using HFSS simulator. (b) VSWR of the proposed antenna. (c) 2D radiation pattern of the proposed antenna at 5 GHz. (d) 3D radiation pattern and gain of the proposed antenna at 5 GHz. (e) The peak gain of the proposed antenna.
(a)(b)(c)(d)(e)We notice that the obtained parameters of the proposed antenna are good compared to those of the other previous studies.According to simulations, there are resonant frequencies at 2.7 GHz, 5 GHz, and 7.75 GHz.Good agreement between simulation and mathematical model results is observed.Figure14(b) shows the VSWR of the proposed antenna in (4–6) GHz frequency band. The simulated VSWR is 1.1 at 5 GHz.Figure14(e) shows the peak gain of the proposed antenna in (4–6) GHz frequency band.The gain of the proposed antenna is 9 dB at 5 GHz.


## 3.1. Design of Proposed Fractal Antenna Using Galerkin Method

Galerkin’s method is applied to the proposed fractal antenna. The value of the equivalent impedance of the antenna versus frequency is then calculated.The resonant frequency of the antenna is the frequency at which the true part of the impedance is the greatest and the imaginary part is the least.Figures4(a)–4(c) show the resonant frequencies of the proposed fractal antenna 2.68 GHz, 5.15 GHz, and 7.72 GHz.Figure 4
(a) The first resonant frequency of the proposed fractal antenna (f = 2.68 GHz) . (b) The second resonant frequency of the proposed fractal antenna (f = 5.15 GHz). (c) The third resonant frequency of the proposed fractal antenna (f = 7.72 GHz).
(a)(b)(c)


## 3.2. Return Loss

After computing the impedance matrix [Z] of the proposed antenna, equation (12) is used to compute the return loss S11 versus frequency [23, 24].(12)s11=z11z22+z02−z12z21z11−z01z22+z02−z12z21.Figures5(a)–5(c) show the value of the return loss of the proposed fractal antenna.Figure 5
(a) Return loss of the proposed antenna (f = 2.68 GHz and S11 = −13 dB). (b) Return loss of the proposed antenna (f = 5.15 GHz and S11 = −18 dB). (c) Return loss of the proposed antenna (f = 7.72 GHz and S11 = −24 dB).
(a)(b)(c)The return loss of the proposed fractal antenna is −13 dB, −18 dB, and −24 dB at 2.68 GHz, 5.15 GHz, and 7.72 GHz, respectively. As a result, the proposed fractal antenna resonates at three frequency bands, thus increasing the number of resonant frequencies in which the antenna can radiate or receive energy (i.e., improving the bandwidth of the antenna).


## 3.3. Metamaterial Structure Design

The metamaterial structure consists of Split-Ring Resonators (SRRs) with straight metallic wires. The SRR cell consists of two rings made of a non-magnetic metal that is folded concentrically with a slit in each of them, as shown in Figure6 [12].Figure 6
Metamaterial structure.Table2 shows the dimensions of the proposed metamaterial structure which is used in this work. This structure gives negative electrical permittivity and negative magnetic permeability in the (4.8–5.8) GHz frequency band. The thickness of the substrate (h2) is 0.7 mm.Table 2
The dimensions of the proposed metamaterial structure in mm.
agwda1r1h2εr50.60.50.55.20.40.72.2TheS-parameters (reflection coefficient S11 and transmission coefficient S21) of the metamaterial structure are analyzed using HFSS software; then, the effective parameters of the metamaterial structure (such as the electrical permittivity (ε), magnetic permeability (μ), impedance (z), and the refractive index (n)) are calculated using Matlab program [16], as shown in Figure 7. As it can be seen, the real parts of the electrical permittivity, magnetic permeability, and refractive index reach a negative value in the (4.8–5.8) GHz band.Figure 7
(a) The impedance, (b) electrical permittivity, (c) magnetic permeability, and (d) refractive index, of the proposed metamaterial structure.
(a)(b)(c)(d)


## 3.4. Metamaterial Fractal Antenna Design

Figure8 shows the proposed metamaterial fractal antenna.Figure 8
The proposed metamaterial fractal antenna.The proposed metamaterial structure is placed in the substrate of the proposed antenna under the patch, at a height ofh2 = 0.7 mm above the ground plane, as shown in Figure 9.Figure 9
The substrate of the proposed metamaterial fractal antenna. (a) Front view. (b) Top view.
(a)(b)


## 3.5. Effective Permittivity and Effective Permeability of the Metamaterial Substrate

The addition of metamaterial structure in the substrate of the proposed fractal antenna will produce a negative electrical permittivity (ε2 < 0) and negative magnetic permeability (μ2 <0) in the (4.8–5.8) GHz frequency band, as shown in Figure 7.Therefore, the effective electrical permittivity (εeff) and effective magnetic permeability (μeff) of the equivalent substrate must be calculated using following equations [26, 27], as shown in Figure (10):(13)1εeff=∑n=1N1εnhnht,1μeff=∑n=1N1μnhnht.Figure 10
The double substrate and its equivalent single layer model.After the development of Galerkin’s method for computing the impedance and return loss, we obtain the frequencies as shown in Figure11.Figure 11
Frequencies of the metamaterial fractal antenna.The return loss versus frequency of the proposed metamaterial fractal antenna is depicted in Figure12; the figure shows that the return loss is −35 dB at 5 GHz.Figure 12
Return loss of the fractal antenna with loaded metamaterial structure.Figure13 shows a comparison of the return loss of the proposed antenna with and without metamaterial structure. The return loss has been improved by about 17 dB, using the proposed metamaterial structure.Figure 13
Comparison of return loss of the proposed antenna with and without metamaterial.To prove the obtained results of the designed model, HFSS simulator is used to design the proposed antenna.Figure14(a) illustrates the return loss of the proposed antenna using HFSS simulator.Figure 14
(a) Return loss of the metamaterial fractal antenna using HFSS simulator. (b) VSWR of the proposed antenna. (c) 2D radiation pattern of the proposed antenna at 5 GHz. (d) 3D radiation pattern and gain of the proposed antenna at 5 GHz. (e) The peak gain of the proposed antenna.
(a)(b)(c)(d)(e)We notice that the obtained parameters of the proposed antenna are good compared to those of the other previous studies.According to simulations, there are resonant frequencies at 2.7 GHz, 5 GHz, and 7.75 GHz.Good agreement between simulation and mathematical model results is observed.Figure14(b) shows the VSWR of the proposed antenna in (4–6) GHz frequency band. The simulated VSWR is 1.1 at 5 GHz.Figure14(e) shows the peak gain of the proposed antenna in (4–6) GHz frequency band.The gain of the proposed antenna is 9 dB at 5 GHz.


## 4. Comparing the Results of This Research with a Number of Recent reference Studies

Table3 shows a comparison of the results of this research with some reference studies that used fractal geometry and metamaterial structures in the design of printed antennas.Table 3
Comparing the results of this research with a number of recent reference studies.
Ref.YearResonant frequency (GHz)Gain (dB)Return loss (dB)Antenna structureDesign method820227.5172.2−35Fractal patch antenna with a COVID-19 shapeHFSS2820223.52.6−40CSRR in the patch and 3∗ 3 metamaterial structures in the groundHFSS2920217.2810−19Fractal patch antenna and four SRR structures around the patchFEKO3020212.45.32.67−35Sierpinski fractal antenna and ground plane with slotsCST120212.5, 5.16.7, 7.78.6, 9.511.35.6−13.8, −10.84−13.48, −23.24−12.88, −14.17−26.09Fractal antenna and metamaterial structure under the patchHFSS620212.46−26Minkowski fractal antenna and metamaterial structure under the patchHFSSThis work20222.7, 5,7.759−14, −35, −20Fractal patch antenna and metamaterial structure under the patchHFSS and mathematical modelWe notice that the obtained parameters of the proposed antenna are good compared to those of the other previous studies.


### 4.1. Fabrication and Measured Results

The proposed metamaterial fractal antenna has been fabricated and tested in order to validate experimentally the design approach. The fractal antenna is fabricated on substrate (εr = 4.4) with dimensions 123.4 × 83.6 × 1.57 mm3 as shown in Figure 15(a).Figure 15
Fabricated proposed antenna. (a) Fractal antenna. (b) Metamaterial structure.
(a)(b)The metamaterial structure is fabricated on another substrate with the sameεr and dimensions as shown in Figure 15(b).The resonant frequency and return loss of the proposed antenna are measured using spectrum analyzer 8593E as shown in Figure16(a).Figure 16
Measurement setup for the proposed antenna.
(a)(b)Figure16(b) shows that the resonant frequency of the proposed fractal antenna is 5.7 GHz, and the return loss of the proposed fractal antenna is −41.99 dBTo prove the obtained measured results, HFSS simulator and mathematical model are used to design the proposed antenna. Figure17 shows the return loss versus frequency of the proposed antenna using HFSS simulator, while Figure 18 illustrates the return loss versus frequency of the proposed antenna using the mathematical model.Figure 17
The return loss versus frequency of the proposed antenna using HFSS simulatorf=5.6GHz;S11=−44dB.Figure 18
The return loss versus frequency of the proposed antenna using mathematical modelf=5.6GHz;S11=−40dB.Table 4
Comparison between simulation, mathematical model, and measured results.
Resonant frequency (GHz)Return loss (dB)Design method5.6−44HFSS software5.6−40Mathematical model5.74−41.99Measured results


## 4.1. Fabrication and Measured Results

The proposed metamaterial fractal antenna has been fabricated and tested in order to validate experimentally the design approach. The fractal antenna is fabricated on substrate (εr = 4.4) with dimensions 123.4 × 83.6 × 1.57 mm3 as shown in Figure 15(a).Figure 15
Fabricated proposed antenna. (a) Fractal antenna. (b) Metamaterial structure.
(a)(b)The metamaterial structure is fabricated on another substrate with the sameεr and dimensions as shown in Figure 15(b).The resonant frequency and return loss of the proposed antenna are measured using spectrum analyzer 8593E as shown in Figure16(a).Figure 16
Measurement setup for the proposed antenna.
(a)(b)Figure16(b) shows that the resonant frequency of the proposed fractal antenna is 5.7 GHz, and the return loss of the proposed fractal antenna is −41.99 dBTo prove the obtained measured results, HFSS simulator and mathematical model are used to design the proposed antenna. Figure17 shows the return loss versus frequency of the proposed antenna using HFSS simulator, while Figure 18 illustrates the return loss versus frequency of the proposed antenna using the mathematical model.Figure 17
The return loss versus frequency of the proposed antenna using HFSS simulatorf=5.6GHz;S11=−44dB.Figure 18
The return loss versus frequency of the proposed antenna using mathematical modelf=5.6GHz;S11=−40dB.Table 4
Comparison between simulation, mathematical model, and measured results.
Resonant frequency (GHz)Return loss (dB)Design method5.6−44HFSS software5.6−40Mathematical model5.74−41.99Measured results


## 5. Conclusions

A mathematical approach was proposed in this work for calculating the return loss and resonant frequencies of a multiband printed fractal antenna with metamaterial substrate. This model was generated by applying Galerkin’s method to calculate the impedance of the proposed structure. The mathematical approach was developed to study the effect of metamaterial structure for improving the return loss of the antenna.The return loss of the proposed fractal antenna has been improved using metamaterial structure. As a result, the radiation energy of the fractal antenna with metamaterial was greater than the energy of the fractal antenna without metamaterial.In order to validate the proposed methodology, metamaterial fractal antenna was fabricated and tested. The simulation, mathematical model, and measured results were almost the same.In future research, we propose to develop a mathematical approach to study the performance of the proposed fractal antenna using fractal distribution of metamaterial structures [7, 28, 29].



---

*Source: 5412200-2022-10-13.xml*