Main Content

Metasurface Antenna Modeling

This example shows how to design, model, and analyze a metasurface antenna based on reference [1]. Metasurface antennas take advantage of the periodic boundary of each radiation unit cell to group the radiation unit in a more compact space.

You construct an unit radiation element with Antenna Toolbox™ objects and functions. You also use infiniteArray and rectangularArray objects to construct the large antenna array, on which you perform pattern and current analysis at the designed frequency.

Construct Unit Radiator

The unit radiator in [1] is a cross I-beam structure with an optimized probe feed in one arm. The unit cell has a dimension of 3 cm and gap of 0.5 mm. Specify the geometry of the antenna.

length = 14.75*1e-3;
width = 3*1e-3;
s = 0.25*1e-3;
thickness = 1.5*1e-3;
viaDia = 0.5*1e-3*2;
feedWidth = viaDia/2;  
gndLength = length + 2*s;
gndWidth = gndLength;

Use the customAntennaGeometry object to construct the top cross I-beam radiator and place it into a reflector to form the unit radiator structure. In this simplified geometry, you remove the substrate from the design. Removing the substrate increases the radiation frequency to 9.14 GHz, which is higher than the value in [1].

s1 = antenna.Rectangle(Length=length, Width=width);
s2 = antenna.Rectangle(Length=width, Width=length);
s3 = antenna.Rectangle(Length=feedWidth, Width=feedWidth);
[~] = translate(s3,[-5*1e-3+feedWidth/2 0 0]);

pr1 = getShapeVertices(s1);
pr2 = getShapeVertices(s2);
pr3 = getShapeVertices(s3);

radiator = customAntennaGeometry(Boundary={pr1,pr2,pr3}, Operation="P1+P2+P3");
radiator.FeedLocation = [-5*1e-3 0 0];
radiator.FeedWidth = feedWidth;

ant = reflector(Exciter=radiator, GroundPlaneLength=gndLength, GroundPlaneWidth=gndWidth, ...
    Spacing=thickness, EnableProbeFeed=true);

figure
show(ant)

Figure contains an axes object. The axes object with title reflector antenna element, xlabel x (mm), ylabel y (mm) contains 5 objects of type patch, surface. These objects represent PEC, feed.

Mesh and Analyze Unit Structure

Generate a mesh with a maximum edge length of 0.1 m.

figure;
mesh(ant,MaxEdgeLength=0.1);

Figure contains an axes object and an object of type uicontrol. The axes object with title Metal mesh, xlabel x (m), ylabel y (m) contains 2 objects of type patch, surface. These objects represent PEC, feed.

Calculate the impedance as a function of frequency.

freq = 9.14e9;
figure;
impedance(ant,freq*[0.8:0.01:1.1]);

Figure contains an axes object. The axes object with title Impedance, xlabel Frequency (GHz), ylabel Impedance (ohms) contains 2 objects of type line. These objects represent Resistance, Reactance.

Implement Periodic Boundary and Analyze Antenna

To include the periodic boundary effect on the unit radiator, specify the radiator as the element in the infiniteArray object. The object constructs an infinite antenna array with the designed radiator. Numerically, the object uses a Green function with a method of moments (MOM) and models the coupling effect between adjacent unit radiators.

infArray = infiniteArray(Element=ant);
figure
show(infArray)

Figure contains an axes object. The axes object with title Unit cell of reflector in an infinite array, xlabel x (mm), ylabel y (mm) contains 7 objects of type patch, surface. These objects represent PEC, feed, Air, unit cell.

figure
current(infArray,freq)

Figure contains an axes object. The axes object with title Current distribution, xlabel x (m), ylabel y (m) contains 5 objects of type patch.

Analyze Finite Array

To analyze the coupling effect on a finite array, construct a 5-by-5 element rectangular array with the designed radiator.

array = rectangularArray(Element=ant);
array.Size = [5 5];
array.RowSpacing = gndLength;
array.ColumnSpacing = gndWidth;
figure
show(array)

Figure contains an axes object. The axes object with title rectangularArray of reflector antennas, xlabel x (mm), ylabel y (mm) contains 102 objects of type patch, surface. These objects represent PEC, feed.

Plot the pattern and current of the array as functions of the frequency.

figure
pattern(array,freq)

Figure contains 2 axes objects and other objects of type uicontrol. Axes object 1 contains 102 objects of type patch, surface. Hidden axes object 2 contains 42 objects of type surface, line, text, patch.

figure
current(array,freq)

Figure contains an axes object. The axes object with title Current distribution, xlabel x (m), ylabel y (m) contains 78 objects of type patch.

Assign Substrate to Infinite Array

Specify Teflon™ as the substrate of the unit cell element and visualize the geometry of the infinite array.

infArray.Element.Substrate = dielectric("Teflon");
figure
show(infArray) 

Figure contains an axes object. The axes object with title Unit cell of reflector in an infinite array, xlabel x (mm), ylabel y (mm) contains 7 objects of type patch, surface. These objects represent PEC, feed, Teflon, unit cell.

Mesh Unit Cell

Mesh the unit cell, specifying a maximum edge length of lambda/40, where lambda is the free-space wavelength at 9.14 GHz. The software captures the dielectric effect in the infinite array analysis in the equivalent periodic Green's function, so use the MOM over the top layer metal mesh only.

lambda = 3e8/9.14e9;
figure
mesh(infArray,MaxEdgeLength=lambda/40) 

Figure contains an axes object and an object of type uicontrol. The axes object with title Metal mesh, xlabel x (m), ylabel y (m) contains 2 objects of type patch, surface. These objects represent PEC, feed.

Compute Pattern

With the meshed geometry, perform the full-wave radiation pattern analysis at 9.14 GHz.

figure
pattern(infArray,9.14e9);

Figure contains 2 axes objects and other objects of type uicontrol. Axes object 1 contains 6 objects of type patch, surface. This object represents Teflon. Hidden axes object 2 contains 17 objects of type surface, line, text, patch. This object represents Teflon.

Conclusion

The current and pattern modeling result obtained from the finite array are qualitatively similar to the results reported in [1]. For other types of metasurface geometry, you can construct the unit cell by using an antenna catalog element, as the Infinite Array of Microstrip Patch Antenna on Teflon Substrate example shows, or using a PCB stack.

References

[1] Badawe, Mohamed El, Thamer S. Almoneef, and Omar M. Ramahi. “A True Metasurface Antenna.” Scientific Reports 6, no. 1 (January 13, 2016): 19268. https://doi.org/10.1038/srep19268.

See Also

Related Topics