Efficient Analysis of Complex FSS Structure Using the WCIP Method

A rigorous full wave technique based on the Transverse Wave Concept Iterative Procedure (WCIP) is used to design a complex Frequency Selective Surface (FSS). These surfaces include a periodically arrangement of identical circuit. There are used as filters and reflector antenna as well as deep-space exploration for multi-frequencies operations. A simple FSS structure is studied in first stage to validate our approach. In second stage two different complex structures are studied. The good agreement between simulated and published data justify the design procedure.


Introduction
Frequency Selective Surfaces (FSSs), which find widespread applications as filters for microwaves and optical signals, have been the subject of extensive studies in recent years [1][2][3][4].These surfaces include a periodically arranged metallic patch elements or aperture elements within a metallic screen and exhibit total reflection (patches) or Transmission (apertures) in the neighborhood of the element resonance [1].Their performances depend on the substrate characteristics, element type, dimensions and the spacing between elements.
The response parameters are predicted by analyzing the surface using different techniques [5][6][7].However, the small dimensions of the circuit produce some problems in result precisions.Thus, the coupling conditions between the different elements must be taken into account.Then, the efficiency of used method, their memory consumption and time requirement are usually made these methods unsuitable for optimization.
This paper presents the analysis of simple and complex passive FSS by the iterative method (WCIP).The WCIP technique takes the advantage of simplicity in its procedure based on Fast Modal Transform (FMT) in passage between spatial and spectral domain [8,9].In addition, there is no matrix inversion was required and the convergence was insured independently of the circuit complexity.Further, there is unlimited shapes of circuit are imposed [10].The simulation results are validated with those calculated with HFSS commercial code and recently published experimental results.

Theory: WCIP Formulation
The general Frequency Selective Surface structure is depicted in Figure 1.The circuit interface is constituted of two sub domains: metal and dielectric.It is deposed on homogeneous dielectric substrate with thickness h and permittivity r  .
WCIP method is based on the full wave transverse wave formulation and the on collection of information at the interfaces.A multiple reflection procedure is started using initial conditions and stopped once convergence which is achieved.Two related operators incidented waves and scattered waves in the spatial domain and in the spectral domain governs the iterative procedure.They are: the scatting operator S  and the reflection .
 We consider the printed circuit, the wave concept is  introduced to express the boundary conditions on the interface Ω (Figure 2).
The incident waves i and the scattered waves i A B are calculated from the tangential electric and magnetic fields i E and i H as: where i indicates the medium 1 or 2 corresponding to a given interface .0i  Z is the characteristic impedance of the same medium i and i J being the surface current density vector given as: with n being the outward vector normal to the interface.Thus, the tangential electric and magnetic fields can be calculated from: The scattered waves are related to the incident waves as: S  is a scattering operator defines in the spatial and it accounts for the boundary conditions.The scattered waves i B will be reflect to generate the incident waves for the next iteration but after adding the incident source waves : being the reflection operator and it is defined in the s rator pectral domain. 

Scattering Ope
Replacing ( 3) in ( 7) results in: The metal domain scattering operator M S is given in e terms of the metallic domain generator In the dielectric domain, the boundary conditions be satisfied on the interface are: In the lumped elements domain, the boundary to be verified is given by: ( ) Then, the total scatting operator S  is given as:

Reflection Operator Dete des where rmination
The modes are decoupled in the domain of mo each mode is characterized by its own reflection coefficient, the need to pass to spectral domain is necessary.
To enable this operation, a transform known as the Fast Modal Transform FMT defined and to go back to spatial domain, FMT  is will be used.The reflec efficient in the s tion co pectral domain is given by: where imn Y  is the admittance of the mn mode at the medium i and α stands for the modes TE or TM.
When no closing ends exist, imn Y  can be calculated by [5,6]: being the propagation constant of the medium i and mn y it is given by:

Fast Mo
The 1 FMT FMT  pair permits to go from spatial doto the spectral main domain and back to the spatial domain [7].It is summarized in the following two equations. where

Applications
our method, we consider at first the In order to validate structure of Figure 3.The dashed lines are the hypothetical periodic walls assumed to reduce the analysis to that of the repeated unit cell with dimensions of a = 3 mm, b = 3 mm, h 1 = 4 mm, h 2 = 3.125 mm, ε r2 = 2.6, ε r1 = 1.The TE01 mode is used as excitation and the iterative process is stopped after 200 iteration.The results of our method are compared to those calculated with HFSS and depicted in Figures 4 and 5 for the reflected and the transmitted respectively coefficient.In the two cases, a good agreement is obtained between results.For the s onator circu nd example is represented in Figures 8 and  10 econd application, we consider res its optimized by the genetic algorithm [8] included in WR90 waveguide (a = 22.86 mm, b = 10.16 mm) The TE01 mode is used as excitation and the iterative process is stopped after 200 iteration.This circuit is used as a filter.Two different configurations of the filter are used.Figure 6 show the first form of the filter.The transmission parameter is depicted in Figure 7.As shown in Figure 7, the WCIP results are in agreement with published data [11].
The seco and its response is depicted in Figures 9 and 11.In this figures, the WCIP results are compared to those of     ublished data [11].Then, we can conclude that the two , the Wave concept Iterative Method nc p results are agreed in pass band.

Conclusions
All over this paper has been outlined and used to study a complex Frequency Selective Surface.A simple FSS structure and two complex are studied.The comparison of our simulated data with those from commercial code and recent published data allowed the validation of the proposed method.
boundary conditions on the metal domain Two domains characterizing interface  of a lo ed FSS are: the dielectric domain and the m tal domain.They can be represented using Heaviside unit steps as:

2
medi are permittivity of the vacuum, the relati dal Transform FMT um i and the permeability of the vacuum respectively.

Figure 9 .
Figure 9. Transmission response as function of frequency.

Figure 11 .
Figure 11.Transmission response as function of freque y.