International

OPEN ACCESS J ournal

Of Modern Engineering Research (IJ MER)

| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 24 |
Numerical Investigation of Multilayer Fractal FSS

R. H. C. Manioba
1
, A. F. Santos
2
, A. V. Lovato
3
, N. M. Oliveira-Neto
4
,
D. B. Brito
5
, A. L. P. S. Campos
6
, A. G. dAssuno
7

1, 2, 3, 4
(Collegiate of Information Systems Course, State University of Southwest Bahia UESB, Jequi-BA,
Brazil)
5
(Institute of Computing, Federal University of Alagoas, Macei-AL, Brazil)
6, 7
(Communication Engineering Department, Federal University of Rio Grande do Norte, Natal-RN, Brazil)









I. INTRODUCTION
A periodic surface can be defined basically as a set of identical elements arranged in two or three
dimensions forming an infinite array. A periodic array formed by conductive patch elements or opening
elements is known as Frequency Selective Surface (FSS). FSS are periodic structures that can provide
frequency filtering to incoming electromagnetic waves, and their frequency response is entirely determined by
the geometry of the structure in one period called a unit cell [1, 2].
Traditional FSS has a long history of development and have been investigated over the years for a
variety applications, e.g., frequency filters or diplexers in high performance reflector antenna systems,
advanced radome designs, and smart surfaces for stealth applications [3].
The question of operating bandwidth is one of important problems in FSS theory, in some applications
a multiband frequency response is desired; however in other applications a wideband frequency response is
more preferred. FSS structures have been successfully proven as a mean to increase the communication
capabilities of satellite applications. In space mission (Voyager, Galileo and Cassini), for example, the use of
dual antenna system with FSS reflectors has made it possible to share the main reflector among different
frequency bands, therefore FSS with multiband responses have been studied by several researchers. But, in
some applications, such as separating successive frequency bands, FSS are required to have a wide
transmission or reflection bands and rapid rolloffs [4-10].
Experimental and numerical investigations on cascading or multilayer FSS have been done, using full
and simple (approximate methods) [11-14]. This paper presents a numerical investigation of multilayer or
cascaded FSS using a simple but efficient method used to predict the frequency response (transmission
coefficient) for these structures. In this work was analyzed a multilayer FSS structure, using conducting Koch
fractal elements. This structure was first proposed in [15] and consists of two FSS, called structure 1 and
structure 2 mounted on a dielectric isotropic layer, separated by an air gap layer. The dielectric substrate used
was the RT-Duroid 3010, with 1.27 mm of height and relative permittivity equal to 10.2 and unit cell
periodicity (T
x
= T
y
) equal to 10 mm. We used the commercial software Ansoft Designer
TM
to obtain the
individual scattering parameters of each one FSS structure. To validate the results a comparison with results
from other numerical method is performed.

II. FSS THEORY
Through the use of the spectral domain approach to analyzing the response of an FSS, which is
assumed to have infinitesimal thickness, one can determine the electromagnetic fields on a plane z = z
n
, see
Fig. 1 and Fig. 2, These electric fields are expressed in terms of a discrete spectrum of plane waves known as
Floquet harmonics.
Abstract: Numerical investigations are presented for a multilayer frequency selective surface with Koch
fractal (levels 1 and 2) conducting patch elements. The structure investigated is obtained using two FSS
screens separated by an air gap layer. For the proposed investigation were used three different values an
air gap height. The results obtained using the numerical method were compared with other technique and
using the commercial software Ansoft Designer
TM
. A good agreement was observed in terms of the
bandwidth.
Keywords: FSS, Numerical Method, Multilayer Structure, Fractal, Wideband

Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 25 |

Fig. 1: Front (superior) view of a FSS.


Fig. 2: FSS Cascading, side view.

Cwik in [16] has defined a set of scattering parameters that incorporates the vector nature of electromagnetic
fields as follows [1]:
( , )
11
( , , , , )
( , , , )
( , )
xm yn xi yj L
inc
V k k k k z
S m n i j
P i j
+
=
( , )
12
( , , , , )
( , , , ) ,
( , )
xm yn xi yj R
inc
V k k k k z
S m n i j
P i j

