Design and Application of a Metamaterial Superstrate on a Bio-Inspired Antenna for Partial Discharge Detection through Dielectric Windows

The adaptation of dielectric windows as metamaterial superstrate over a bio-inspired Printed Monopole Antenna (PMA) was evaluated in order to improve the detection sensitivity of Ultra High Frequency (UHF) sensors designed for Partial Discharge (PD) measurement. For this purpose, rectangular and circular Split Ring Resonators (SRR) structures were designed and evaluated aiming to achieve a metamaterial superstrate that improves the characteristics of the bio-inspired PMA as the gain, bandwidth, and radiation pattern. Measurements of the PMA with metamaterial superstrate were carried out in an anechoic chamber and compared to the simulations performed. The results show that the metamaterial superstrate insertion did not impact the original operating bandwidth, covering most of the characteristic frequency range of PD activity. Moreover, this insertion resulted in a mean gain enhancement of 0.7 dBi regarding the reference PMA, resulting in an antenna with better sensitivity for PD detection (mean gain of 3.61 dBi). The PMA-metamaterial set PD detection sensitivity was evaluated through laboratory tests with a point-to-plane PD generator setup and in field with measurements from a 230 kV current transformer. The developed PMA-metamaterial set was able to detect, successfully, the activity of PD for both tests, being classified as an optimized sensor for PD detection through dielectric windows.


Introduction
High voltage equipment insulating systems are subjected to stressful conditions, such as intense electric fields, chemical reactions, mechanical efforts, temperature variations, and several environmental phenomena. Due to the occurrence of these stresses, low magnitude electrical discharges that partially short-circuit the insulating material may occur. This phenomenon is defined as partial discharge (PD) [1].
The continuous action of the PDs can induce a significant degradation of the insulating material, which can eventually lead to a full dielectric breakdown and, consequently, equipment failure. Therefore, the continuous monitoring of PD activity in high voltage equipment is crucial to prevent the development of dielectric problems [2].
The most traditional method used for PD detection is the one established by the IEC 60270 standard [3], which allows the measurement of current pulses emitted during the occurrence of PDs. However, this method is considered highly invasive, since it requires an electrical connection with the monitored equipment through a coupling capacitor.
In order to overcome the practical limitations of the IEC 60270 standard method for online monitoring, researchers have been studying several alternative methodologies for PD detection [4][5][6][7][8]. Among these methods, the radiometric ones are highlighted due to their non-invasive nature, since there is no need of an electrical connection between the sensor and the monitored equipment. The radiometric or Ultra High Frequency (UHF) method consists of the detection of electromagnetic waves emitted by PD current pulses, which propagate through the insulation system of high voltage equipment in the UHF frequency range (300 MHz to 3 GHz). The UHF operation frequency range isolates the detection system from most electromagnetic interference in substations, such as power electronics switching and corona discharges, since the signals generated by these phenomena have components with significant energy for frequencies between 200 and 300 MHz [9]. Therefore, the UHF method is characterized as a non-invasive and interference-resistant technique, making it attractive for the continuous monitoring of high voltage equipment.
In order to improve the detection sensitivity of UHF sensors, researchers studied the spectral features of the PD radiation, concluding that, generally, most of the radiated energy is concentrated between 300 MHz and 1.5 GHz [9,10]. One of the pioneering papers regarding the application of UHF sensors in PD detection is presented in [11], which reports the detection and location of PDs in a 420 kV Gas Insulated Substation (GIS), attesting the efficiency of the UHF method. Over the years, other researches have proved the efficiency and reliability of the UHF method in GISs [12,13], and its application was extended to other high voltage equipment, mainly power transformers and high voltage cable connections [10,14,15].
For power transformers and GIS, the allocation of UHF sensors is usually made through dielectric windows, as presented in Figure 1. Generally, a dielectric window can be defined as an opening covered by dielectric material on the equipment tank, providing a propagation way between the PD radiated waves and the installed UHF sensor [16]. The most traditional method used for PD detection is the one established by the IEC 60270 standard [3], which allows the measurement of current pulses emitted during the occurrence of PDs. However, this method is considered highly invasive, since it requires an electrical connection with the monitored equipment through a coupling capacitor.
In order to overcome the practical limitations of the IEC 60270 standard method for online monitoring, researchers have been studying several alternative methodologies for PD detection [4][5][6][7][8]. Among these methods, the radiometric ones are highlighted due to their non-invasive nature, since there is no need of an electrical connection between the sensor and the monitored equipment. The radiometric or Ultra High Frequency (UHF) method consists of the detection of electromagnetic waves emitted by PD current pulses, which propagate through the insulation system of high voltage equipment in the UHF frequency range (300 MHz to 3 GHz). The UHF operation frequency range isolates the detection system from most electromagnetic interference in substations, such as power electronics switching and corona discharges, since the signals generated by these phenomena have components with significant energy for frequencies between 200 and 300 MHz [9]. Therefore, the UHF method is characterized as a non-invasive and interference-resistant technique, making it attractive for the continuous monitoring of high voltage equipment.
In order to improve the detection sensitivity of UHF sensors, researchers studied the spectral features of the PD radiation, concluding that, generally, most of the radiated energy is concentrated between 300 MHz and 1.5 GHz [9,10]. One of the pioneering papers regarding the application of UHF sensors in PD detection is presented in [11], which reports the detection and location of PDs in a 420 kV Gas Insulated Substation (GIS), attesting the efficiency of the UHF method. Over the years, other researches have proved the efficiency and reliability of the UHF method in GISs [12,13], and its application was extended to other high voltage equipment, mainly power transformers and high voltage cable connections [10,14,15].
For power transformers and GIS, the allocation of UHF sensors is usually made through dielectric windows, as presented in Figure 1. Generally, a dielectric window can be defined as an opening covered by dielectric material on the equipment tank, providing a propagation way between the PD radiated waves and the installed UHF sensor [16].  [14,17].
The UHF sensors applied in dielectric windows for PD detection can assume several structures, such as loop electrode, disk electrode, conical coupler, and antennas (monopole, dipole, microstrip, and others) [18,19]. Among these sensors, printed monopole antennas (PMA) are highlighted for having desirable characteristics for practical applications, such as low cost, attractive radiation patterns, wideband, and ease of installation/construction [20]. However, PMA with traditional radiating element (patch) shapes designed to operate in UHF range would assume relatively large dimensions, limiting their use for practical applications in dielectric windows.
The operating frequency of PMA is directly related to the perimeter of the patch. The greater the perimeter, the lower the operating frequency. Therefore, techniques used for PMA miniaturization mainly seek optimized geometries that maximize the patch perimeter/area ratio. These optimized geometries can be addressed in the shapes of living beings developed in order to provide greater efficiency regarding the ability to survive. Hence, aiming to achieve a better radiating efficiency, bio-inspired antennas use the shapes of plants or animals as basis for their design [21][22][23][24][25][26].   [14,17].
The UHF sensors applied in dielectric windows for PD detection can assume several structures, such as loop electrode, disk electrode, conical coupler, and antennas (monopole, dipole, microstrip, and others) [18,19]. Among these sensors, printed monopole antennas (PMA) are highlighted for having desirable characteristics for practical applications, such as low cost, attractive radiation patterns, wideband, and ease of installation/construction [20]. However, PMA with traditional radiating element (patch) shapes designed to operate in UHF range would assume relatively large dimensions, limiting their use for practical applications in dielectric windows.
PMA have a simple structure and great ease of construction and installation. In addition, features such as large bandwidth make this type of antenna suitable for PD application [26]. The general structure of a PMA is presented in Figure 2, in which W0 and L0 are the antenna length and width, respectively, L and W are the length and width of the patch, respectively, Lg and Wg are the length and width of the truncated ground plane, respectively, Wf is the width of the transmission line, g is the distance between the patch and the ground plane, and h is the substrate thickness. According to [21], the electric current density in a PMA is greater at the edges of the patch. Therefore, an increase in the patch perimeter would also increase the wavelength and, consequently, decrease the lower operating frequency. Hence, the lower frequency of the PMA operating band can be approximated by the patch perimeter, as presented in Equation (1) [21].
where p is the patch perimeter and ε re f is the relative permittivity of the dielectric, approximated as [21]:

