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.
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 Figure 1(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].

(a)

(b)
The return loss S11 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].
Figure 1 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:
Equations (1) and (2) can be written as follows:
and are the voltages resulting from the propagating modes in the first and second waveguide. is the discontinuity impedance. is the propagating modes in the discontinuity region.
Therefore, equation (3) becomes as follows:
is the propagating modes in the first waveguide. is the propagating modes in the second waveguide. 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:(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).
is the wavelength. is the speed of light in a vacuum. is the frequency of the antenna. is the thickness of the substrate, . is the relative dielectric constant of the substrate. 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.
The antenna is designed on substrate ( = 2.2) with dimensions 123.4 × 83.6 × 1.57 mm3.
Table 1 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 Figure 2.

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 Figure 3(b).

(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 are X1 = L/3.5 = 10.97 mm and Y1=/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.
Figures 4(a)–4(c) show the resonant frequencies of the proposed fractal antenna 2.68 GHz, 5.15 GHz, and 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].
Figures 5(a)–5(c) show the value of the return loss of the proposed fractal antenna.

(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 Figure 6 [12].

Table 2 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.
The S-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 (), and the refractive index ()) 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.

(a)

(b)

(c)

(d)
3.4. Metamaterial Fractal Antenna Design
Figure 8 shows the proposed metamaterial fractal antenna.

The proposed metamaterial structure is placed in the substrate of the proposed antenna under the patch, at a height of h2 = 0.7 mm above the ground plane, as shown in Figure 9.

(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 () and effective magnetic permeability () of the equivalent substrate must be calculated using following equations [26, 27], as shown in Figure (10):

After the development of Galerkin’s method for computing the impedance and return loss, we obtain the frequencies as shown in Figure 11.

The return loss versus frequency of the proposed metamaterial fractal antenna is depicted in Figure 12; the figure shows that the return loss is −35 dB at 5 GHz.

Figure 13 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.

To prove the obtained results of the designed model, HFSS simulator is used to design the proposed antenna.
Figure 14(a) illustrates the return loss of the proposed antenna using HFSS simulator.

(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.
Figure 14(b) shows the VSWR of the proposed antenna in (4–6) GHz frequency band. The simulated VSWR is 1.1 at 5 GHz.
Figure 14(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
Table 3 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.
We 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).

(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 Figure 16(a).

(a)

(b)
Figure 16(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 dB
To prove the obtained measured results, HFSS simulator and mathematical model are used to design the proposed antenna. Figure 17 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.


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].
Data Availability
There are no datasets in this research.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The experiments on the fabricated fractal antenna were applied at the communication laboratory in Tishreen University by assistance of Dr. Mohammad Asaad.