<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JEMAA</journal-id><journal-title-group><journal-title>Journal of Electromagnetic Analysis and Applications</journal-title></journal-title-group><issn pub-type="epub">1942-0730</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jemaa.2011.311071</article-id><article-id pub-id-type="publisher-id">JEMAA-8336</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Efficient Analysis of Complex FSS Structure Using the WCIP Method
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>assi</surname><given-names>Aroussi</given-names></name></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Latrach</surname><given-names>Latrach</given-names></name></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Noureddine</surname><given-names>Sboui</given-names></name><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>Gharsallah</given-names></name></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Abdelhafidh</surname><given-names>Gharbi</given-names></name></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Henry</surname><given-names>Baudrand</given-names></name></contrib></contrib-group><author-notes><corresp id="cor1">* E-mail:<email>noureddine.sboui@fst.rnu.tn(NS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>04</day><month>11</month><year>2011</year></pub-date><volume>03</volume><issue>11</issue><fpage>447</fpage><lpage>451</lpage><history><date date-type="received"><day>August</day>	<month>2nd,</month>	<year>2011</year></date><date date-type="rev-recd"><day>September</day>	<month>1st,</month>	<year>2011</year>	</date><date date-type="accepted"><day>September</day>	<month>25th,</month>	<year>2011.</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  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.
 
</p></abstract><kwd-group><kwd>FSS</kwd><kwd> WCIP</kwd><kwd> Wave</kwd><kwd> 2D-FFT Algorithm</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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-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 [<xref ref-type="bibr" rid="scirp.8336-ref1">1</xref>]. Their performances depend on the substrate characteristics, element type, dimensions and the spacing between elements.</p><p>The response parameters are predicted by analyzing the surface using different techniques [5-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.</p><p>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 [<xref ref-type="bibr" rid="scirp.8336-ref10">10</xref>]. The simulation results are validated with those calculated with HFSS commercial code and recently published experimental results.</p></sec><sec id="s2"><title>2. Theory: WCIP Formulation</title><p>The general Frequency Selective Surface structure is depicted in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The circuit interface is constituted of two sub domains: metal and dielectric. It is deposed on homogeneous dielectric substrate with thickness h and permittivity<img src="3-9801114\bb97d05e-3b13-4360-87cf-a872efdb732f.jpg" />.</p><p>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 <img src="3-9801114\fab400a1-0743-40ec-95f0-127c8412dcb0.jpg" /> and the reflection<img src="3-9801114\e358ae36-d238-46b1-89de-36ee9076ac23.jpg" />.</p><p>We consider the printed circuit, the wave concept is</p><p>introduced to express the boundary conditions on the interface Ω (<xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><p>The incident waves <img src="3-9801114\310b17cd-e23b-480d-b568-33b7c3aa45ef.jpg" /> and the scattered waves <img src="3-9801114\a3008fd9-2aa6-445f-8c9e-edac027e0746.jpg" /> are calculated from the tangential electric and magnetic fields <img src="3-9801114\536b3cfa-75f5-407f-9721-818ee76231ba.jpg" />and <img src="3-9801114\e6e5575a-9132-4d59-9429-0e7c2cb03964.jpg" />as:</p><p><img src="3-9801114\7586cee2-db78-459f-8d6e-d491d3e69232.jpg" /></p><disp-formula id="scirp.8336-formula88772"><label>(1)</label><graphic position="anchor" xlink:href="3-9801114\e2f5a4bd-5d17-4e2f-bd19-95aea54038ca.jpg"  xlink:type="simple"/></disp-formula><p>where i indicates the medium 1 or 2 corresponding to a given interface<img src="3-9801114\4c926894-a74d-426f-8ddb-e15b1ef8e516.jpg" />. <img src="3-9801114\60344917-2332-49cc-99bb-89d29bd2c56c.jpg" />is the characteristic impedance of the same medium i and <img src="3-9801114\c47c61d3-0078-4ea5-bd74-0b67d257e49f.jpg" />being the surface current density vector given as:</p><disp-formula id="scirp.8336-formula88773"><label>(2)</label><graphic position="anchor" xlink:href="3-9801114\121f0fcc-0a8b-4ac7-90be-699bd2b43277.jpg"  xlink:type="simple"/></disp-formula><p>with n being the outward vector normal to the interface. Thus, the tangential electric and magnetic fields can be calculated from:</p><disp-formula id="scirp.8336-formula88774"><label>(3.1)</label><graphic position="anchor" xlink:href="3-9801114\8e48ca76-3327-413d-aa79-1bbafe213704.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.8336-formula88775"><label>(3.2)</label><graphic position="anchor" xlink:href="3-9801114\7c976df5-ff9d-4bce-abfc-7a05877719a2.jpg"  xlink:type="simple"/></disp-formula><p>The scattered waves are related to the incident waves as:</p><disp-formula id="scirp.8336-formula88776"><label>(4)</label><graphic position="anchor" xlink:href="3-9801114\a85132cc-15b0-44f9-949f-21cb88fe892a.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-9801114\cdbdceb3-998e-4135-abb4-a0b55efadc80.jpg" />is a scattering operator defines in the spatial and it accounts for the boundary conditions. The scattered waves <img src="3-9801114\d3de6dd9-9cf4-4455-b76f-99b6ab4f82da.jpg" /> will be reflect to generate the incident waves for the next iteration but after adding the incident source waves<img src="3-9801114\062a20c7-78df-45f5-b943-1f327cbe6290.jpg" />:</p><disp-formula id="scirp.8336-formula88777"><label>(5)</label><graphic position="anchor" xlink:href="3-9801114\6f372a8c-8fe6-4323-a584-feada9a705fb.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-9801114\b16af42d-52d6-4431-b952-ab7fe5611c24.jpg" />being the reflection operator and it is defined in the spectral domain.</p><sec id="s2_1"><title>2.1. Scattering Operator <img src="3-9801114\73d47a4a-b11f-47c1-8bfd-5cbe239ea0a9.jpg" /> Determination</title><p>Two domains characterizing the interface <img src="3-9801114\5741cacf-a96e-43eb-8ecc-cbf8a234fa5a.jpg" /> of a loaded FSS are: the dielectric domain and the metal domain. They can be represented using Heaviside unit steps as:</p><disp-formula id="scirp.8336-formula88778"><label>(6.1)</label><graphic position="anchor" xlink:href="3-9801114\491bf7df-ca64-425f-99e1-57cf72c36526.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.8336-formula88779"><label>(6.2)</label><graphic position="anchor" xlink:href="3-9801114\3da70e75-23c8-4baa-96a7-a57b8475df80.jpg"  xlink:type="simple"/></disp-formula><p>The boundary conditions on the metal domain <img src="3-9801114\80e9e191-5ecf-4653-8c59-3d92c7ae4028.jpg" /> are:</p><disp-formula id="scirp.8336-formula88780"><label>(7)</label><graphic position="anchor" xlink:href="3-9801114\a0f7b73b-21d2-456d-9a83-66da3bbb947b.jpg"  xlink:type="simple"/></disp-formula><p>Replacing (3) in (7) results in:</p><disp-formula id="scirp.8336-formula88781"><label>(8)</label><graphic position="anchor" xlink:href="3-9801114\a258122d-732e-4e95-be87-621dd0097fc4.jpg"  xlink:type="simple"/></disp-formula><p>The metal domain scattering operator <img src="3-9801114\67366929-50bc-475b-96d0-bf0656a88521.jpg" />is given in the terms of the metallic domain generator <img src="3-9801114\f3fdb3af-3c67-40ca-9b0e-7f42b43983c8.jpg" /> by:</p><disp-formula id="scirp.8336-formula88782"><label>(9)</label><graphic position="anchor" xlink:href="3-9801114\a2da20fc-843c-48c3-8622-ab3cd54ee870.jpg"  xlink:type="simple"/></disp-formula><p>In the dielectric domain, the boundary conditions be satisfied on the interface are:</p><disp-formula id="scirp.8336-formula88783"><label>(10)</label><graphic position="anchor" xlink:href="3-9801114\c7788591-a891-48f5-9216-88e26ee078e4.jpg"  xlink:type="simple"/></disp-formula><p>Using (3) and (10), and defining<img src="3-9801114\41092055-6f26-4d2e-9fb9-b1e8114f1dee.jpg" />. Thus, the dielectric domain scattering operator <img src="3-9801114\08c66a5a-902d-4088-b639-d218593cf600.jpg" /> can be given terms of the dielectric generator <img src="3-9801114\ae647ac7-d2f7-4aef-868a-b1c490e28220.jpg" /> as:</p><disp-formula id="scirp.8336-formula88784"><label>(11)</label><graphic position="anchor" xlink:href="3-9801114\4cb9467d-072b-44d5-9ae6-020de7980aa9.jpg"  xlink:type="simple"/></disp-formula><p>In the lumped elements domain, the boundary to be verified is given by:</p><disp-formula id="scirp.8336-formula88785"><label>(12)</label><graphic position="anchor" xlink:href="3-9801114\e4834e12-c96b-4122-b9f9-634bed10f0b1.jpg"  xlink:type="simple"/></disp-formula><p>Then, the total scatting operator <img src="3-9801114\f25f2e76-a496-4b47-9063-1fe13add6f79.jpg" /> is given as:</p><disp-formula id="scirp.8336-formula88786"><label>(13)</label><graphic position="anchor" xlink:href="3-9801114\5b97cf31-5370-4525-87ad-d568649251b4.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2"><title>2.2. Reflection Operator Determination</title><p>The modes are decoupled in the domain of modes where 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 <img src="3-9801114\82ff49a4-6783-444c-8609-20499dfc2193.jpg" /> defined and to go back to spatial domain, <img src="3-9801114\a783d958-8d53-43d8-9599-65a855601f09.jpg" />is will be used.</p><p>The reflection coefficient in the spectral domain is given by:</p><disp-formula id="scirp.8336-formula88787"><label>(14)</label><graphic position="anchor" xlink:href="3-9801114\aa7a6567-3a4e-4059-a6cd-6a705360ec04.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="3-9801114\d6e4c5e8-e322-4072-a3d4-6f946ee3f3fb.jpg" /> is the admittance of the mn mode at the medium i and α stands for the modes TE or TM.</p><p>When no closing ends exist, <img src="3-9801114\0349ce66-b644-4ada-ac9f-d572ca137422.jpg" />can be calculated by [5,6]:</p><disp-formula id="scirp.8336-formula88788"><label>(15)</label><graphic position="anchor" xlink:href="3-9801114\91e57439-3a27-4262-896c-38c17eada873.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.8336-formula88789"><label>(16)</label><graphic position="anchor" xlink:href="3-9801114\a8797da9-195b-435e-9312-3134753a4320.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-9801114\0215b87b-ec90-481a-8d39-3224a021a9e4.jpg" />being the propagation constant of the medium i and it is given by:</p><disp-formula id="scirp.8336-formula88790"><label>(17)</label><graphic position="anchor" xlink:href="3-9801114\fe2232ab-9f6f-420d-a904-616ca862e541.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-9801114\0649e68f-826b-4ad9-bcfc-e037ce04458e.jpg" />are permittivity of the vacuum, the relative of the medium i and the permeability of the vacuum respectively.</p></sec><sec id="s2_3"><title>2.3. Fast Modal Transform FMT</title><p>The <img src="3-9801114\d66e5c62-cf1c-407f-ba7f-8981cd6ccb73.jpg" /> pair permits to go from spatial domain to the spectral domain and back to the spatial domain [<xref ref-type="bibr" rid="scirp.8336-ref7">7</xref>]. It is summarized in the following two equations.</p><disp-formula id="scirp.8336-formula88791"><label>(18)</label><graphic position="anchor" xlink:href="3-9801114\e7473e03-c3a3-477f-b62a-724f9dea4d26.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.8336-formula88792"><label>(19)</label><graphic position="anchor" xlink:href="3-9801114\6cbbc937-eb2c-42f9-8b1a-d9e499daa34f.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="3-9801114\609de078-7ae5-4cef-866f-f02623cede2d.jpg" />and<img src="3-9801114\2b9f45a0-13fb-4489-89df-32e0c50980c5.jpg" />.</p></sec></sec><sec id="s3"><title>3. Applications</title><p>In order to validate our method, we consider at first the structure of <xref ref-type="fig" rid="fig3">Figure 3</xref>. 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<sub>1 </sub>= 4 mm, h<sub>2 </sub>= 3.125 mm, ε<sub>r</sub><sub>2 </sub>= 2.6, ε<sub>r</sub><sub>1</sub> = 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.</p><p>For the second application, we consider resonator circuits optimized by the genetic algorithm [<xref ref-type="bibr" rid="scirp.8336-ref8">8</xref>] 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. <xref ref-type="fig" rid="fig6">Figure 6</xref> show the first form of the filter. The transmission parameter is depicted in <xref ref-type="fig" rid="fig7">Figure 7</xref>. As shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>, the WCIP results are in agreement with published data [<xref ref-type="bibr" rid="scirp.8336-ref11">11</xref>].</p><p>The second example is represented in Figures 8 and 10 and its response is depicted in Figures 9 and 11. In this figures, the WCIP results are compared to those of</p><p>published data [<xref ref-type="bibr" rid="scirp.8336-ref11">11</xref>]. Then, we can conclude that the two results are agreed in pass band.</p></sec><sec id="s4"><title>4. Conclusions</title><p>All over this paper, the Wave concept Iterative Method 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.</p></sec><sec id="s5"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.8336-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">R. Mitter, C. H. Chan and T. Cwik, “Technique for An- alyzing Frequencey Selective Surfaces,” Proceedings of the IEEE, Vol. 76, No. 12, 1988, pp. 1593-1615.</mixed-citation></ref><ref id="scirp.8336-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">B. A. Munk, “Frequency Selective Surfaces: Theory and Design,” Wiley, New York, 2000.  
doi:10.1002/0471723770</mixed-citation></ref><ref id="scirp.8336-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">D. J. D. H. Werver, A. Monorchio, L. Lanuzza and M. J. Wilhelm, “The Design Synthesis of Multi-Band Artificial Magnetic Conductors Using High Impedance Frequency Selective Surfaces,” IEEE Transactions on Antennas and Propagation, Vol. 53, No. 1, 2005, pp. 8-17.  
doi:10.1109/TAP.2004.840540</mixed-citation></ref><ref id="scirp.8336-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">D. Hokim and J. I. Choi, “Design of a Multi-Band Frequency Selective Surface,” ETRI Journal, Vol. 28, No. 4, 2006, pp. 506-508.</mixed-citation></ref><ref id="scirp.8336-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">N. Bliznyuk and N. Engheta, “Numerical Study of Polarization Dependant Focusing for a Bi-Layer Planar FSS Reflective at Millimeter Wave Length,” Microwave and Optical Technology Letters, Vol. 40, No. 5, 2004, pp. 361-365.</mixed-citation></ref><ref id="scirp.8336-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">A. K. Bhattacharyya, “A Numerical Model for Multilayered Microstrip Phased-Array Antennas,” IEEE Transactions on Antennas and Propagation, Vol. AP-44, No. 10, 199, pp. 1386-13936. doi:10.1109/8.933485</mixed-citation></ref><ref id="scirp.8336-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">M. Bozzi, L. Perregrini, J. Weinzierl and C. Winnewisser, “Efficient Analysis of Quasi-Optical Filters by a Hybrid MoM/BI-RME Method,” IEEE Transactions on Antennas and Propagation, Vol. 49, No. 7, 2001, pp. 1054-1064.</mixed-citation></ref><ref id="scirp.8336-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">L. Latrach, N. Sboui, A. Gharsallah, A. Gharbi and H. Baudrand, “A Design and Modelling of Microwave Active Screen Using a Combination of the Rectangular and Periodic Waveguide Modes,” Journal of Electromagnetic Waves and Applications, Vol. 23, No. 11-12, 2009, pp. 1639-1648.</mixed-citation></ref><ref id="scirp.8336-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">N. Sboui, L. Latrach, A. Gharsallah, H. Baudrand and A. Gharbi, “A 2D Design and Modeling of Microstrip Structures on Inhomogeneous Substrate,” RF and Microwave Computer-aidedvEngineering, Vol. 19, No. 3, 2009, pp. 346-353. doi:10.1002/mmce.20354</mixed-citation></ref><ref id="scirp.8336-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Sboui, N., A. Gharsallah, A. Gharbi and H. Baudrand, “Global Modelling of Microwave Active Circuits by an Efficient Iterative Procedure,” International Journal of RF and Microwave Computer-Aided Engineering, Vol. 148, No. 3, 2001, pp. 209-212.</mixed-citation></ref><ref id="scirp.8336-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">M. Ohira, H. Deguchi, M. Tsuji and H. Shigesawa, “Novel Waveguide Filters with Multiple Attenuation Poles Using Frequency Selective Surfaces,” IEEE Transactions on Microwave Theory and Techniques, Vol. 53, No. 11, 2005, pp. 3320-3326.</mixed-citation></ref></ref-list></back></article>