Metamaterials: Concept and Design
According to [28], from Maxwell equations, it is possible to classify the materials in four different classes according to their values of magnetic permeability (µ) and electrical permittivity (ε), as presented in Figure 3. According to [21], the electric current density in a PMA is greater at the edges of the patch. Therefore, an increase in the patch perimeter would also increase the wavelength and, consequently, decrease the lower operating frequency. Hence, the lower frequency of the PMA operating band can be approximated by the patch perimeter, as presented in Equation (1) [21].
where p is the patch perimeter and εre f is the relative permittivity of the dielectric, approximated as [21]:

Metamaterials: Concept and Design
According to [28], from Maxwell equations, it is possible to classify the materials in four different classes according to their values of magnetic permeability (µ) and electrical permittivity (ε), as presented in Figure 3. The first quadrant of Figure 3 (ε > 0, µ > 0) represents the materials in which the interaction between the electric and magnetic vectors (E and H, respectively) results in the propagation of a vector wave K with direction defined by the right hand rule, being denominated as RHM (Right Handed Materials). The RHM represents most of the materials commonly found in nature and presents a straight wave propagation mode [40].
The second quadrant (ε < 0, µ > 0) describes the electric plasmas in which evanescent waves propagate with higher efficiency [41]. The materials that are represented by this quadrant are known as ENG (Epsilon Negative). Similarly, the fourth quadrant (ε > 0, µ < 0) also presents greater propagation of evanescent waves, but for magnetic media. Therefore, the materials in this quadrant are called MNG (Mu Negative). Finally, the third quadrant (ε < 0, µ < 0) represents the metamaterials, or Double Negative (DNG) materials. Due to its double negative characteristic, the refractive index, calculated from Maxwell's equations, also assumes negative values. The negative refractive index condition reflects in an interaction between E and H that results in the propagation The first quadrant of Figure 3 (ε > 0, µ > 0) represents the materials in which the interaction between the electric and magnetic vectors (E and H, respectively) results in the propagation of a vector wave K with direction defined by the right hand rule, being denominated as RHM (Right Handed Materials). The RHM represents most of the materials commonly found in nature and presents a straight wave propagation mode [40].
The second quadrant (ε < 0, µ > 0) describes the electric plasmas in which evanescent waves propagate with higher efficiency [41]. The materials that are represented by this quadrant are known as ENG (Epsilon Negative). Similarly, the fourth quadrant (ε > 0, µ < 0) also presents greater propagation of evanescent waves, but for magnetic media. Therefore, the materials in this quadrant are called MNG (Mu Negative). Finally, the third quadrant (ε < 0, µ < 0) represents the metamaterials, or Double Negative (DNG) materials. Due to its double negative characteristic, the refractive index, calculated from Maxwell's equations, also assumes negative values. The negative refractive index condition reflects in an interaction between E and H that results in the propagation of a K wave vector governed by the left hand rule. Thus, metamaterials can also be referred as LHM (Left Handed Materials) and have a reverse propagation mode regarding RHM. Besides the presentation of the characteristics described for ENG and MNG materials, the LHM have a negative refractive index that enables these materials to be able to concentrate incident electromagnetic waves in the medium in which they are applied. Hence, LHM can be defined as converging lenses that increase the focus of electromagnetic waves, both incident and evanescent waves, in a given dielectric medium.
The use of LHM as electromagnetic wave lenses has a useful practical application in microwave engineering, since a higher concentration of electromagnetic waves on the surface of an antenna, for example, can result in a considerable increase in its gain [42][43][44].
However, materials such as ENG, MNG, and LHM are not found in nature, requiring the manufacture of resonant structures that emulate such materials.
The behavior of an ENG-like material can be emulated from Capacitive Load Strips (CLS) that, for certain frequency ranges, may have higher electron concentrations and an increase in their effective mass, resulting in a negative electrical permittivity [29]. The basic structure of a CLS is presented in Figure 4. However, materials such as ENG, MNG, and LHM are not found in nature, requiring the manufacture of resonant structures that emulate such materials.
The behavior of an ENG-like material can be emulated from Capacitive Load Strips (CLS) that, for certain frequency ranges, may have higher electron concentrations and an increase in their effective mass, resulting in a negative electrical permittivity [29]. The basic structure of a CLS is presented in Figure 4. For MNG-type materials, the emulation of their characteristics can be obtained from the use of a Split Ring Resonator (SRR). For specific frequency ranges, the electromagnetic interactions between the characteristic inductance of the rings and the capacitance present between their spacing results in high energy densities that characterize the material as MNG [30].
Therefore, from the joint application of CLS and SRR, it is possible to obtain a metamaterial that can be used as superstrate of antennas, resulting in the gain enhancement of these structures [36,39,[45][46][47].
A more detailed physic-mathematical explanation regarding the functioning of SRR and CLS structures, as well as the interaction between them, can be found in [29][30][31]41,48,49].
The metamaterial operating frequency can be defined from the resonance frequency of the SRR geometry adopted. Although the SRR can be designed with several shapes [50][51][52], the circular and rectangular shapes represent the most adopted geometries in practice due to their lower complexity and ease of construction. Therefore, in this paper, only circular and rectangular SRRs were evaluated for the design of the proposed metamaterial.
The general structures of circular and rectangular SRR are presented in Figure 5, respectively. For MNG-type materials, the emulation of their characteristics can be obtained from the use of a Split Ring Resonator (SRR). For specific frequency ranges, the electromagnetic interactions between the characteristic inductance of the rings and the capacitance present between their spacing results in high energy densities that characterize the material as MNG [30].
Therefore, from the joint application of CLS and SRR, it is possible to obtain a metamaterial that can be used as superstrate of antennas, resulting in the gain enhancement of these structures [36,39,[45][46][47].
A more detailed physic-mathematical explanation regarding the functioning of SRR and CLS structures, as well as the interaction between them, can be found in [29][30][31]41,48,49].
The metamaterial operating frequency can be defined from the resonance frequency of the SRR geometry adopted. Although the SRR can be designed with several shapes [50][51][52], the circular and rectangular shapes represent the most adopted geometries in practice due to their lower complexity and ease of construction. Therefore, in this paper, only circular and rectangular SRRs were evaluated for the design of the proposed metamaterial.
The general structures of circular and rectangular SRR are presented in Figure 5, respectively.
structures, as well as the interaction between them, can be found in [29][30][31]41,48,49]. The metamaterial operating frequency can be defined from the resonance frequency of the SRR geometry adopted. Although the SRR can be designed with several shapes [50][51][52], the circular and rectangular shapes represent the most adopted geometries in practice due to their lower complexity and ease of construction. Therefore, in this paper, only circular and rectangular SRRs were evaluated for the design of the proposed metamaterial. For a circular SRR, the resonance frequency can be approximated by Equations (3) and (4) [53].
where c represents the speed of the light, r1 and r2 represent the external and internal mean radius length, respectively, and Scirc is the existing gap in each ring with external and internal resonance frequencies defined by fC1 and fC2, respectively. Similarly, for the rectangular SRR [53]: For a circular SRR, the resonance frequency can be approximated by Equations (3) and (4) [53].
where c represents the speed of the light, r 1 and r 2 represent the external and internal mean radius length, respectively, and S circ is the existing gap in each ring with external and internal resonance frequencies defined by f C1 and f C2 , respectively. Similarly, for the rectangular SRR [53]: where l 1 and l 2 represent the external and internal side length, respectively, and W m is the width of each resonance ring with external and internal resonance frequencies defined by f R1 and f R2 , respectively. In order to verify if the designed SRR-CLS set represents a metamaterial (ε < 0, µ < 0), the ε and µ values can be estimated from the refraction index (η) and impedance (z) values, as presented in Equations (7) and (8), respectively [54].
where z and η can be defined in terms of coefficients of reflection (S 11 ) and transmission (S 21 ) of the designed set, as presented in Equations (9) and (10), respectively [54].
Re ln e inkl − Im ln e inkl , (10) in which, the S 11 and S 21 values are represented by its respective magnitude and phase components, Re[] and Im[] represent the real and imaginary components, respectively, k represents the wave number, l the longest metamaterial unit cell length, and e inkl is defined as [54]:

Material and Methods
The applied methodology was divided in three different stages: computational procedures, laboratory experimental tests, and field application.
In the computational procedures, the software High Frequency Structure Simulator (HFSS) from ANSYS Electronics Desktop was used for all the simulation and design processes of the metamaterial (individually and positioned as superstrate). In all the simulations, the substrate is a low-cost fiberglass (FR-4) with ε r = 4.4, h = 1.6 mm and loss tangent (δ) equal to 0.02. Moreover, the simulations were performed for the frequency of 200 MHz up to 1600 MHz with a discrete sweep type with 201 points. The solution frequency was centered in 1 GHz (central frequency for the designed antenna) with a maximum number of passes and Delta S equals to 6 and 0.02, respectively.
The simulated PMA and metamaterial bandwidths were defined as the entire frequency range at which the reflection coefficient values were below the −10 dB threshold. In addition, the gain values were also extracted for the individual PMA and the set PMA-metamaterial in order to evaluate the gain enhancement provided by the metamaterial application as superstrate on the bio-inspired designed PMA. Lastly, the radiation patterns were obtained through simulations in order to evaluate if the metamaterial superstrate application resulted into distortions in the original bio-inspired PMA patterns.
In the experimental procedures, the designed PMA and the set PMA-metamaterial were built and subjected to comparative measurements of reflection coefficient and gain in an anechoic chamber in order to reduce the effect of external interferences and signal reflections during the tests. Then, the set PMA-metamaterial was used in a PD generator setup in order to estimate the PD detection sensitivity of the set regarding IEC 60270 standard method.
Finally, the PMA-metamaterial set was subjected to field application for the detection of partial discharges in a 230 kV Current Transformer (CT) from a Brazilian substation, Mussuré II (Chesf-Companhia Hidro Elétrica do São Francisco).
In the following subsections, the computational procedures, laboratory experimental tests, and field application are detailed.

Bio-Inspired PMA Design
From the PMA model presented in Figure 2 and Equations (1) and (2), the bio-inspired PMA based on the Jatropha mollissima (Pohl) Baill leaf was obtained in [26]. The final dimensions of the designed bio-inspired PMA in [26] and a photograph of the Jatropha mollissima (Pohl) Baill leaf are presented in Figure 6. Finally, the PMA-metamaterial set was subjected to field application for the detection of partial discharges in a 230 kV Current Transformer (CT) from a Brazilian substation, Mussuré II (Chesf-Companhia Hidro Elétrica do São Francisco).
In the following subsections, the computational procedures, laboratory experimental tests, and field application are detailed.