=
( , )
21
( , , , , )
( , , , ) ,
( , )
xm yn xi yj L
inc
V k k k k z
S m n i j
P i j
+ +
=
( , )
22
( , , , , )
( , , , ) ,
( , )
xm yn xi yj R
inc
V k k k k z
S m n i j
P i j
+
=
where z
R
is the above interface and z
L
is the below interface in Fig. 1, ie, z
R
= 0 and z
L
= z
p+q
, V
(,)
are Floquet
voltage waves given by [1]:
( , ) ( , )
( , , , , ) ( , , , , ) ( , ),
xm yn xi yj R
V k k k k z f m n i j z P m n

=


Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 26 |
With,
2 2
( , ) ( ) ,
xm yn mn
P m n k k Y = +
Where:
,
mn
mn
Y
j

e
=

The terms K
xm
and K
yn
(wave numbers of Floquet special harmonics) are defined as:
2
,
inc
xm x
x
k m k
T
t
= +
2
,
inc
yn y
y
k n k
T
t
= +
The incident harmonics are functions of elevation angle and the azimuthal angle, which are given as [1]:
0
cos ,
inc
x inc inc
k k senu | =
0
,
inc
y inc inc
k k sen sen u | =
Referring to the coordinate system of Fig. 1, can be seen that the notation (, ) associated with the
definition for the voltage waves indicates the direction, both +z or -z, of the incident and scattered energy in a
FSS system. Regarding the scattering parameters, it is more convenient to express the potential vector, used in
the definition of the Floquet voltage waves in terms of the electric fields transform.
The potential vector, expressed in terms of total electric fields, can be written as [1]:

( , )
2 2
( )
( , , , , ) ,
total total
yn x xm y
xm yn
j k E k E
f m n i j z
k k



=
+

The total electric field can be found as follows [1]:
( , ) ( , ) ( , )
( , ) ( ) ( ),
total S ref
x y xm yn x y x y xi xm yj yn
E k k E E k k k k o o
+ +
= +


( , ) ( , ) ( , )
( , ) ( ) ( ),
total S trans
x y xm yn x y x y xi xm yj yn
E k k E E k k k k o o

= +



Where E
(ref, trans)
are the reflected and transmitted electric fields, calculated in z = 0 and z = z
p+q
. We can write
the scattering parameters in terms of the electric field evaluated at the reference planes above (z
R
) and below
(z
L
) of a FSS:
11
2 2 2 2
( )
( , , , ) ,
( )( )
total total
yn x xm y
mn
ij
xm yn xi yj
j k E k E
Y
S m n i j
Y
k k k k
+

=
+ +


12
2 2 2 2
( )
( , , , ) ,
( )( )
total total
yn x xm y
mn
ij
xm yn xi yj
j k E k E
Y
S m n i j
Y
k k k k

=
+ +


21
2 2 2 2
( )
( , , , ) ,
( )( )
total total
yn x xm y
mn
ij
xm yn xi yj
j k E k E
Y
S m n i j
Y
k k k k
+ +

=
+ +


22
2 2 2 2
( )
( , , , ) ,
( )( )
total total
yn x xm y
mn
ij
xm yn xi yj
j k E k E
Y
S m n i j
Y
k k k k
+

=
+ +


Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 27 |
III. NUMERICAL METHOD
Since scattering parameters (scattering matrix) of finite dimension for a FSS are computed, many
analytical procedures are available for obtaining multilayer composite representation. It is possible to make
direct use of the scattering matrix along with the following equation to obtain the representation of a system
with two FSS (multilayer or cascaded FSS), each FSS (individual structure) is viewed as a subsystem [1]. The
scattering matrix for the multilayer FSS, using two FSS, is given by [17]:

(1) (1) (2) (1) (1) (2)
11 12 11 21 12 12
(2) (1) (1) (2) (1) (2) (1) (2) (2)
21 21 22 11 21 21 22 12 22
,
( )
C
S S PTS PS S PTS
S
S P S S PTS PS S PS PTS S
( +
=
(
+ +



Where
) 1 (
11
S ,
) 1 (
12
S ,
) 1 (
21
S and
) 1 (
22
S are the scattering parameters which represent the first subsystem, namely,
the first FSS,
) 2 (
11
S ,
) 2 (
12
S ,
) 2 (
21
S and
) 2 (
22
S are the scattering parameters associated with the second subsystem,
and
C
S is the scattering matrix for the system composed by two FSS (multilayer FSS).
And T is given by [17]:
(2) (1) 1
11 22
( ) , T I S PS P

=

P is a diagonal matrix that has
) 2 ( kd j
e

as its elements, k is the wave number [45], in this work k = k
0
(because
the multilayer FSS is formed by two single FSS separated by an air gap layer), and I is an identity matrix.
To validate the results, that will be presented, a comparison with results from other numerical method called
One Mode Interaction Technique that can be seen in [17] and obtained by Method of Moments (MoM)
using the commercial software Ansoft Designer
TM
is performed. The numerical results were computed using
the commercial software MATLAB
TM
(MATrix LABoratory).
Using the One Mode Interaction (OMI) technique, the transmission (C
T
) and reflection (C
R
) coefficients for the
multilayer FSS, in this case formed by two single FSS, are [18]:

1
1 2
( 2 )
1 2
,
1
T j kd
TT
C
R R e

=

1
1
2
( 2 ) 1 2
1 ( 2 )
1 2
,
1
j kd
R j kd
T R
C R e
R R e

= +

Where d
1
is the spacing between the two structures. T
1
, R
1
are the transmission and reflection coefficients for
the first FSS, and T
2
, R
2
are the transmission and reflection coefficients for the second FSS.

IV. MULTILAYER FSS STRUCTURE
The multilayer FSS used in this work consists of two FSS screens called structure 1 and structure 2,
respectively, each one using Koch fractal patch elements printed on a dielectric substrate separated by an air
gap layer. The first structure with Koch fractal level 1 is defined as structure 1 and the second FSS screen with
Koch fractal level 2 is defined as structure 2, as can be seen in Table 1 with their respective resonant
frequencies (f) and bandwidths (BW). The process to generate the Koch fractal elements used in this work can
be seen in [14]. The element shapes and the multilayer structure considered in our investigation are shown in
Fig. 3 and Fig. 4, respectively.

Table 1 Identification of the elements
Structure Element type f (GHz) BW (GHz)
1
Koch fractal level
1
9.115
1.951
2
Koch fractal level
2
8.333
1.618
Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 28 |

Fig. 3: Koch fractal shape elements.


Fig. 4: Multilayer Structure.

V. RESULTS
In order to validate the results obtained with the numerical method presented before, they are
compared with results obtained with the OMI technique and the results obtained using the commercial
software Ansoft Designer
TM
, the latter uses the full wave method, MoM, to compute its results.
Fig. 5 illustrates the results obtained with the numerical method for the case with an air gap height
equal to 2.0 mm. For the multilayer FSS with an air gap height equal to 2.0 mm, the results obtained using the
numerical method presents a bandwidth, for 20 dB insertion loss reference, equals to 3.420 GHz, a good
agreement is observed between the results.
For the multilayer FSS with an air gap height equal to 4.0 mm, the results obtained using the
numerical method presents a bandwidth, for 20 dB insertion loss reference, equal to 3.360 GHz,
approximately. Fig 6. illustrates a comparison between the results.
Using an air gap height equal to 6.0mm, the numerical method presents a bandwidth, for the 20 dB insertion
loss reference, equal to 3.290 GHz, approximately.

Fig. 5: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 2.0 mm.
Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 29 |

Fig. 6: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 4.0 mm.

Fig. 7: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 6.0 mm.

VI. CONCLUSION
In this paper, a multilayer FSS was investigated using a simple numerical method. This multilayer
structure was formed by two FSS screens separated by an air gap layer, each one FSS screen using a
conducting patch element with fractal geometry in the unit cell. The results were obtained for the numerical
method and compared with other numerical technique called One Mode Interaction. Moreover, both were
compared with the results obtained using the commercial software Ansoft Designer
TM
for different values of
spacing between the FSS screens. The numerical method is efficient and can be used in conjunction with other
methods, like full wave methods such like: Finite Element Method, Wave Concept Iterative Procedure, and
Equivalent Circuit Model, among others. For all the cases considered, the frequency responses have
characteristics of a high-pass filter with a very large rejection band. A good agreement between the results was
observed in terms of bandwidth.

REFERENCES
[1] T. K. Wu, Frequency selective surface and grid array (Jonh Wiley & Sons, New York, E.U.A., 1995).
[2] B. A. Munk, Frequency Selective Surfaces Theory and Design (Jonh Wiley & Sons, New York, E.U.A., 2000).
[3] F. Bayatpur and K. Sarabandi, Single-Layer High-Order Miniaturized-Element Frequency-Selective Surfaces, IEEE
Transactions on Microwave Theory and Techniques, 56 (4), 2008, 774-781.
[4] J. Romeu and Y. R. Samii, Fractal FSS: A novel dual-band Frequency selective surface, IEEE Transactions on
Antennas and Propagation, 48 (7), 2000, 1097-1105.
Numerical Investigation of Multilayer Fractal FSS
| IJMER | ISSN: 22496645 | www.ijmer.com | Vol. 4 | Iss. 5| May. 2014 | 30 |
[5] J. P. Gianvittorio, Y. Rahmat-Samii and J. Romeu, Fractal FSS: Various Self-Similar Geometries Used for Dual-
Band and Dual-Polarized FSS, IEEE Antennas and Propagation International Symposium, Boston, MA, 2001, 640-
643.
[6] J. C. Zhang; Y. Z. Yin, and S. F. Zheng, Double Screen FSSs with Multi-Resonant Elements for Multiband,
Broadband Applications, Journal of Electromagnetic Waves and Applications, 23 (16), 2009, 2209-2218.
[7] A. L. P. S. Campos, T. L. Silva and A. G. Neto, Multiband Frequency Selective Surfaces with Simple Modification
of a Rectangular Patch Element, Microwave and Optical Letters, 55 (12), 2013, 2206-2209.
[8] F. C. G. S. Segundo, A. L. P. S. Campos and A. G. Neto, A Design Proposal for Ultrawide Band Frequency
Selective Surface, Journal of Microwaves, Optoelectronics and Electromagnetics Applications, 12 (2), 2013, 2943-
2946.
[9] S. Baysakhia, R. Sivasamy, M. Kanagasabai and S. Periaswany, Novel Compact UWB Frequency Selective Surface
for Angular and Polarization Independent Operation, Progress In Electromagnetics Research Letters, 40, 2013, 71-
79.
[10] M. V. Narayana, I. Govardhani, A. Vikranth, S. K. Nizamuddin, C. Venkatesh, V. V. V. S. Krishna and K.
Rajkamal, Design of a Frequency selective surface with multiple four legged Slots, International Journal of Modern
Engineering Research (IJMER), 2012, 2 (3), 807-810.
[11] H. Oraizi and M. Afsahi, Analysis of Planar Dielectric Multilayers as FSS by Transmission Line Transfer Matrix
Method (TLTMM), Progress in Electromagnetics Research, 74, 2007, 217-240.
[12] T. L. Silva, A. L. P. S. Campos, A. G. dAssuno and R. H. C. Manioba, A comparative Study Two Numerical
Techniques to Analyze Double Screen Frequency Selective Surface, Microwave and Optical Technology Letters, 55
(9), 2013, 981-985.
[13] A. L. P. S. Campos and T. L. Silva, Spectral domain analysis of double screen frequency selective surfaces, Journal
of Microwaves, Optoelectronics and Electromagnetics Applications, 11(1), 2012, 81-92.
[14] M. Titaouine, N. Raveu, A. G. Neto and H. Baudrand, Dual-band and Enhanced Band FSS Characterization Using
WCIP Method, Microwave and Optical Technology Letters, 52 (4), 2010, 836-839.
[15] R. H. C. Manioba, A. G. dAssuno and A. L. P. S. Campos, Improving Stop-Band Properties of Frequency
Selective Surface with Koch Fractal Elements, 2010 International Workshop on Antenna Technology (iWAT),
Lisbon, 2010, 57-60.
[16] T. Cwik and R. Mittra, The cascade connection of planar periodic surfaces and lossy dielectric layers to form an
arbitrary periodic screen, IEEE Transactions on Antennas and Propagation, 35 (12), 1987, 1397 1405.
[17] C. Wan and J. A. Encinar, Efficient computation of generalized scattering matrix for analyzing multilayered
periodic structures, IEEE Transactions on Antennas and Propagation, 43 (11), 1995, 1233 1242.
[18] S. W. Lee, G. Zarrillo and C. L. Law, Simple Formulas for Transmission Through Periodic Metal Grids or Plates,
IEEE Transactions on Antennas and Propagation, 30 (5), 1982, 904 909.