Bio-Inspired PMA Design
From the PMA model presented in Figure 2 and Equations (1) and (2), the bio-inspired PMA based on the Jatropha mollissima (

Metamaterial Unit Cell Simulations
Through Equations (3)-(6) and dimensional fine adjustments, the resonance frequencies of the designed metamaterials unit cells were adjusted to present the DNG feature for the bandwidth within the PD activity frequency range (300-1500 MHz).
The final dimensions for the designed rectangular and circular SRR-CLS metamaterial unit cells are presented in Figure 7, respectively.

Metamaterial Unit Cell Simulations
Through Equations (3)-(6) and dimensional fine adjustments, the resonance frequencies of the designed metamaterials unit cells were adjusted to present the DNG feature for the bandwidth within the PD activity frequency range (300-1500 MHz).
The final dimensions for the designed rectangular and circular SRR-CLS metamaterial unit cells are presented in Figure 7, respectively.

Metamaterial Unit Cell Simulations
Through Equations (3)-(6) and dimensional fine adjustments, the resonance frequencies of the designed metamaterials unit cells were adjusted to present the DNG feature for the bandwidth within the PD activity frequency range (300-1500 MHz).
The final dimensions for the designed rectangular and circular SRR-CLS metamaterial unit cells are presented in Figure 7, respectively. In order to characterize the models presented in Figure 7 as metamaterials, each model was simulated individually before the application as superstrate on the PMA. For the unit cells evaluation in electric and magnetic terms (ε and µ), special boundary conditions were established as presented in Figure 8, in which each metamaterial unit cell model is surrounded by perfect electric conductor walls (PEC) and perfect magnetic conductor walls (PMC) in parallel and orthogonal orientations regarding the Z axis, respectively. The unit cell model excitation is performed by two In order to characterize the models presented in Figure 7 as metamaterials, each model was simulated individually before the application as superstrate on the PMA. For the unit cells evaluation in electric and magnetic terms (ε and µ), special boundary conditions were established as presented in Figure 8, in which each metamaterial unit cell model is surrounded by perfect electric conductor walls (PEC) and perfect magnetic conductor walls (PMC) in parallel and orthogonal orientations regarding the Z axis, respectively. The unit cell model excitation is performed by two waveports. In this way, it is possible to acquire the S 11 and S 21 values used for the calculation of ε and µ, as presented in Equations (5)- (9). From the ε and µ values, it was possible to classify the unit cells as metamaterials (ε < 0 and µ < 0). waveports. In this way, it is possible to acquire the S11 and S21 values used for the calculation of ε and µ, as presented in Equations (5)- (9). From the ε and µ values, it was possible to classify the unit cells as metamaterials (ε < 0 and µ < 0).

Metamaterial Superstrate Simulations
After the characterization as metamaterials, the simulated metamaterial unit cells were positioned in a 3 × 3 array as superstrates over the bio-inspired PMA as presented in Figure 9.

Metamaterial Superstrate Simulations
After the characterization as metamaterials, the simulated metamaterial unit cells were positioned in a 3 × 3 array as superstrates over the bio-inspired PMA as presented in Figure 9.

Metamaterial Superstrate Simulations
After the characterization as metamaterials, the simulated metamaterial unit cells were positioned in a 3 × 3 array as superstrates over the bio-inspired PMA as presented in Figure 9. The superstrate height was adjusted so that the reflection coefficient values of the bio-inspired PMA were minimally impacted. The optimal superstrate height was defined as 135 and 125 mm for the rectangular and circular SRRs, respectively. In addition, the final array arrangement was also adjusted so that the gain values provided by the metamaterial superstrates were maximized. These adjustments consisted in some physical parameter variations, such as, the interconnection or not between the CLS of each unit cell, the adjustment in the horizontal and vertical spacing between the cells, the number of cells, and the position of the array regarding the PMA. Then, the horizontal and vertical spacing between the SRR unit cells were defined, respectively, as 10 and 15 mm, for the rectangular SRR superstrate, and as 8 and 14 mm, for the circular SRR superstrate. Moreover, the CLS of both unit cell designs were separated by a 5 mm gap.
In addition, in order to attest that the gain improvement was actually coming from the inserted metamaterial structures and not from the dielectric material adopted as superstrate (FR-4), the bio-inspired PMA was also simulated with only the dielectric material as superstrate at similar height. After the application of the simulation settings presented at the beginning of this section in Figure 9 models, the following meshing presented in Figure 10 was obtained. The superstrate height was adjusted so that the reflection coefficient values of the bio-inspired PMA were minimally impacted. The optimal superstrate height was defined as 135 and 125 mm for the rectangular and circular SRRs, respectively. In addition, the final array arrangement was also adjusted so that the gain values provided by the metamaterial superstrates were maximized. These adjustments consisted in some physical parameter variations, such as, the interconnection or not between the CLS of each unit cell, the adjustment in the horizontal and vertical spacing between the cells, the number of cells, and the position of the array regarding the PMA. Then, the horizontal and vertical spacing between the SRR unit cells were defined, respectively, as 10 and 15 mm, for the rectangular SRR superstrate, and as 8 and 14 mm, for the circular SRR superstrate. Moreover, the CLS of both unit cell designs were separated by a 5 mm gap.
In addition, in order to attest that the gain improvement was actually coming from the inserted metamaterial structures and not from the dielectric material adopted as superstrate (FR-4), the bio-inspired PMA was also simulated with only the dielectric material as superstrate at similar height. After the application of the simulation settings presented at the beginning of this section in Figure 9 models, the following meshing presented in Figure 10 was obtained.

Reflection Coefficient and Gain Measurements
The reflection coefficient values were obtained through an ENA E5062A network analyzer from Agilent Technologies. For the gain measurements, an experimental setup composed by a reference antenna (Hyperlog 30100X: 380 MHz-10 GHz), with a 4.5 dBi mean gain, and the antenna under test (AUT) positioned beyond the far field distance (R) as presented in Figure 11.

Reflection Coefficient and Gain Measurements
The reflection coefficient values were obtained through an ENA E5062A network analyzer from Agilent Technologies. For the gain measurements, an experimental setup composed by a reference antenna (Hyperlog 30100X: 380 MHz-10 GHz), with a 4.5 dBi mean gain, and the antenna under test (AUT) positioned beyond the far field distance (R) as presented in Figure 11.
The selected far field distance was defined regarding the maximum operation frequency of PD main activity range in oil, i.e., 1500 MHz [9,10]. Therefore, from the final bio-inspired PMA dimensions and the far field equations presented in [20], the value of R was defined as 1.0 m.

Reflection Coefficient and Gain Measurements
The reflection coefficient values were obtained through an ENA E5062A network analyzer from Agilent Technologies. For the gain measurements, an experimental setup composed by a reference antenna (Hyperlog 30100X: 380 MHz-10 GHz), with a 4.5 dBi mean gain, and the antenna under test (AUT) positioned beyond the far field distance (R) as presented in Figure 11. The selected far field distance was defined regarding the maximum operation frequency of PD main activity range in oil, i.e., 1500 MHz [9,10]. Therefore, from the final bio-inspired PMA dimensions and the far field equations presented in [20], the value of R was defined as 1.0 m.
From the measured values of transmission and reflection coefficients and reference antenna gain (GR = 4.5 dBi), the designed bio-inspired PMA gain (GD) was calculated according to the adapted Friis equation [20]: where, S21aut and S21ref represent the transmission coefficient of the antenna under test and reference antenna, respectively, and S11aut and S11ref represent the reflection coefficient of the antenna under test and reference antenna, respectively. For the application of Equation (12) in the gain estimation, it was necessary to assess the performance of two measurements using the schematic presented in Figure 8. The first measurement was performed using two identical reference antennas (Hyperlog 30100X), one positioned at the AUT distance and the other at the reference antenna position. In this way, it was possible to obtain the reference transmission coefficient (S21ref). The second measurement was performed with the PMA Anechoic Chamber Reference Antenna AUT Support Support R Figure 11. Schematic of the experimental setup used in the gain measurement test.
From the measured values of transmission and reflection coefficients and reference antenna gain (G R = 4.5 dBi), the designed bio-inspired PMA gain (G D ) was calculated according to the adapted Friis equation [20]: where, S 21aut and S 21ref represent the transmission coefficient of the antenna under test and reference antenna, respectively, and S 11aut and S 11ref represent the reflection coefficient of the antenna under test and reference antenna, respectively. For the application of Equation (12) in the gain estimation, it was necessary to assess the performance of two measurements using the schematic presented in Figure 8. The first measurement was performed using two identical reference antennas (Hyperlog 30100X), one positioned at the AUT distance and the other at the reference antenna position. In this way, it was possible to obtain the reference transmission coefficient (S 21ref ). The second measurement was performed with the PMA positioned at the AUT distance, allowing the obtaining of the S 21aut . A photograph of the schematic shown in Figure 11 is presented in Figure 12. positioned at the AUT distance, allowing the obtaining of the S21aut. A photograph of the schematic shown in Figure 11 is presented in Figure 12. The PMA sensitivity was tested by checking the criteria established by [27] for PD detection applications, which recommends that PMAs must present a mean gain higher than 2 dBi for their respective operating bandwidth. The PMA sensitivity was tested by checking the criteria established by [27] for PD detection applications, which recommends that PMAs must present a mean gain higher than 2 dBi for their respective operating bandwidth.

PD Detection Sensitivity Tests
The experimental procedure applied for the generation and detection of PD consisted of a gradual increase of voltage, by means of a regulating transformer, until PD activity was generated in the PD source and detected by the bio-inspired PMA and the IEC 60270 standard method.
The bio-inspired PMA was positioned at the selected far field distance (1.0 m) from the PD source. In order to represent typical internal partial discharge signals presented in high voltage equipment, an oil cell composed by a point-to-plane electrode configuration spaced by a dielectric (polyamide disk), was used as PD generator. The standard method is represented by a coupling capacitor (1000 pF) and a resonant circuit (LDM-5) as presented in Figure 13. The PMA sensitivity was tested by checking the criteria established by [27] for PD detection applications, which recommends that PMAs must present a mean gain higher than 2 dBi for their respective operating bandwidth.

PD Detection Sensitivity Tests
The experimental procedure applied for the generation and detection of PD consisted of a gradual increase of voltage, by means of a regulating transformer, until PD activity was generated in the PD source and detected by the bio-inspired PMA and the IEC 60270 standard method.
The bio-inspired PMA was positioned at the selected far field distance (1.0 m) from the PD source. In order to represent typical internal partial discharge signals presented in high voltage equipment, an oil cell composed by a point-to-plane electrode configuration spaced by a dielectric (polyamide disk), was used as PD generator. The standard method is represented by a coupling capacitor (1000 pF) and a resonant circuit (LDM-5) as presented in Figure 13. The PD pulses were acquired using a Keysight oscilloscope DSO90604A with 6 GHz bandwidth, 20 GSa/s sampling rate, 70 ps rise time, and four analog channels.
In order to assess the antenna PD detection sensitivity in terms of apparent charge, all the calibration procedures established by the IEC 60270 standard were executed before the voltage  The PD pulses were acquired using a Keysight oscilloscope DSO90604A with 6 GHz bandwidth, 20 Gsa/s sampling rate, 70 ps rise time, and four analog channels.
In order to assess the antenna PD detection sensitivity in terms of apparent charge, all the calibration procedures established by the IEC 60270 standard were executed before the voltage application tests. For this, the LDC-5 calibrator from Doble Lemke was connected to the PD generator terminals and the pulses were collected by the LDM-5 system. The PD generator calibration apparent charge results are presented in Table 1. Since the results presented in Table 1 show an approximately linear relation between the apparent charge and the measured voltage values, it is possible to correlate the apparent charge PD values, generated by the oil cell, with the voltage values measured by the IEC 60270 standard method. Then, the antenna PD detection sensitivity in terms of apparent charge can be evaluated. Calibration procedures for radiometric free-space, such as presented in [55], can also be used for this purpose.

PD Detection in a 230 kV Current Transformer
The equipment selected for the field application tests was a 230 kV Current Transformer (CT), presented in Figure 14.

500
584 Since the results presented in Table 1 show an approximately linear relation between the apparent charge and the measured voltage values, it is possible to correlate the apparent charge PD values, generated by the oil cell, with the voltage values measured by the IEC 60270 standard method. Then, the antenna PD detection sensitivity in terms of apparent charge can be evaluated. Calibration procedures for radiometric free-space, such as presented in [55], can also be used for this purpose.

PD Detection in a 230 kV Current Transformer
The equipment selected for the field application tests was a 230 kV Current Transformer (CT), presented in Figure 14. Previous routine tests performed by the maintenance personnel, such as power factor tests, indicated CT insulating system degradation, pointed out the possibility of high partial discharges activity in this equipment. Therefore, this CT was selected for the field application, since a high partial discharge activity was expected. Such as for the PD sensitivity tests, the PD pulses were acquired in time domain. However, a Tektronix oscilloscope TDS5104B with 1 GHz bandwidth, 5 GSa/s sampling, and four analog channels, was used.
According to [1,56], a PD signal can be identified by the pulse location regarding a reference voltage signal. Usually, the PD signals occur in the phase values corresponding to the semicycle transitions of a sinusoidal reference voltage. In addition, the PD pulses present a relative periodicity Previous routine tests performed by the maintenance personnel, such as power factor tests, indicated CT insulating system degradation, pointed out the possibility of high partial discharges activity in this equipment. Therefore, this CT was selected for the field application, since a high partial discharge activity was expected. Such as for the PD sensitivity tests, the PD pulses were acquired in time domain. However, a Tektronix oscilloscope TDS5104B with 1 GHz bandwidth, 5 Gsa/s sampling, and four analog channels, was used.
According to [1,56], a PD signal can be identified by the pulse location regarding a reference voltage signal. Usually, the PD signals occur in the phase values corresponding to the semicycle transitions of a sinusoidal reference voltage. In addition, the PD pulses present a relative periodicity with a phase displacement of 180 • between each pulse. Therefore, in order to evaluate if the signals emitted by the 230 kV CT and detected by the antenna were originated by PD, a signal generator was used in order to produce a sinusoidal reference voltage.
Moreover, the bio-inspired antenna without the metamaterial superstrate was used as reference in the field applications for comparative purposes. For this, both antennas were positioned simultaneously at the same distance from the CT. In addition, cables with the same length and attenuation were used. The signals were obtained simultaneously by triggering the oscilloscope.

Results
The analysis of the results is divided into three subsections. The first is about the simulation results regarding the bio-inspired PMA, metamaterial unit cell, and superstrate. The second subsection is regarding the practical results regarding reflection coefficient and gain measurements, as well the PD detection sensitivity. Finally, Section 3 is regarding the field application of the designed bio-inspired PMA-metamaterial set on the PD detection in a 230 kV CT.

Simulation Results
As previously mentioned, the bio-inspired PMA presented in [26] meets the size, bandwidth, and gain values requirements for PD detection through dielectric windows in high voltage equipment. However, although the PMA mean gain is higher than the 2 dBi criteria threshold, it can be verified in Figure 15 that the bio-inspired PMA presented maximum gain values lower to the recommended threshold for the frequency values near to the lower operating frequency (487 MHz). Therefore, in order to explore the dielectric windows for the PMA gain enhancement, metamaterial unit cells were simulated for the application as dielectric superstrates. simultaneously at the same distance from the CT. In addition, cables with the same length and attenuation were used. The signals were obtained simultaneously by triggering the oscilloscope.

Results
The analysis of the results is divided into three subsections. The first is about the simulation results regarding the bio-inspired PMA, metamaterial unit cell, and superstrate. The second subsection is regarding the practical results regarding reflection coefficient and gain measurements, as well the PD detection sensitivity. Finally, the third section is regarding the field application of the designed bio-inspired PMA-metamaterial set on the PD detection in a 230 kV CT.

Simulation Results
As previously mentioned, the bio-inspired PMA presented in [26] meets the size, bandwidth, and gain values requirements for PD detection through dielectric windows in high voltage equipment. However, although the PMA mean gain is higher than the 2 dBi criteria threshold, it can be verified in Figure 15 that the bio-inspired PMA presented maximum gain values lower to the recommended threshold for the frequency values near to the lower operating frequency (487 MHz). Therefore, in order to explore the dielectric windows for the PMA gain enhancement, metamaterial unit cells were simulated for the application as dielectric superstrates. The reflection coefficient results for the metamaterial unit cells with circular and rectangular SRR are presented, respectively, in Figure 16. The reflection coefficient results for the metamaterial unit cells with circular and rectangular SRR are presented, respectively, in Figure 16.
in the field applications for comparative purposes. For this, both antennas were positioned simultaneously at the same distance from the CT. In addition, cables with the same length and attenuation were used. The signals were obtained simultaneously by triggering the oscilloscope.

Results
The analysis of the results is divided into three subsections. The first is about the simulation results regarding the bio-inspired PMA, metamaterial unit cell, and superstrate. The second subsection is regarding the practical results regarding reflection coefficient and gain measurements, as well the PD detection sensitivity. Finally, the third section is regarding the field application of the designed bio-inspired PMA-metamaterial set on the PD detection in a 230 kV CT.

Simulation Results
As previously mentioned, the bio-inspired PMA presented in [26] meets the size, bandwidth, and gain values requirements for PD detection through dielectric windows in high voltage equipment. However, although the PMA mean gain is higher than the 2 dBi criteria threshold, it can be verified in Figure 15 that the bio-inspired PMA presented maximum gain values lower to the recommended threshold for the frequency values near to the lower operating frequency (487 MHz). Therefore, in order to explore the dielectric windows for the PMA gain enhancement, metamaterial unit cells were simulated for the application as dielectric superstrates. The reflection coefficient results for the metamaterial unit cells with circular and rectangular SRR are presented, respectively, in Figure 16. Both the designed metamaterial unit cells presented bandwidths within the PD activity frequency range and for the frequency values that presented lower gain results for the bio-inspired PMA. The circular and rectangular unit cells presented operating bands up to 1033 and 1117 MHz, respectively.
Besides the reflection coefficient, the transmission coefficient values were also obtained as presented in Figure 17. Therefore, through the magnitude and phase values ( Figure 18 and Figure 19) extracted from both coefficients (reflection and transmission), it was possible to calculate the permittivity (ε) and permeability (µ) values through Equations (5)- (9). The real values of µ and ε for both metamaterial unit cells are presented in Figure 20.
PMA. The circular and rectangular unit cells presented operating bands up to 1033 and 1117 MHz, respectively.
Besides the reflection coefficient, the transmission coefficient values were also obtained as presented in Figure 17. Therefore, through the magnitude and phase values ( Figure 18; Figure 19) extracted from both coefficients (reflection and transmission), it was possible to calculate the permittivity (ε) and permeability (µ) values through Equations (5)- (9). The real values of µ and ε for both metamaterial unit cells are presented in Figure 20.

respectively.
Besides the reflection coefficient, the transmission coefficient values were also obtained as presented in Figure 17. Therefore, through the magnitude and phase values ( Figure 18; Figure 19) extracted from both coefficients (reflection and transmission), it was possible to calculate the permittivity (ε) and permeability (µ) values through Equations (5)- (9). The real values of µ and ε for both metamaterial unit cells are presented in Figure 20.

respectively.
Besides the reflection coefficient, the transmission coefficient values were also obtained as presented in Figure 17. Therefore, through the magnitude and phase values ( Figure 18; Figure 19) extracted from both coefficients (reflection and transmission), it was possible to calculate the permittivity (ε) and permeability (µ) values through Equations (5)- (9). The real values of µ and ε for both metamaterial unit cells are presented in Figure 20.  From Figure 20, it can be verified that both metamaterial unit cells presented double negative characteristic of ε and µ for almost the entire range of PD activity, with the exception of the bands corresponding to 1190-1350 MHz for the rectangular SRR and 1100-1300 MHz for the circular SRR, where the values of µ assume positive magnitudes. Therefore, both the SRR unit cells were classified as metamaterials applicable as superstrate for the designed bio-inspired PMA.
After the classification as metamaterials, the SRR unit cells were positioned as superstrates, as From Figure 20, it can be verified that both metamaterial unit cells presented double negative characteristic of ε and µ for almost the entire range of PD activity, with the exception of the bands corresponding to 1190-1350 MHz for the rectangular SRR and 1100-1300 MHz for the circular SRR, where the values of µ assume positive magnitudes. Therefore, both the SRR unit cells were classified as metamaterials applicable as superstrate for the designed bio-inspired PMA.
After the classification as metamaterials, the SRR unit cells were positioned as superstrates, as presented in Figure 9 and described in Section 4. The first analysis regarding the set PMA-Metamaterial is concerning the influence on the reflection coefficient due the application of the superstrate over the bio-inspired PMA. The results for this analysis are presented in Figure 21. From Figure 20, it can be verified that both metamaterial unit cells presented double negative characteristic of ε and µ for almost the entire range of PD activity, with the exception of the bands corresponding to 1190-1350 MHz for the rectangular SRR and 1100-1300 MHz for the circular SRR, where the values of µ assume positive magnitudes. Therefore, both the SRR unit cells were classified as metamaterials applicable as superstrate for the designed bio-inspired PMA.
After the classification as metamaterials, the SRR unit cells were positioned as superstrates, as presented in Figure 9 and described in the Material and Methods section. The first analysis regarding the set PMA-Metamaterial is concerning the influence on the reflection coefficient due the application of the superstrate over the bio-inspired PMA. The results for this analysis are presented in Figure 21. From Figure 21, it can be verified that the superstrate application, with or without the presence of the metamaterial unit cells, resulted in shifts in the lower operating frequency of the bio-inspired PMA. For the FR4 superstrate (without metamaterial), a displacement of 42 MHz can be observed, resulting in a 5.2% operating bandwidth reduction regarding the reference bio-inspired PMA (487-1497 MHz). For the metamaterial unit cells, the displacements verified for the lower operating frequency were 56 and 77 MHz for the rectangular and circular SRRs, respectively. In addition, the metamaterial application also resulted in an increase of the upper frequency for both SRR geometries, 22 MHz for the rectangular and 64 MHz for the circular, resulting in a 3.9% and 1.8% operating bandwidth reduction, respectively, regarding the reference bio-inspired PMA.
Moreover, the application of superstrates resulted in higher reflection coefficient values regarding the reference bio-inspired PMA for the frequency range of 1050-1300 MHz, resulting in an entire operating bandwidth with reflection coefficient values under the -10 dB threshold. Therefore, the superstrate addition had a positive influence, since the coefficient reflection values were improved, and no significant bandwidth changes were observed.  From Figure 21, it can be verified that the superstrate application, with or without the presence of the metamaterial unit cells, resulted in shifts in the lower operating frequency of the bio-inspired PMA. For the FR4 superstrate (without metamaterial), a displacement of 42 MHz can be observed, resulting in a 5.2% operating bandwidth reduction regarding the reference bio-inspired PMA (487-1497 MHz). For the metamaterial unit cells, the displacements verified for the lower operating frequency were 56 and 77 MHz for the rectangular and circular SRRs, respectively. In addition, the metamaterial application also resulted in an increase of the upper frequency for both SRR geometries, 22 MHz for the rectangular and 64 MHz for the circular, resulting in a 3.9% and 1.8% operating bandwidth reduction, respectively, regarding the reference bio-inspired PMA.
Moreover, the application of superstrates resulted in higher reflection coefficient values regarding the reference bio-inspired PMA for the frequency range of 1050-1300 MHz, resulting in an entire operating bandwidth with reflection coefficient values under the −10 dB threshold. Therefore, the superstrate addition had a positive influence, since the coefficient reflection values were improved, and no significant bandwidth changes were observed.
In order to assess the superstrate's influence on the bio-inspired PMA radiation patterns, the lower, central, and upper frequency values were simulated and presented in Figures 22-24. In order to assess the superstrate's influence on the bio-inspired PMA radiation patterns, the lower, central, and upper frequency values were simulated and presented in Figure 22    From Figures 22-24, it can be verified that the superstrate addition without SRRs (only FR4) did not significantly change the radiation patterns regarding the reference bio-inspired PMA. However, radiation pattern distortions were verified for the upper operating frequencies for both SRR superstrates, as presented in Figure 24c,d.
Moreover, from the obtained radiation patterns, it can be verified that the antenna with metamaterial is suitable for dielectric window application, since the most of maximum radiation values are concentrated in the direction of coupling of the PMA to the dielectric window, i.e., 0° direction as presented in the Figure 25. In addition, considering that the PMA coupled to the dielectric window is back shielded against external interferences, the symmetrical back lobes presented in the radiation patterns also represent an advantage in the application for PD detection, since the back lobe allows the detection of the PD radiated signals reflected by the shielding, enhancing the PMA sensitivity of detection [26].  From Figures 22-24, it can be verified that the superstrate addition without SRRs (only FR4) did not significantly change the radiation patterns regarding the reference bio-inspired PMA. However, radiation pattern distortions were verified for the upper operating frequencies for both SRR superstrates, as presented in Figure 24c,d.
Moreover, from the obtained radiation patterns, it can be verified that the antenna with metamaterial is suitable for dielectric window application, since the most of maximum radiation values are concentrated in the direction of coupling of the PMA to the dielectric window, i.e., 0° direction as presented in the Figure 25. In addition, considering that the PMA coupled to the dielectric window is back shielded against external interferences, the symmetrical back lobes presented in the radiation patterns also represent an advantage in the application for PD detection, since the back lobe allows the detection of the PD radiated signals reflected by the shielding, enhancing the PMA sensitivity of detection [26]. From Figures 22-24, it can be verified that the superstrate addition without SRRs (only FR4) did not significantly change the radiation patterns regarding the reference bio-inspired PMA. However, radiation pattern distortions were verified for the upper operating frequencies for both SRR superstrates, as presented in Figure 24c,d.
Moreover, from the obtained radiation patterns, it can be verified that the antenna with metamaterial is suitable for dielectric window application, since the most of maximum radiation values are concentrated in the direction of coupling of the PMA to the dielectric window, i.e., 0 • direction as presented in the Figure 25. In addition, considering that the PMA coupled to the dielectric window is back shielded against external interferences, the symmetrical back lobes presented in the radiation patterns also represent an advantage in the application for PD detection, since the back lobe allows the detection of the PD radiated signals reflected by the shielding, enhancing the PMA sensitivity of detection [26].
Finally, from the extracted radiation patterns, the maximum gain values along the operating frequency for each simulated structure were obtained, as presented in Figure 26.
The application of both designed superstrates resulted in maximum gain values almost increased in double regarding the reference bio-inspired PMA for the frequency range of 500-1000 MHz. In addition, the circular SRR superstrate also demonstrated significant gain improvement for the frequency range of 1300-1500 MHz, probably due to the strong DNG feature presented for this frequency range, as presented in Figure 14.
Before the superstrate application, the reference bio-inspired PMA presented gain values higher than 2 dBi only for frequencies above 600 MHz. With the metamaterial addition, the gain values were higher than 2 dBi for the entire antenna operating frequency (500-1500 MHz).
The calculated mean gain values for each simulated structure are presented in Table 2.  Finally, from the extracted radiation patterns, the maximum gain values along the operating frequency for each simulated structure were obtained, as presented in Figure 26. Finally, from the extracted radiation patterns, the maximum gain values along the operating frequency for each simulated structure were obtained, as presented in Figure 26. The application of both designed superstrates resulted in maximum gain values almost increased in double regarding the reference bio-inspired PMA for the frequency range of 500-1000 MHz. In addition, the circular SRR superstrate also demonstrated significant gain improvement for the frequency range of 1300-1500 MHz, probably due to the strong DNG feature presented for this frequency range, as presented in Figure 14a.  Table 2, it can be noticed that the metamaterial application as superstrate resulted in a maximum mean gain improvement higher than 1 dBi, regarding the reference bio-inspired PMA, for both SRR geometries chosen. Although the circular SRR structure presented a slightly higher mean gain value (4.07 dBi) in comparison to the rectangular geometry (3.96 dBi), the reflection coefficient results for the rectangular SRR presented lower bandwidth impacts for PD application than the circular one. In addition, the rectangular SRR presented lower radiation patterns distortions than the circular one for the upper frequencies, mainly in the back lobe direction, as presented in Figures 24 and 25. Therefore, since symmetrical back lobes are desirable for PD detection through dielectric windows, the rectangular SRR was considered the fittest structure for the experimental PD tests presented in the following section.

Laboratory Experimental Results
As mentioned in the previous subsection, the simulated reflection coefficients, radiation patterns, and mean gain appointed the rectangular SRR metamaterial superstrate as the fittest structure for the experimental PD tests. Therefore, from the models presented in Figures 6b and 9a, the bio-inspired PMA without and with the metamaterial superstrate were built, as presented in Figure 27.  Table 2, it can be noticed that the metamaterial application as superstrate resulted in a maximum mean gain improvement higher than 1 dBi, regarding the reference bio-inspired PMA, for both SRR geometries chosen. Although the circular SRR structure presented a slightly higher mean gain value (4.07 dBi) in comparison to the rectangular geometry (3.96 dBi), the reflection coefficient results for the rectangular SRR presented lower bandwidth impacts for PD application than the circular one. In addition, the rectangular SRR presented lower radiation patterns distortions than the circular one for the upper frequencies, mainly in the back lobe direction, as presented in Figure 24; Figure 25. Therefore, since symmetrical back lobes are desirable for PD detection through dielectric windows, the rectangular SRR was considered the fittest structure for the experimental PD tests presented in the following section.

Laboratory Experimental Results
As mentioned in the previous subsection, the simulated reflection coefficients, radiation patterns, and mean gain appointed the rectangular SRR metamaterial superstrate as the fittest structure for the experimental PD tests. Therefore, from the models presented in Figure 6b and Figure 9a, the bio-inspired PMA without and with the metamaterial superstrate were built, as presented in Figure 27. The first experimental analysis for the manufactured antennas is regarding the reflection coefficients. In Figure 28 a comparison between the measured and simulated reflection coefficients for the bio-inspired PMA with and without the rectangular SRR metamaterial superstrate is presented. The first experimental analysis for the manufactured antennas is regarding the reflection coefficients. In Figure 28 a comparison between the measured and simulated reflection coefficients for the bio-inspired PMA with and without the rectangular SRR metamaterial superstrate is presented. Both measured and simulated reflection coefficient results presented similar behavior. For the bio-inspired PMA without metamaterial, the measured operating bandwidth started at 472 MHz up to values above the upper PD frequency of interest (1500 MHz). With the metamaterial, the lower frequency was displaced in 68 MHz (close to the 77 MHz simulated value) regarding the reference bio-inspired PMA, starting from 540 MHz until values above 1500 MHz. In addition, better reflection coefficient performance for the frequency range of 1.1-1.5 GHz was verified for the measured curve with metamaterial, as predicted by the simulations. The slight differences between the simulated and measured reflection coefficient behaviors can be attributed to PMA constructive aspects, such as minor discrepancies between the simulated FR4 and the one used in practice, as well the non-consideration of the SMA connectors during the simulation.
Through the procedures described in Section 4, the bio-inspired PMA gains with and without the metamaterial superstrates were measured. The comparison between the measured and simulated gains are presented in Figure 29. Both measured and simulated reflection coefficient results presented similar behavior. For the bio-inspired PMA without metamaterial, the measured operating bandwidth started at 472 MHz up to values above the upper PD frequency of interest (1500 MHz). With the metamaterial, the lower frequency was displaced in 68 MHz (close to the 77 MHz simulated value) regarding the reference bio-inspired PMA, starting from 540 MHz until values above 1500 MHz. In addition, better reflection coefficient performance for the frequency range of 1.1-1.5 GHz was verified for the measured curve with metamaterial, as predicted by the simulations. The slight differences between the simulated and measured reflection coefficient behaviors can be attributed to PMA constructive aspects, such as minor discrepancies between the simulated FR4 and the one used in practice, as well the non-consideration of the SMA connectors during the simulation.
Through the procedures described in the Material and Methods section, the bio-inspired PMA gains with and without the metamaterial superstrates were measured. The comparison between the measured and simulated gains are presented in Figure 29. From Figure 29, it can be observed that, although there is a discrepancy in the gain values, the behavior of the measured and simulated curves are in agreement.
For both measured and simulated curves for the bio-inpired PMA without metamaterial, an increase in the gain values up to frequencies around 1.1 GHz can be verified, followed by a gain decrease up to 1.2 GHz and, again, a gain increase up to 1.4 GHz, followed by a new decrease   From Figure 29, it can be observed that, although there is a discrepancy in the gain values, the behavior of the measured and simulated curves are in agreement.
For both measured and simulated curves for the bio-inpired PMA without metamaterial, an increase in the gain values up to frequencies around 1.1 GHz can be verified, followed by a gain decrease up to 1.2 GHz and, again, a gain increase up to 1.4 GHz, followed by a new decrease tendency up to frequency values above 1.5 GHz. For both measured and simulated curves for the bio-inspired PMA with metamaterial, the same trend was verified, except for the frequencies above 1.4 GHz, in which the bio-inspired PMA with metamaterial presented a gain increase trend up to frequencies above 1.5 Ghz.
Besides the constructive aspects mentioned before, the differences between the measured and simulated values can be attributed to the conception of the simulated environment, which neglects some practical limitations faced during the measurements, such as the cables attenuation, the alignment between the reference antenna and the test antenna, and the losses in the far-field measuring system.
Even so, the gain enhancement provided by the metamaterial superstrate was according to that predicted by the simulations, presenting significant gain improvement for mostly the bandwidth of the bio-inspired PMA. The measured mean gain for the bio-inspired PMA without metamaterial was equal to 2.92 dBi, close to the simulated mean gain value presented in Table 2 (2.86 dBi). In addition, the bio-inspired PMA without metamaterial only presented gain values above the 2 dBi threshold for the frequencies above 700 MHz. With the metamaterial superstrate, the calculated mean gain was 3.61 dBi, resulting in an enhancement of 0.7 dBi regarding the bio-inspired PMA. Moreover, the application of the metamaterial superstrate resulted in an antenna with gain values above 2 dBi from 575 MHz up to 1.5 GHz, resulting in an antenna with high sensitivity for PD detection.
From the measured reflection coefficient and gain results presented, the bio-inspired PMA-metamaterial set can be considered fit for PD detection. Therefore, it was submitted to high voltage PD detection sensitivity tests.
Through the setup presented in Figure 13, the PD activity was detected simultaneously by the bio-inspired PMA-metamaterial set and the IEC 60270 standard method for the application of 13.4 kV on the oil cell. Figure 30 presents one sample of the simultaneously detected PD pulses.
predicted by the simulations, presenting significant gain improvement for mostly the bandwidth of the bio-inspired PMA. The measured mean gain for the bio-inspired PMA without metamaterial was equal to 2.92 dBi, close to the simulated mean gain value presented in Table 2 (2.86 dBi). In addition, the bio-inspired PMA without metamaterial only presented gain values above the 2 dBi threshold for the frequencies above 700 MHz. With the metamaterial superstrate, the calculated mean gain was 3.61 dBi, resulting in an enhancement of 0.7 dBi regarding the bio-inspired PMA. Moreover, the application of the metamaterial superstrate resulted in an antenna with gain values above 2 dBi from 575 MHz up to 1.5 GHz, resulting in an antenna with high sensitivity for PD detection.
From the measured reflection coefficient and gain results presented, the bio-inspired PMA-metamaterial set can be considered fit for PD detection. Therefore, it was submitted to high voltage PD detection sensitivity tests.
Through the setup presented in Figure 13, the PD activity was detected simultaneously by the bio-inspired PMA-metamaterial set and the IEC 60270 standard method for the application of 13.4 kV on the oil cell. Figure 30 presents one sample of the simultaneously detected PD pulses. From Figure 30, it can be verified that the bio-inspired PMA-metamaterial set was able to detect all the PD pulses detected by the IEC 60270 standard method. In addition, the voltage measured at the antenna's terminals presented values reduced by one-half compared to the standard method. This result is expected, since the radiated waves are susceptible to higher losses, such as reflections and propagating attenuations, than those in the standard method, directly connected with the oil cell. Nonetheless, the antenna sensitivity could be attested, since, according to Table 1, the detected PD pulses presented apparent charge values between 15 and 60 pC. This range of apparent charges is similar to the pC values detected in practice for the PD low intensity pulses related to the inception of an insulating problem in transformers oil, for example in [57][58][59].
In order to evaluate the frequency range of the detected PD pulses, a spectrogram was extracted from the PMA detected PD pulses shown in Figure 30 and presented in Figure 31.
From Figure 31, it can be verified that the four PD pulses observed in Figure 23 were highlighted in the spectrogram, presenting higher concentrations of spectral energy for frequencies up to 1500 MHz. The obtained frequency range is in agreement with the values reported in [9,10,60] as the main range of interest for the monitoring of partial discharges in oil (300-1500 MHz). and propagating attenuations, than those in the standard method, directly connected with the oil cell. Nonetheless, the antenna sensitivity could be attested, since, according to Table 1, the detected PD pulses presented apparent charge values between 15 and 60 pC. This range of apparent charges is similar to the pC values detected in practice for the PD low intensity pulses related to the inception of an insulating problem in transformers oil, for example in [57][58][59].
In order to evaluate the frequency range of the detec Figure 31. Spectrogram of the PD pulses sample detected by the bio-inspired PMA-metamaterial set presented in Figure 30.
From Figure 31, it can be verified that the four PD pulses observed in Figure 23 were highlighted in the spectrogram, presenting higher concentrations of spectral energy for frequencies up to 1500 MHz. The obtained frequency range is in agreement with the values reported in [9,10,60] as the main range of interest for the monitoring of partial discharges in oil (300-1500 MHz).
Therefore, besides the sensitivity (gain), the antenna operating bandwidth was also attested in the high voltage PD tests, since the bio-inspired PMA-metamaterial set was able to successfully detect all the main frequency range of energy concentration for PD pulses.
Hence, the bio-inspired PMA-metamaterial set was approved in the laboratory PD experimental tests, being able to be applied in PD measurements on substations, as presented in the following subsection.

PD Detection in a 230 kV Substation
Applying the methodology presented in Subsection 4.6, PD measurements on a 230 kV CT were performed and the comparison between the PD pulses detected by the bio-inspired PMA-metamaterial set and by the reference bio-inspired PMA are presented in Figure 32. Therefore, besides the sensitivity (gain), the antenna operating bandwidth was also attested in the high voltage PD tests, since the bio-inspired PMA-metamaterial set was able to successfully detect all the main frequency range of energy concentration for PD pulses.
Hence, the bio-inspired PMA-metamaterial set was approved in the laboratory PD experimental tests, being able to be applied in PD measurements on substations, as presented in the following subsection.

PD Detection in a 230 kV Substation
Applying the methodology presented in Section 4.6, PD measurements on a 230 kV CT were performed and the comparison between the PD pulses detected by the bio-inspired PMA-metamaterial set and by the reference bio-inspired PMA are presented in Figure 32. From Figure 32, it can be verified that the detected signals presented PD features, since they occurred in the half-cycle transitions of the sinusoidal reference voltage and with a phase displacement approximately equal to 180° between each other. In addition, the bio-inspired From Figure 32, it can be verified that the detected signals presented PD features, since they occurred in the half-cycle transitions of the sinusoidal reference voltage and with a phase displacement approximately equal to 180 • between each other. In addition, the bio-inspired PMA-metamaterial set presented higher detection sensitivity than the reference bio-inspired PMA for all of the detected PD pulses, presenting a signal with amplitudes almost doubled regarding the reference bio-inspired PMA.
As expected, according to previous measurements performed by the substation maintenance personnel, a considerable level of PD was detected in the 230 kV CT, attesting the sensitivity of the developed antenna in the practical detection of PD.
Although the PMA presented high sensitivity for PD detection, it was not possible to estimate the PD apparent charge levels inside the equipment and detected by the antennas, since the calibration procedure on the 230 kV CT could not be performed due to practical limitations (coupling capacitor parallel connection). Therefore, at this point of the research, the CT operational condition diagnosis (low, medium, or critical level of PD) could not be estimated.
Nonetheless, the application of the developed bio-inspired PMA-metamaterial set in the CT PD measurement was done in order to attest its practical sensitivity in the PD detection. Hence, as future works, calibration, classification, and localization techniques will be applied over the proposed antenna in order to fully characterize it as a sensor applicable in the PD continuous monitoring and diagnosis, specially through dielectric windows improved by the application of metamaterials on its surface.

Discussion
Although there are more competitive antenna types that can provide higher gain than PMA, such as hornet and Vivaldi antennas, their use in PD detection is compromised due to their constructive aspects (mainly their length), since they are not suitable for installation in dielectric windows in power transformers or GIS. Despite the fact that PMA are not as efficient as another types of antennas, they still represent a good option for PD detection, since its large bandwidth, gain, radiation pattern and dimensions are still suitable for PD detection, as verified and shown on other studies in the literature [25][26][27][61][62][63][64][65][66].
As presented in this paper, the distance between the antenna and the metamaterial superstrate (dielectric window) was adopted in order to reduce the impact on the antenna's reflection coefficient. Still, even with the distance between the bio-inspired antenna and the metamaterial superstrate (dielectric window), the final structure has a volumetric dimension of 200 mm × 200 mm × 135 mm to accommodate an antenna with an operating bandwidth of 500-1500 MHz and average gain of 3.61 dBi, which represents an acceptable size for PD application and good PD detection sensitivity features as demonstrated experimentally during the paper and through comparisons with other results obtained in the literature [67][68][69][70][71][72][73]. If another type of antennas were used for this same bandwidth, their higher length would require a much larger dielectric window, becoming an impractical application. The dimensions of some more directional antennas (commercial and from other papers) with lower operating frequency equivalent to the proposed antenna (500 MHz) are presented in the Table 3. As can be seen from Table 3, although the good results of bandwidth and gain for PD detection application, all the other antennas present a much larger length than the proposed PMA-dielectric window set (135 mm). Moreover, by coupling the presented antennas over dielectric windows, we would have even larger dimensions for the antenna-dielectric window set than those presented in Table 3, becoming impractical for the application on PD detection in power transformers or GIS. Therefore, this reinforces the advantage of using PMA for the achievement of dielectric windows with smaller dimensions and higher ease of coupling on the structure of high voltage equipment.

Conclusions
In this paper, the representation of dielectric windows as a metamaterial superstrate over a bio-inspired antenna (based on the Jatropha mollissima (Pohl) Baill leaf) was proposed in order to improve the detection sensitivity of UHF sensors designed for PD measurement. During the simulation tests, both the selected SRR shapes (rectangular and circular) proved themselves efficient in the enhancement of the antenna detection sensitivity, presenting a mean gain improvement higher than 1 dBi compared to the structure without metamaterial. However, the rectangular SRR was selected for the laboratory tests due to the better relation gain/bandwidth/radiation pattern achieved. The designed bio-inspired PMA-metamaterial set presented potential to be used as a UHF sensor applied through dielectric windows for PD detection, since it meets the bandwidth (300-1500 MHz) and gain (>2 dBi) requirements according to the laboratory tests carried out in an anechoic chamber. This potential was attested in the PD detection tests, in which the developed sensor was able to detect several PD signals generated in laboratory, by the point-to-plane electrode configuration, and in practice, by a 230 kV CT. Although the PD apparent charge levels could not be estimated in practice through the signals detected by the antenna, the evaluation of the induced voltage in the antenna terminals can provide information about the PD activity evolution inside the equipment. Therefore, the proposed bio-inspired PMA-Metamaterial set can provide continuous, non-invasive, and relatively low-cost monitoring of PD activity by means of dielectric windows in high voltage insulating systems.
The conclusions obtained from the results of this article can be summarized as follows: (1) Both the SRR-CLS designed unit cells presented the double negative behavior (µ < 0, ε < 0) for almost all the main frequency range of PD activity (300-1500 MHz). Therefore, both the SRR unit cells were classified as metamaterials applicable as superstrate for the designed bio-inspired PMA; (2) The application of a metamaterial superstrate resulted in displacements in the bio-inspired PMA operating frequency range, resulting in a 3.9% and 1.8% bandwidth reduction for the circular and rectangular SRRs, respectively; (3) Although there was a bandwidth reduction regarding the bio-inspired PMA, the application of the both metamaterial superstrates resulted in higher reflection coefficient values for the frequency range of 1100-1500 MHz, in which some reflection coefficient values were lower than the threshold of −10 dB, improving the antenna power transmission/reception capacity for this frequency range; (4) Radiation patterns distortions were observed at the upper bio-inspired PMA operating frequencies for both metamaterial superstrates, in which the rectangular SRR presented lower radiation patterns distortions than the circular one for the upper frequencies, mainly at the back lobe direction; (5) The application of both designed superstrates resulted in maximum gain values almost increased in double regarding the reference bio-inspired PMA for the frequency range of 500-1000 MHz, resulting in a mean gain improvement higher than 1 dBi; (6) The practical application of the rectangular SRR superstrate resulted in a UHF sensor with an operating bandwidth (540-1500 MHz) that covers 80% of the main frequency range of PD activity (300-1500 MHz); (7) The measured mean gain for the bio-inspired PMA-metamaterial set was equal to 3.61 dBi, resulting in an enhancement of 0.7 dBi regarding the bio-inspired PMA (2.92 dBi) and representing a good detection sensitivity for PD application;